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Abstract 



The work to be described falls into three parts: (1) the design, construction and testing 
of a continuous-wave (CW) microwave NMR spectrometer; (2) an NMR study of the 
hyperfine splittings of holmium trifluoride, supplemented by magnetometry; and (3) 
theoretical analysis. 

(1) The computer-controlled CW spectrometer was designed to supplement the Man- 
chester pulsed microwave spectrometer in situations where rapid nuclear relaxation 
makes spin-echo spectroscopy difficult. Its operating range is 4-8 GHz. Resonator 
designs and modulation strategies will be discussed in the light of practical experience. 

(2) Both CW and pulsed NMR have been used to study the field dependence of the 
hyperfine splittings of 165 Ho in H0F3 and, as a dilute substituent, in YF3. The low 
site symmetry results in a singlet crystal-field ground state for the Ho 3+ ion, giving 
Van Vleck paramagnetism and enhanced nuclear magnetism at low temperatures. The 
measurements were made at temperatures in the range 1.5 to 4.2 K and in fields of up 
to 8 T. This work has revealed, for the first time, distinct spectra from the two subtly 
inequivalent rare-earth sites in the orthorhombic unit cell. Because of the non-colinear 
spin structure of H0F3, the NMR and magnetometry measurements give independent 
and complimentary information about the ionic moments. 

(3) The measured hyperfine splittings have been interpreted in terms of a 15-parameter 
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crystal-field Hamiltonian appropriate to the C\h site symmetry. This work has entailed 
a substantial effort to clarify the notational confusion that exists in the literature. A 
computer program has been developed to automate conversion between notational con- 
ventions prior to diagonalization of the 136-dimensional electronic-nuclear Hamiltonian 
comprising the Zeeman, crystal-field and hyperfine interactions. 

Our experimental results are in fair agreement with calculated magnetizations and hy- 
perfine splittings based on crystal-field parameters derived from high-resolution optical 
spectroscopy. However, there are discrepancies which suggest the need for refinement 
of the crystal-field parameters. The existing parameters predict that the hyperfine 
splitting will vary dramatically with the orientation of the applied field but technical 
difficulties have prevented an investigation of this effect during the time scale of the 
present work. 
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1. Introduction 



The rare-earth elements and compounds containing them exhibit a wide variety of 
magnetic properties. Moreover, they have proved amenable to detailed study. Since 
the 1950s, when rare earths of reasonable purity became available, they have been 
intensively studied revealing not only many physically interesting phenomena, but also 
industrially useful materials. In this work we shall be concerned solely with insulating 
rare-earth compounds. 

Table [0| shows the rare-earth or lanthanide group of the periodic table. The elements 
La to Eu are often called light rare earths and the elements Gd to Lu heavy rare earths. 
Yttrium is not a rare earth but is often given honorary rare earth status because it 
is chemically very similar. Most rare earths are triply ionised in solids; in this state 
only the 4/ shell is partially filled. The radius of the 4/ shell is several times smaller 
than typical interionic separations and the 4/ electrons do not take part in chemical 
bonding. However, the 4/ electrons have large angular momentum and dominate the 
magnetic properties of the ion. 



57 


58 


59 


60 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


71 


La 


Ce 


Pr 


Nd 


Pm 


Sm 


Eu 


Gd 


Tb 


Dy 


Ho 


Er 


Tm 


Yb 


Lu 



Table 1.1: The rare-earth elements. 
The 'outer electrons' shield the 4/ shell from its surroundings and the rare earths 
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retain their free-ion character in solids to a high degree. As in the free ions, the L-S 
or Russell- Saunders coupling scheme is a good approximation in solids. In general L, 
S and J are fairly good quantum numbers, and spin-orbit coupling gives highest J as 
ground state for heavy rare earths. 



1.1 Rare-earth magnetism 

Rare-earth ions maintain their free-ion character in solids more closely then any other 
elements. To a good first approximation the crystal-field and Zeeman interactions 
leave the spin-orbit coupling intact; their principal effect is to lift the degeneracy of the 
J-manifolds. 

Bethe's seminal paper Q showed that open-shell energy levels of an ion in a crystalline 
environment were associated with the symmetry of the site. He also considered the 
terms required to describe the interaction of the ion with its environment for differ- 
ent point group symmetries. This was the birth of crystal-field theory. Condon and 
Shortley p5[ provided the basic techniques for a perturbation approach to crystal-field 



calculations. Then, in a classic series of papers, Racah |30|, |61|, |62| developed the pow- 
erful tensor operator notation and discussed calculation of the matrix elements. All 
of this work was based in group theory and on the work of Wigner, especially the 
Wigner-Eckart theorem 

From the original well-structured works of Bethe, Condon and Shortly, and Racah, 
crystal-field theory has been developed in an ad hoc manner. Abragam and Bleaney O 
note that it is unfortunate that the pioneer work of Stevens [74] and of Elliott and 



Stevens 33, 35, 34 1 was not expressed in the more rational formalism of Racah. How- 
ever, Elliott and Stevens were successful in accounting, in detail, for paramagnetic 
resonance data on some rare-earth salts. In appendix A we review parameter con- 
ventions for crystal fields and the inter-relationships between them. Having identified 
the more coherent conventions, appendix B describes the transformation of parame- 
ters under coordinate rotation. This is useful in helping to relate the crystal fields of 
inequivalent ions in sites with surroundings related by rotation. 
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Ab initio calculation of crystal-field and free-ion parameters remains intractable, so 
the parameters must be determined experimentally. Data from NMR spectroscopy 
alone are not usually sufficient to determine crystal- field parameters; the parameters 
are usually obtained from optical or neutron spectroscopy. Crystal-field parameters 
so obtained provide a basis for the detailed analysis of NMR data. For example Car- 
boni |2(]] used the computed ground state of Ho 3+ in Ho(OH)3 to include J-mixing in 
his analysis of NMR data and to refine the hyperfine parameters of Ho 3+ . 



1.2 Rare-earth trifluorides 

The rare earths form many fluorides, the structures and chemisty of which are reviewed 
by Greis and Haschke [pq j. Industrially, the fluorides are used for the manufacture of 
arc carbons with well balanced light emission. Of these the anhydrous trifluorides have 
been most intensively studied. They are chemically stable in air and moisture at room 
temperature. As a result, rare-earth trifluorides have proved a useful intermediate in 
the preparation of high-purity rare-earth metals [f7q1 - Pure metals are produced from 
the trifluorides by reduction with calcium metal in an inert atmosphere. 

In this work we consider only the trifluorides that are isostructural with YF3: in 
particular H0F3 and YF3. There have been several X-ray and neutron structure deter- 
minations for rare-earth trifluorides which are in fairly good agreement. There is no 
argument about the C\h site symmetry of the rare-earth ions. Using just the structural 
and symmetry data we deduce that there are two inequivalent rare-earth sites in the 
unit cell (chapter 5). 

H0F3 is a Van Vleck paramagnet with strongly enhanced nuclear magnetism. It orders 
antiferromagnetically at Tjv = 0.53 K pl[ , due mainly to the dipole-dipole interaction. 
The ordered state has been studied by neutron diffraction |13T]. Of the other trifluo- 



rides, only TbF3 has been studied in the ordered state [41]. Like H0F3, TbF3 orders 
antiferromagnetically, but at a considerably higher temperature: Tjv = 3.95 K EI ], 
In this work we have used NMR to study the field dependence of the hyperfine split- 
tings of Ho 3+ in paramagnetic H0F3, and as a dilute substituent in YF3. We have 
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also measured the magnetization in fields of up to 10 T along the principal crystal- 
lographic directions. The experimental data are compared with calculations based on 
crystal-field parameters derived from optical spectroscopy. 



1.3 CW NMR 

Rare-earth NMR is usually performed by pulsed techinques, principally spin-echo. 
Whilst pulsed NMR techniques have become extremely sophisticated for chemical sys- 
tems, NMR of rare-earths is constrained by the technical difficulties associated with 
the very high frequencies (up to 7 GHz) and fast relaxation. As part of this work I have 
built a new continuous-wave (CW) microwave spectrometer. Chapter 3 discusses the 
CW technique and compares it with pulsed techniques. Various detection strategies and 
sample-cell designs are considered. Chapter 4 describes the new spectrometer and how 
it has been used. The majority of the NMR data in this thesis were taken using the CW 



spectrometer; the rest were taken using the Manchester pulsed spectrometer [21, 49 1. 



The CW spectometer operates over the range 4-8 GHz, which has allowed study of 
the hyperfine interection of 165 Ho in H0F3. It can also be used for studies of 141 Pr. 
The microwave system is very simple and could fairly easily be modified to extend the 
operating range down to 2 GHz, giving access to the resonances of 159 Tb, 169 Tm and 
m Yb. 
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2. Theory 



In order to understand the properties of rare-earth compounds we need to consider 
both the interactions within the rare-earth ion and the interactions of the ion with 
its surroundings. In atomic physics, interactions are traditionally arranged in order 
of descending strength and each is treated as a perturbation on the previous one. 
The perturbation approach has proved remarkably successful for rare-earth ions. The 
ground LS-terms of rare-earth ions are given, to a good first approximation, by Hund's 
rule. In the heavy rare earths the spin-orbit coupling results in the ground J-manifold 
having maximum J. Whilst we adopt the perturbation approach to arrive at the ground 
J-manifold, the crystal-field, Zeeman and hyperfine interactions are treated together. 
This is unavoidable for the the electronic crystal-field and Zeeman interactions which 
may be of similar strengths. The hyperfine interaction is weak compared with the 
electronic interactions but cannot be treated as a perturbation when the electronic 



levels are degenerate, or nearly degenerate: see section 2.1 



The hyperfine interaction is of particular importance to this work as NMR is our 
principal experimental technique. Section |2.6| discusses hyperfine spectra on the basis 
of an effective nuclear Hamiltonian; typical parameter values are given. Discussion of 
NMR techniques is deferred to chapter 3. In this work we consider only insulators, for 
which the exchange interaction is very weak. The theory described is applicable to the 
majority of the heavy rare-earth ions, but we shall concentrate on the Ho 3+ ion. 
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Data from bulk magnetization and low-temperature heat capacity measurements can 
be related to the ionic energy levels and eigenstates. Low-temperature heat capacity 
measurements are easiest to interpret for systems where we can use the approximation 
of a two-level system, i.e. where a low-lying first excited state is well separated from the 



next excited state. In section 2.7 we consider the Schottky anomaly in the two-level 



approximation. Section 2.5 relates the magnetization to the thermal average of the 
electronic angular momentum. 



2.1 The Hamiltonian 

The Hamiltonian for a rare-earth ion in a crystalline environment consists of free-ion, 
crystal-field, Zeeman and hyperfine terms. We treat the didole-dipole and exchange 
interactions in the molecular-field approximation by including an effective molecular 
field in the Zeeman term. In sections |2.1.1| and 2.1.2 we briefly describe the contribu- 



tions to the Hamiltonian which determine the ground J-manifold in the L-S coupling 
approximation . 

We may write the total Hamiltonian as: 

H = Ho + Hi + H so + H e i + Hhf, (2.1) 

where Ho is the central-field approximation to the Coulomb interaction; Hi is the non- 
central part of the Coulomb interaction; H so is the spin-orbit interaction; H e i is the 
interaction of the electrons with their surroundings, both magnetic and electric; and 
Hhf is the hyperfine interaction. 

For rare-earth ions it is usually a good approximation that 

Ho » Hi > Hso > H e i. (2.2) 

We shall consider the terms in sequence; each term as a perturbation on the previous 
one. The term H e \ includes both the crystal-field and Zeeman interactions: H e i = 
H c f + H z . At high fields the strength of the Zeeman interaction can approach that 
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of the crystal field, so they cannot be treated separately. The hyperfine interaction is 
much weaker than the crystal-field and Zeeman interactions and can usually be treated 
as a perturbation (H e i 3> TChf)- However, this approximation can break down in the 
region of electronic level crossings. This situation has been investigated by Han j5| . 

With current computing power it is not necessary to treat the hyperfine interaction as 
a perturbation on H e i\ therefore we diagonalise H e \ and Hhf together: see appendix C. 
The combined electronic- nuclear Hamiltonian is (21 + 1)(2J + l)-dimensional. 

2.1.1 Coulomb Hamiltonian, Ho + Hi 

The dominant contribution to the Coulomb Hamiltonian may be reprsented by a 
central-field Hamiltonian Ho which determines the arrangements of electons in shells, 
the configurations. In the rare earths, all shells apart from the 4/ are completely full or 
completely empty and so have no net angular momentum or magnetic moment. Omit- 
ting the labels of closed shells, the ground configuration of the holmium ion Ho 3+ is 
4/ 10 . The first excited configuration of holmium is > 50000 K above the ground con- 
figuration. There can be admixture of different configurations, usually referred to as 
the configuration interaction |34| , but this effect is very small and we shall not consider 
it further. In this work we assume the 4/ 10 ground configuration of holmium. 

The Coulomb repulsion Hi couples the individual electronic angular and orbital mo- 
menta such that L and S are good quantum numbers. This splits the configuration 
into non-degenerate L5-terms, a situation known as L-S or Russell- Saunders coupling. 

2.1.2 Spin-orbit coupling, H so 

The principal effect of the spin-orbit interaction is to couple L and S to form a total 
angular momentum J, with J as a good quantum number. Here we are assuming L-S 
coupling which is a good approximation in rare-earths because Hi 3> H so - Thus the 
spin-orbit interaction may be written as 

Hso = (L.S, (2.3) 
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where £ is the coupling parameter. In Ho 3+ the first excited manifold ^I-j is ~ 7500 K 
above the 5 Is ground manifold [29|. The spin-orbit coupling shifts the energies of the 



LS-terms slightly, and also mixes different L5-terms: see section 2.1.3 



In this work we use the L-S states as the basis for electronic calculations. They are 
labelled | tLSJMj), where r is and extra index to identify different terms with the 
same L and S. We use Mj to denote the eigenvalues of J z . In heavy rare-earth ions, 
the spin-orbit interaction is such that L and S couple to give the maximum possible 
J = L + S, and conversely for light rare-earths. The ground state of Ho 3+ is 5 /s, i.e. 
L = 6, 5 = 2, J = 8 and Mj = —8 ... 8. 

2.1.3 Intermediate coupling 

Spin-orbit coupling admixes states of different L and S but with the same J. The 
effect of this is that J is a much better quantum number than either L or S. Although 
the L5-terms are well separated, the admixture is a significant effect and is usually 
allowed for by the use of modified operator-equivalent coefficients (J || atk \\ J) and 
Lande g-factors g r The values of (J \\ \\ J) for the ground manifold of Ho 3+ 
are modified by ~8%, and this is known as intermediate coupling. Estimates for the 



operator-equivalent coefficients for some of the rare-earths are given by Dieke 2S] and 
for Ho 3+ by Rajnak and Krupke [)63[|. 



Carboni |2(| discusses the admixture of L5-terms by the spin-orbit coupling. If we 
write the ground state | £o) as linear combination of L-S states (our chosen basis) then 

| Co) =]>> I LiSiJi), (2.4) 

i 

where the summation is over all LSJ-manifolds i. For Ho 3+ in Y(OH)3 and in LaCl3 
the ground state is 

| Co) ~ 0.967 5 I + 0.116 3 A"i -0.221 3 K 2 - 0.031 3 L - 0.006 3 M 

+0.011 ^1 + 0.033 l L 2 . (2.5) 

For both compounds the Lande g-factor is g r = 1.2417 with an uncertainty estimated 
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at ±0.00005 @. This is about 0.7% smaller than the pure L-S value of 1.25. 



2.1.4 J-mixing 

In the free ion, J is a rigorously good quantum number if we ignore the hyperfine 
interaction. In general, the crystal-field and Zeeman interactions have matrix ele- 
ments between different J-manifolds so J ceases to be a good quantum number, an 
effect known as J-mixing. The splitting between the lower manifolds in heavy rare 
earths is much larger than in light rare earths, so J-mixing is less significant in the 
heavy rare earths. In holmium the ground and first excited manifolds are separated 



by ~ 7500 K [32]. With such a large separation we expect J to remain a fairly good 
quantum number. In high precision studies J-mixing can be significant, see for example 
Carboni [2C]. J-mixing has not been included in this work. 



2.2 Crystal-field Hamiltonian, 7i c f 

The crystal-field interaction is the interaction of the aspherical electronic charge dis- 
tribution of the ion with the inhomogeneous electric field produced by the surrounding 
ions. Firstly, we make the approximation that the charges producing the crystal field 
do not overlap with the 4/ electrons. Then Laplace's equation, V 2 V = 0, will hold 
for the potential experienced by the 4/ electrons and we may expand the crystal-field 
interaction in terms of spherical harmonics, Yk q (6, tp): 

n 

W C / = EE(^) fc ^W^). (2-6) 

i=l kq 

where (rj, 9i,(fi) is the position of the i th electron and the Al are the coefficients of the 
expansion. The number of parameters in the expansion is limited by the requirements 
that k must be even, and k < 21. For / electrons / = 3, so at most 27 parameters are 
required: we need only the terms for k = 2, 4, 6 and q = —k . . . k. 

In writing the expansion we have assumed a coordinate system (r, 9, (p) relating the 
ion to its surroundings. This choice of coordinate system will affect the parameters 
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A q k of the expansion and, possibly, the number required. The number of parameters 
required depends on the site symmetry in terms of the chosen coordinate system. This 
point is discussed in appendix A where the minimum number of parameters required 
for different site symmetries are given. At this stage we take the opportunity to define 
the crystal-field axes, shown in figure [O]. Appendix B describes the transformation of 
crystal-field parameters corresponding to rotation of the coordinate system. 



Z 



Figure 2.1: Definition of crystal- field axes. 



Application of the Wigner-Eckart theorem allows equation |2.6| to be recast in terms 
of a set of angular momentum operators. Adopting the convention of Morrison and 
Leavitt | |52|] we may write the crystal-field Hamiltonian as 

n cf = £ £ JU*Cf, (2.7) 

fc=2,4,6g=-fc...fc 

where the M k are parameters and the C q are tensor operators which may be written 
in terms of the components of J. The Mq are real and the M q are complex for q ^ 0, 
also M k = (-l) q M k *. The operators C k are given by C* = (J \\ at \\ J)T k , where 



the T k are listed in table [Aq and the (J \\ at \\ J) are operator- equivalent coefficients. 



Values of (J \\ \\ J) for the ground manifold of Ho 3+ are given in table 2J There 



are many other conventions for parametrizing the crystal field: see appendix A. 
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The principal effect of the crystal-field interaction is to lift some or all of the degeneracy 
in Mj. It can also cause some J-mixing (see section |2.1.4| above). The crystal-field 
model represents the independent interaction of equivalent electrons with an effective 
potential. In using it to fit experimental data, any other effects that may be expressed 
in this form will automatically be included. The principal component of the interaction 
comes from the electrostatic potential due to surrounding ions. However, other mech- 
anisms can make significant contributions to an interaction of this form, eg., covalency 



and the configuration interaction [54]. 



2.3 Zeeman interaction 

The Zeeman interaction is the interaction of the electronic moment with a magnetic 
field B. The Hamiltonian is 

U z = g jf J. B B.J (2.8) 

where g d is the Lande g-factor. It is convenient to include all effective fields in B and 
we write: 

B = B app + Bdip + B m , (2.9) 



where B app is the applied field; B^p is the dipolar field (see section 2.5); and B m is 
the molecular field used to model the exchange interaction. 



Following McCausland and Mackenzie ]50|] the exchange field is considered to interact 
with the projected spin a = {g } — 1) J. The projected spin is related to J by a scaling 
factor so the interaction can be rewritten in terms of an equivalent molecular field: 

B m = X, (2.10) 

where X is the exchange field. It is hard to calculate the exchange field and B m is 
usually treated as a parameter. The exchange interaction is much weaker in insulating 
materials than in metals where the conduction electrons interact with the ionic spins. 
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The exchange field is considered by Bunbury et al |18|] who, using the Scott's data [fi 
conclude that B m « — S^ p /4 in Ho(OH)3. 

2.4 Hyperfine interaction 



In this section we follow the notation of McCausland and Mackenzie [5C]. The hyperfine 
interaction may be split into dipolar and quadrupolar terms: 

H-hf = Kfifd + Khfq, (2-11) 

where Tihfd is the magnetic dipole Hamiltonian and T~Chfq is the electric quadrupole 
Hamiltonian. There are higher order interactions but their effect is not significant in 
our measurements so they are neglected. Both the dipole and quadrupole terms have 
intra-ionic and extra-ionic components, and it is convenient to separate them: 

Whf = ti'hfd + K'hfd + H'hfq + K'hfq> (2.12) 

where the intra-ionic terms are denoted by a single prime and extra-ionic terms by a 
double prime. 

2.4.1 Dipolar hyperfine interaction 

The intra-ionic part of the dipolar hyperfine interaction is given by 

H' hfd = hA(J.I), (2.13) 

where A is the dipolar coupling coefficient expressed as a frequency, hence the inclusion 
of Planck's constant h. For holmium A has been deduced from EPR experiments by 



Bleaney [1C] and refined by Carboni [20|. The dipolar and other hyperfine parameters 
for the rare-earth are tabulated by Han pC| . When the hyperfine interaction is con- 
sidered as a perturbation on the electronic Hamiltonian it is customary to express the 
intra-ionic dipolar parameter as a' = A J (see section ^1]). 
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The extra-ionic dipole interaction is the Zeeman interaction of the nuclear dipole mo- 
ment with the magnetic field at the ion. It is thus analogous to equation |2.8| : 

K fd = g im B".I, (2.14) 

where <j is the Lande g- factor for the nucleus (see table |2.1| ). The magnetic field B" 
may be written as 

B" = B app + B dip (2-15) 

where B app and B dip are the applied and dipolar contributions respectively. In metals 
there is also a contribution from the conduction electrons; here we are concerned only 
with insulators. 



2.4.2 Quadrupolar hyperfine interaction 

The intra-ionic part of the quadrupolar hyperfine interaction is 

2 



K'hfq — hC 



3(j./) 2 + |(j./)- j(j + + 1; 



(2.16) 

where C is the quadrupole coupling coefficient in the notation of Bunbury at al IS], 



expressed as a frequency. When the hyperfine interaction is considered as a perturba- 
tion on the electronic Hamiltonian it is customary to write the intra-ionic quadrupolar 
parameter as Pq = CJ(2J — 1) (see section |2.6| ). 



The extra-ionic quadrupolar interaction is the interaction of the nuclear quadrupole 
moment with the electric-field gradient at the nucleus. In this way it is analogous 
to the quadrupolar crystal-field interaction of the electrons. In the absence of di- 
rect information, the electric field gradient (EFG) at the nucleus can be estimated 
from the crystal-field parameters. First, the EFG 'seen' by the ion can be estimated 
from the crystal-field parameters using the electronic antishielding factor and the elec- 
tronic quadrupole moment. Then, 'crystal-field parameters' for the nucleus can be 
estimated from the extra-ionic EFG using the nuclear antishielding factor and the nu- 



clear quadrupole moment. Details of this procedure are given by Bunbury et al [18] for 
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the Bq term (using the crystal-field notation of Baker, Bleaney and Hayes Han 4C] 
added B\ to cope with orthorhombic symmetry. 

Here, we present a more general approach which includes all quadrupolar terms so 
that the extra-ionic quadrupolar interaction may be expressed with arbitrary coordi- 
nate axes and for any site symmetry. We write the extra-ionic quadrupole interaction 
similarly to the crystal-field interaction: 

K k = E KK ( 2 - 17 ) 

q=~2...2 

where the M 2 are the paremeters; and the C 2 are tensor operators like the but in 



the components of I. Following the procedure of Bunbury et al [fl^] and of Han |4C] 

2 

we can determine the M q : 



" Qjle 9 

where the are the crystal-field parameters (equation 2.7); Qj and Qj are defined 



below; 7 n and j e are the nuclear and electronic antishielding factors respectively, in 



2 

the notation of Edmonds |pl| . The operators C„ are 



i[3/ 2 2 -I(/ + l)] 

)^-[I z I ±+ I ±IZ ] 

^4 (2-19) 



where the (J || ai || J) are the operator-equivalent coefficients and I± = (I x ± I y ) as 





= (J\ 


«2 


r 2 


= T(J 


«2 


r 2 


= (J\ 


«2 



usual. Putting these operators into equation 2.17 gives 



K S q = %^ J W^W J){mI\wI-i{i + i)]-mI^[i z i + + i + i z ]+ 



MlAl z l„ + + mA\ + M^A . (2.20) 
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We define Qj and Qj [] as 



For calculation we use just the single parameter Q ex t, where 

Q ext = ^LllL. (2.22) 

Qj 7e 

It is commonly assumed that the antishielding factors are isotropic and there is no 
experimental evidence to the contrary. They are, however, host dependent. 

The extra-ionic quadrupole interaction may alternatively be written in terms of param- 
eters in the Baker, Bleaney and Hayes p| notation (extended to complex parameters): 



K k = rf-{B 2 [^-I(I+l)]+Bl\[I z I + +I + I Z }+ 

B?hl z I- + I-I z ] + Bl\l\ + B 2 2 \p_ \ . (2.23) 



In the case of orthorhombic symmetry, equation 2.23 reduces to the equation given 
by Han |4(| (our B\ = B 2 2 real, and equal to Han's B 2 ). For details of parameter 
conventions for crystal-fields see appendix A. 



2.5 Dipolar field 



The dipolar field is the field at a particular site resulting from all the other dipoles in 
the sample. We restrict this discussion to a paramagnetic sample but otherwise follow 
McCausland and Mackenzie [pCj ]. The dipolar field may be written as 



B 



dip 



E 



4-7T 



3(Pi- r i)ri _ Vi_ 
' i ' i 



(2.24) 



where the summation is over all other dipole moments /x,, at positions r« relative to the 

site of interest. Clearly, it is not possible to compute such a sum over a macroscopic 

1 These quantities are denned and used by Han [^o). However his definition (El equa- 
tions 2.34) includes extra factors of e/4 in both terms. Obviously this has no effect on the ratio 
Qi/Qji but the values obtained do not agree with the values for Qi and Qj given in Han's 
tables ( ]4C[| tables 1.1 and 2.2) which appear use the definitions given here. 
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sample. The normal procedure is to split the dipolar field into three components: 



Bdip — Bi n t + Bl + Bdm, (2.25) 

where Bi n t is the field resulting from dipole moments in a sphere centred on the site of 
interest, the 'Lorentz sphere'; Bl is the Lorentz field; and Bdm is the demagnetizing 



field due to the outer surface of the sample. Bi n t can be calculated from equation 2.24 



and this is referred to as a dipole sum, see section |2.5.l| . We make the approximation 
of uniform magnetization outside the sphere. In doing this we assume that effect on 
the dipolar field of the positions of individual dipole moments in the lattice is not 
significant outside the Lorentz sphere. A uniformly magnetized spherical shell has no 
field at its centre, Bl = —Bdm- Thus, the dipolar field at the centre of a macroscopic 
spherical sample may be calculated by considering a microscopic sphere: Bdi v = Bi nt 
because Bl = —B^ m . 

The power of the Lorentz sphere concept is in decoupling the effects of the local mag- 
netization from the demagnetizing effects of the sample surface. We may write Bl and 

Bdm 



B L =^M and B dm = —/j, HM, (2.26) 



where M is the magnetization at the surface of the Lorentz sphere; N is the second-rank 
demagnetization tensor and M is the average magnetization over the entire sample. 



In general, equation |2.26j is hard to apply. Only if the sample is an ellipsoid is the 
magnetization uniform througout the sample (and hence M = M) and N independent 
of position. Osborn |5(| tabulates demagnetization factors along the three axes of the 
general ellipsoid. Most magnetization measurements are interpreted by approximating 
the actual sample shape to an ellipsoid. The approximation is discussed by Crone- 
meyer J27|| . In the special case of a sphere, N reduces to a scalar N = |. Akishin 
and Gaganov ^ consider the demagnetizing effects in cylindrical and rectangular box 
samples, and list values for the diagonal components of the demagnetization tensor for 
various geometries. 
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2.5.1 Dipole sum 



The field B due to a collection of dipole moments iii at a relative positions r% is given 



by equation 2,24 , By resolving B into cartesian components we obtain the components 
of the dipolar field tensor: B xx , B yy , B zz , B xy = B yx , B xz = B zx and B yz = B zy . B xy 
denotes the component of the field in direction x resulting from moments in direction y. 

If the dipole moments are all the same (/i = fi^ for all i), then the diagonal components 
are 

B xx = y~] ( n ) and similarly for B yy and B zz ; (2.27) 

and the off-diagonal components are 

Bxy = 3xi]ji similarly for B yz and B xz . (2.28) 

i % 

A computer program has been written to compute the components of the dipolar field 
tensor for any arrangement of ions in a unit cell. The unit cell axes are each specified by 
a vector, and do not have to be orthogonal. Summations are performed over a sphere 
of given radius. The program is written in ANSI C, and on an IBM PC compatible 
computer (33 MHz '486DX) the summation runs at ~7000 ions/s. The calculation time 
increases linearly with the number of ions in the sum and hence as the radius cubed. 
For H0F3 the sum was found to converge to ~1% in 10 nm radius and to within 0.01% 
in 50 nm radius (percentage of the largest term). Calculation to 50 nm radius for H0F3 
corresponds to a sum over ~ 10 7 ions and is computed in less than half an hour. 



2.6 Rare-earth NMR 



This section outlines the relation of NMR to the energy levels of the ion. A careful and 



more detailed exposition of NMR in general is given by Slichter [72]; NMR in rare earths 
is considered by McCausland and Mackenzie [pT] ]. Consider the Zeeman interaction of 
a nucleus with moment fx in a magnetic field B. We may write the Hamiltonian: 

U = -fiQfi.B. (2.29) 
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This gives 21+1 equally spaced Zeeman levels: E = — fj, n B Mj where Mj = —I . . . I. 
An NMR experiment normally excites transitions between adjacent levels: the selection 
rule is AM/ = ±1. In rare-earth ions the dominant contribution to the 'field' at the 
nucleus is the 'effective hyperfine field' due to the 4/ electrons. 



If we assume that the hyperfine interaction does not modify J then equation |2.13 
becomes 

H' m ~ hA(J).I = fj, ogil i n B eff .I, (2.30) 
where B eff = hA(J) / '(/ioft fin) ■ 

For a set of (21 + 1) evenly spaced energy levels there would be just a single NMR 
frequency. In practice the nuclear energy levels are not quite evenly spaced because of 
quadrupolar interaction. This results in 21 different NMR frequencies. The quadrupo- 
lar interaction is usually much weaker than the dipolar interaction resulting in a spec- 
trum of 21 resonances where the separation of the resonances is a small fraction of 
their frequencies (see section 2.6.1| ). Unlike most NMR, in rare-earth NMR resolved 



quadrupole structure is the norm. Traditionally, an effective nuclear Hamiltonian is 
derived by treating TLhf a s a perturbation on the electronic Hamiltonian 7i e i (sec- 
tion |2.6.1|) . Although we diagonalise the combined Hamiltonian for the electrons and 



nucleus directly, the perturbation results are instructive. 
2.6.1 Effective nuclear Hamiltonian 

By considering the hyperfine interaction as a perturbation on the electronic Hamilto- 
nian an effective nuclear Hamiltonian is obtained: 

H n = a t I z + P t (l 2 z - ij(J + 1)) + wll (2.31) 

where the z axis is along the direction of (J) and off-diagonal quadrupolar terms 
have been neglected. They only have effect in second-order and it is expected to be 
small. The parameters at and Pt are combined intra- and extra-ionic parameters for the 
dipole and quadrupole interactions respectively. Thus, at = a' + a" and P t = P' + P" . 
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The parameter w is the pseudo-octupole term which arises from cross-coupling of the 
dipolar and quadrupolar terms when the latter is treated in second-order perturbation 
theory |50[|. 



From eqaution 2.31, the transition frequencies for | Mj) Mj — 1) are: 



VMuM.-x = a t + (2Mi - l)p t + (3M| - 3Mj + l)w, 



(2.32) 



where {—I + 1) < M/ < /. The effect of w is to create an asymmetry in the spectrum 
and also to shift the whole spectrum by u>/4. The pseudo-octupole term is usually 
very small (w < 1 MHz). The NMR spectrum of a fully polarized Ho 3+ ion ((J z ) = J) 
is shown in figure |2.2| . To give an idea of the frequencies involved the parameters are 
taken as at = a' Q , Pt = Pq and w = (see table [O]). 

If NMR spectra are obtained in the frequency domain then this model is very convenient 
for interpretation. Each spectrum can be fitted to some lineshape function centred on 



frequencies given by equation 2.32 to obtain at, Pt and w. NMR spectra taken by 
sweeping an applied field cannot be so easily interpreted because there is no simple 
relationship between the frequencies and the applied field. Each field corresponds to a 
different ionic Hamiltonian (7i e i + Ti-hf)- 



at 



2P, 2P r 



2P t 2P, 



2P t _ 



6120.8 6246.2 6371.6 6497.0 6622.4 6747.8 6873.2 v (MHz) 



Ho 3+ , I = 7/2, taking a t = a' Q = 6497 MHz, P t = P^ = 62.7 MHz and w = 0. 



Figure 2.2: NMR spectrum for a Ho 3+ ion. 
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165 Ho 3+ , ground manifold predominantly 5 Ig 
J = 8, I = | 

9j = 1.2417( 2 ), ft = 1.151 
Intra-ionic hyperfine paremeters: 

~ = 6l97 (8) MHz*- 1 ) ^ A = 812.1 (10) MHz 
6502 (6) MHzW =^ 812.8 (8) MHz 

= 62.7 (30) MHz^ =^ C = 0.523 (25) MHz 



Radial averages for the 4/ electrons ^ : 
(r~ s ) = 7.106xl0 31 m~ 3 
(r 2 ) = 2.085xl0~ 21 m 2 
(r 4 ) = 1.081 xl0~ 41 m 4 
(r 6 ) = 1.181xl0~ 61 m 6 



Reduced matrix elements for the ground manifold: 



Reduced matrix 


Russell- S aunders 


Intermediate 


element 


coupling^ 4 ) 


coupling^ 5 ) 


(J a 2 J) 


-2.222xl0~ 3 


-2.040 xlO" 3 


(J «4 J) 


-3.330xl0~ 5 


-3.082xl0~ 5 


(J a 6 J) 


-1.294xl0" 6 


-1.203xl0" 6 



Nuclear and electronic quadrupole moments: 

~Q n = 3.51(2) xlO" 28 m 2 ^ Qi = 1.67(1) x 10~ 29 m 2 

Q e = 5.104xl0~ 22 m 2 <= Qj = 4.254xl0" 24 m 2 ( 7 ) 



Ratio of nuclear to electronic antishielding factors: 

grille = 149(15) ( 6 ) in holmium hydroxide. 

246(10) ( 2 ) in holmium ethylsulphate. 
254 ( 8 ) in holmium aluminium garnet. 

1 Bleaney @. 

2 Carboni @. 

3 Freeman and Desclaux |5?J 

4 Abragam and Bleaney|l[] . 

5 Rajnak and Krupke |6^| . 

6 Bunbury et al @. 

7 Calculated from the (J || a 2 \\ J) of ( 5 ) and the (r 2 ) of < 3 ). 

8 McMorrow 



Table 2.1: Electronic and hyperfine parameters for holmium. 
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2.7 Low temperature heat capacity 



Heat capacity measurements in the helium temperature range can be dominated by the 
contribution from the electronic Schottky anomaly. In general, the electronic internal 
energy for a single ion U e i is given by 

i n E 

U el = -Y,Eie** (2.33) 
Z i=l 

where E% are the electronic energy levels and Z; T is the temperature; and Z is the 
partition function, 



E 

i=i 



-E, 



e~ . (2.34) 



Note that we consider only insulating compounds so there is no conduction electron 
contribution. Differentiating with respect to T at constant field gives the electronic 
Schottky heat capacity C e i, 

1=1 Ki=l ) 

In systems where the ground state and first excited state are well separated from higher 
excited states the heat capacity at low temperatures may be approximated by a two 
level system. For two levels E\ and E2, and writing j3 = ^ujr 1 , equation 2.35 simplifies 
to 

C el = kP 2 -^^. (2.36) 

Differentiating this expression to find the value of (3 corresponding to the maximum in 
C e i yields a transcendental equation in (3. Numerical solution gives maximum in C e i 
when = 2.399 . . . =¥ (E 2 - Ex) » 2.399kT max . 

The analysis of experimental data is not as simple as just finding the heat capacity 
maximum. Unfortunately, in the temperature range we are interested in (up to 20 K) 
the lattice heat capacity rises steeply. It is usually necessary to subtract the lattice 
heat capacity Ci a t oc T 3 before considering the electronic Schottky heat capacity. The 
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electronic Schottky heat capacity varies quite slowly around the maximum so small 
gradients resulting from other contributions may shift the apparent peak significantly. 



2.8 Bulk magnetization 



Given a suitable model for the rare-earth ion that predicts the electronic energy levels 
it is possible to calculate the bulk magnetization. 

M = hb9j 5" < i I J \i > -z- (2-37) 
i=i Z 

where the summation is over all populated states i. The Ei are the energies of the 
states and Z is the partition function (equation 2.34| ). Care must be taken if there 



are inequivalent sites in the unit cell as the individual moments may then differ in 
magnitude and direction. One of the more significant problems encountered when 
comparing predictions with experimental data is the effect of the shape of the sample 



on the demagnetizing field (see section 2.5). 
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3. CW NMR 



NMR was first observed by CW spectroscopy. However, when Hahn [39] discovered 
spin echoes, pulse techniques rapidly became the method of choice and CW NMR is 
now almost forgotten. In this chapter we argue a case for the use of CW NMR with 



some rare-earth systems. First, section |3.1| outlines some of the physical properties 
that affect NMR. 

Rare-earth NMR is technically difficult, principally because of the large hyperfine split- 
tings which give NMR frequencies of up to 7 GHz. This is often compounded by fast 
relaxation rates which mean that fast pulse sequences are required for pulsed NMR. 
This, in turn, requires short, high power pulses. The Manchester pulsed spectrome- 



ter [21, 49 1 provides up to 200 W of microwave power in pulses as short as 30 ns with 
rise and fall times of 5 ns. In spite of technical difficulties, pulsed NMR has been very 
successful in studying the hyperfine splittings of rare-earth compounds. However, for 
some systems fast relaxation can be a crippling problem. The sensitivity of CW NMR 
does not suffer from fast relaxation as badly as pulsed NMR. CW and pulsed NMR 



techniques are compared in section p.2j. 



Section |0] goes on to discuss the design of a CW NMR spectrometer. A peculiarity of 
NMR as used to study magnetism, as opposed to chemical structure, is that we often 
work in regimes where the relationship between the magnetization and the applied field 
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Figure 3.1: Graph illustrating a typical relationship between the line width in fre- 
quency Aui and the corresponding line width in field AB. 

is extremely non-linear. The intra-ionic dipolar interaction is the dominant contribu- 
tion to the hyperfme splitting, so the NMR frequencies are roughly proportional to the 
electronic moment. Hence the relationship between the NMR frequencies and the ap- 
plied field is often extremely non-linear. Partial saturation of the electronic magnetism 
is the norm and in extreme cases such as ferromagnetic materials the magnetization 
may respond only very weakly to applied fields. Line widths in rare-earth NMR vary 
from a few MHz to hundreds of MHz. Given the relationship between the ionic mo- 
ment (fj) and the applied field B for a given system, we can translate the line width 



in frequency, at any given field, into the corresponding line width in field. Figure |3.1 
illustrates this relationship. Typical values of the line width in field range from < 0.1 T 
to ~ 1 T, and the line width obviously increases as the differential susceptibility falls 
with saturation. This tends to make the NMR frequency insensitive to changes in the 
applied field at high fields. 



3.1 Spin dynamics and nuclear susceptibility 

This section outlines some of the factors that affect the design of an NMR experiment 



on rare-earth systems. Slichter [72] gives a much more detailed description of NMR 
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principles and techniques. We note that most of the common 'chemical NMR' tech- 
niques are unsuitable for rare-earth systems. This is principally because of the much 
faster relaxation rates and higher NMR frequencies encountered in rare-earth systems. 

In general we perform NMR on system of nuclear spins that are coupled the each other 
and to the lattice. The lattice acts as a thermal resorvoir and the lattice temperature 
determines the equilibrium populations of the nuclear levels. If the populations are 
temporarily disturbed, by NMR for example, they will then relax back to their equi- 
librium values exponentially with a characteristic time T\, the spin-lattice relaxation 
time. This relaxation time will impart a homogeneous line width to the resonance, with 
a Lorentzian line shape: 

91 H = n , Tprr \2v ( 3 - 1 ) 

7T(1 +T{(U) - UqY) 

where ujq is the resonance frequency. However, spin-lattice relaxation is not usually 
the dominant contribution to the homgeneous line width. In rare-earths T\ typically 
ranges from ~ 10 lis to ~ 10 ms at liquid-helium temperatures. Spin-lattice relaxation 
is intimately linked with saturation of NMR: see section |3.1.1| . 

Spin-spin relaxation is usually the dominant contribution to the lifetime of the 
individual nuclear spin eigenstates, and hence to the homogeneous line width: 
Atu/j ~ 7jr + tjt ~ 7jr, where T2 is the spin-spin relaxation time. In general the line 
shape is not Lorentzian (often it is more closely Gaussian) so the definition of T2 is not 
straightforward. However, for most of this chapter we assume exponential relaxation 
and hence a Lorentzian line shape, for which T2 is simply defined (cf. T\ and equa- 



tion |3.1| ). In section 3.3.1 we compare the derivative line shapes for Lorentzian and 
Gaussian lines. Relaxation processes are discussed more fully by Slichter [72] and by 
Abragam and Bleaney JIJ. 

Inhomogeneous broadening can be caused by spatial variations in the magnetic field 
and the electric field gradients; and by unresolved quadrupole structure. In rare-earth 
NMR the most common source of inhomogeneous broadening is physical inhomogeneity 
in the sample, either deliberate (alloying) or accidental (impurities, interstitials etc.). 
For broadening dominated by random variation of local fields we expect a Gaussian 



35 



line shape but in general the situation is more complex. It is useful to make the distinc- 
tion between microscopic and macroscopic inhomogeneous broadening. The borderline 
between the two is determined by the range of the spin-spin interaction. Typically 
macroscopic inhomogeneous broadening is caused by inhomogeneities in the applied 
field or by domain structure; and microscopic inhomogeneous broadening is caused by 
local impurity or interstitial effects. Inhomogeneous broadening does not limit the life- 
time of the nuclear spin eigenstates. If it is microscopic, it will actually increase T2 
because it reduces the coupling between spins. Significant inhomogeneous broadening 
is a requirement for spin-echo NMR. The spins must dephase and be re-phased by the 
second pulse within a time short or comparable to the spin-spin relaxation time, T2. 
We define the inhomogeneous line width Ao;^ similarly to the homogeneous line width. 

Both homogeneous and inhomogeneous broadening will destroy the phase coherence of 
the precessing nuclear magnetization and it is useful to consider a combined dephasing 
time, 

T 2 ~ -7 • (3-2) 

From this quantity we may restate the requirement for spin echo NMR as T^IT^- In 
chemical NMR there is often insufficient intrinsic inhomogenous broadening for spin 
echo experiments so field gradient coils are used to produce macroscopic inhomogeneous 
broadening (see, for example, Slichter |fT2[). However, in rare-earth systems inhomo- 
geneous broadening usually dominates anyway. When that is not the case, T2 is often 
so short that spin-echo NMR becomes impracticable, even if enough inhomogeneous 
broadening were produced by the application of a large field gradient. Free-precession [] 
or CW NMR are then the only practical options. 

A very important consideration in NMR is the enhancement of the nuclear magnetism 
by its coupling to the electronic magnetization of the ion. We note that 'enhancement' 
has two distinct meanings in this context: there is enhancement of the effective hy- 

perfine field which results in microwave NMR frequencies for the rare-earths, and also 

1 In this work we use 'free precession' to refer to the initial free precession or free-induction 
decay (FID) signal after pulsed excitation. Although a spin-echo signal is also caused by the 
free precession of the nuclear spins we exclude spin echoes from the term free precession. 



36 



enhancement of the transverse RF nuclear susceptibility which affects the strength of 
the NMR signal. Here, we consider enhancement of the transverse RF nuclear sus- 
ceptibility. Further details are given by McCausland and Mackenzie p0[ . First, the 
transverse RF field (ie. perpendicular to the electronic magnetization) 'seen' by the 
nucleus is enhanced by the response of the electrons. For 'small' transverse RF fields, 
the transverse electronic magnetization is proportional to the applied transverse field 
to a good approximation, and we may write the RF field seen by the nucleus as 

B[ = (l + r,)B 1 , (3.3) 

where n, the enhancement factor can be anywhere between 1 and 10 4 ; and B\ is the 
amplitude of the applied transverse RF field. Secondly, the response of the system is 
not just the precessing nuclear magnetization, but the sum of the electronic and nuclear 
magnetizations. The amplitude of the combined transverse magnetization, m'i is given 
by 

m'i = (1 + n)mi, (3-4) 

where mi is the transverse component of the precessing nuclear magnetization. Thus 
the combined effect is to enhance the NMR signal by the factor (1 + n) 2 . We may 
alternatively express this as an enhancement of the transverse nuclear susceptibility by 
the same factor: 

Xs = (l + r/) 2 X n (3.5) 

where Xn is the transverse nuclear susceptibilty, B[/mi; and Xs is the total transverse 
RF susceptibility of the ion resulting from the nuclear susceptibility. The transverse 
RF nuclear susceptibility is discussed further in section 3.1.2| . 



3.1.1 Saturation 

The onset of NMR saturation is determined by the spin-lattice relaxation time T\ and 
by the intensity of the RF field. Saturation is caused by deviation of the populations 
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of the nuclear levels from their equilibrium values. If T\ were infinite then any RF 
excitation would equalise the populations of the nuclear levels and absorption would 
cease. The power absorption for a two-level system on resonance is given by 

f = no ^TTW^' (3 - 6) 

where no is the equilibrium population difference; W is the RF induced transition rate; 
and uj is the frequency. Immediately, we see that for negligible saturation we require 
W <C 2^"- For 'weak' RF excitation W is given by 

W = 7T 7 ^M(^)I(1|4|2)| 2 J (3.7) 

where B x is the RF excitation; g(oS) is the line shape function; and 7„ is the nuclear 
gyromagnetic ratio. This expression can be generalised to a multilevel system by re- 
placing the matrix elements (1 | I x | 2) with {m n \ I x \ m n i) where the m n are the 
eigenstates of the unperturbed nuclear Hamiltonian. Equation does not hold for 
RF excitation that is 'strong' compared with the interactions responsible for the ho- 
mogeneous line width. However, in rare earth solids T\ is usually much longer than T2 
and the NMR will be saturated by 'weak' excitation. 

It is interesting to recall that the first attempt to see NMR, by C G Gorter, failed 
because the sample had an extremely long T\ so the resonance was saturated (see, for 
example, Slichter |7^|). 



3.1.2 Transverse RF nuclear susceptibility 

If we assume that both the spin-lattice and spin-spin relaxation processes are expo- 
nential we may use the Bloch equations (see, for example, Slichter J72| ). These give a 
Lorentzian NMR absorption line shape: 




where oj s is the sample resonance frequency and Aw s is the line width. Here we have 
neglected inhomogeneous broadening which, in general, will result in just a portion of 
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the spin system being excited for any given set of conditions. It is convenient to define 
the quantity x, the deviation from resonance in units of half the line width, as 

<-> 

Then the real and imaginary parts of the transverse RF nuclear susceptibility, 
Xn = x' n ~ i-Xn, around sample resonance are given by 

x ' n = X0 it s TTx^ and (3 ' 10) 
x " = Xo a^TT^' (3 - n) 

where xo is the static nuclear susceptibility. The real part of Xn represents dispersion 
and the imaginary part represents absorption. The sign of x'n i s chosen such that 
positive x'n indicates absorption. 

In most real situations, the sample will not entirely fill the experimental cell. To model 
the effect of the nuclear susceptibility on the whole cell it is convenient to define the 
'filling factor', 



Bldv 

sample ^ 

B^dv 

cell 

where B x is the amplitude of the the RF field. It is often the case that the sample 
is small compared with the cell volume. Then it is more convenient to consider an 
'effective cell volume', 

1 r 9 



B 



v c = ~52 B *dv, (3.13) 



l J cell 



where B\ is the value of B x at the sample (assumed constant over the sample). Then the 
filling factor is simply, q = v s /v c where v s is the sample volume. Taking enhancement 
and the filling factor into account, the effect of the sample RF nuclear susceptibility is 
equivalent to the cell being filled with a material of RF susceptibility given by 

x = q(l+7 1 ) 2 (x'n-iXn), (3-14) 
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where we assume that the NMR is not saturated. Depending on the experimental 
configuration, NMR can be detected by absorption, dispersion or a mixture of the two. 



3.2 CW and pulsed NMR compared 

Both CW and pulsed NMR have been used in this work, and in the case of H0F3 
both techniques give adequate NMR spectra. However, in this chapter we consider 
rare-earth NMR in general. Spin-echo NMR is not possible when broadening is purely 
homogeneous; free precession following a single pulse is then the only viable pulse 
technique. However, spin-echo and free-precession NMR suffer a similar degradation 
in sensitivity as the relaxation times (Ti and T2) decrease. 

The theoretical signal to noise ratios obtainable from CW and pulsed NMR have been 
considered by McCausland and Mackenzie fl50|| . For pulsed NMR the maximum signal 
to noise ratio is 



S\ n /i A 2 (uv a f(m)(l + rj)m 



JV/max kT R T x v c \ Aio ih 



(3.15) 



where v c is the effecitve cell volume; v s is the sample volume; f(m)(l + r])mo is the 
peak value of the RF magnetization; t, is the integration time; and Tr is the effective 



noise temperature of the receiver. Equation 3.15 assumes dominant inhomogeneous 
broadening. The factor A allows for relaxation during the experiment: for spin echo 
A = e* -2 * 12 ^ 2 ) where t\2 is the pulse separation, and for free precession A = e(~* ld / T 2' 
where t%d is the time between the pulse and detection. Ideally A w 1 but in practice 
this may be hard to achieve if T2 is short. This assumes that the resonator Q has 
been optimised so that the ring time is about the same as the pulse length, and the 
excitation power tips all the spins by ir/2. 



For CW NMR the maximum signal to noise ratio is 

1 2 

(3.16) 



n n T c f uv s f(in)(l + r])m ^ 2 



N/max kT R TiT 2 v c { Au ih 

where r c is the ring time of the resonant cell; and the other parameters are as in 
equation 3.15[ Comparing equations 3.15 and 3.1€ shows that CW will be more sensitive 
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if r c > T2, ie. the resonator ring time is longer than the spin-spin relaxation time. This 
implies that we want as high a Q as possible for CW NMR, but high Q can also bring 
technical difficulties: see section |3T 



In fast relaxing systems, T2 is the eventual downfall of pulse techniques. The factor A 2 
can significantly decrease the signal to noise ratio if the spectrometer cannot produce 
a fast enough pulse sequence. Even if the spectrometer can switch fast enough, the 
pulses must be short compared with T^. The shortness of the pulses will impart a 
frequency uncertainty giving an instrumental broadening of Au oc 1/t p where t p is the 
pulse length. To achieve efficient excitation with short pulses, high power is required. 
Assuming that the resonator Q is limited by acceptable ring time, the power required 
goes as T 2 ~ 3 . Typically, the power required becomes prohibitive when T2 is much less 
than 50 ns. Fast relaxation (short T2) also reduces the sensitivity of CW NMR but this 
is not compounded with the increasing technical difficulties associated with fast pulsed 
experiments. 

In general, pulse techniques suffer less from background fluctuations, but there can be 
problems with spurious echo-like signals. Also the spin-echo technique allows straight- 
forward measurements of T2, T 2 * and, with minor modifications, T\. From the expe- 
rience of this work, the main problems with CW NMR seem to be background fluc- 
tuations and the complicated line shapes obtained with some modulation strategies. 
Field-modulated CW NMR gives easily intepreted derivative spectra. Frequency mod- 
ulation is fraught with technical difficulties but may be the only option in systems 
where the hyperfine splitting is insensitive to the applied field. Modulation strategies 



for CW NMR are considered in section 3.3 and in chapter 4. 



Unusually, the problem with pulsed NMR in H0F3 is not that T2 is particularly short, 
but instead that T\ is long (seconds). This limits the repetition rate of pulsed experi- 
ments, thus decreasing the signal to noise ratio. For H0F3 we have obtained the best 
resolved spectra from field-modulated CW NMR: see chapter 5. 

To conclude, pulse techniques appear to be the best approach in systems where the 
relaxation is not too fast, particularly if the spin-echo technique can be used. However, 
when T2 is very short (say < 100 ns) the balance tips in favour of CW NMR. 
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3.3 Considerations in CW spectrometer design 



In this section we discuss some of the considerations applicable to the design of a CW 
microwave NMR spectrometer. As the section progresses we focus on the configuration 
chosen for this work. 

Irrespective of the detection strategy, some form of 'sweep' is required to find NMR. In 
addition, some form of modulation, though not strictly essential, is required to overcome 
1/f noise. In principle, one can sweep (and modulate) the field, the frequency or, in 
exceptional cases, the temperature (see Cowan and Cha [j26|). In practice, sweeping the 
frequency has severe problems arising from the non-uniform frequency response of the 
spectrometer. Field sweep, combined with field modulation, is therefore the method 
normally adopted for CW NMR. 

Even in situations where the hyperfine splitting has a significant field dependence, the 
NMR line width, in field units, is rarely less than 0.1 T and may well be of the order 
of 1 T or more when the magnetization is saturated. Since efficient detection requires 
a modulation amplitude which is a significant fraction of the line width, this places 
severe demands on the modulation source. Under such circumstances, a combination 
of field sweep and frequency modulation may be the best solution. 

Provided that technical problems can be overcome, sensitivity will ultimately be limited 
by noise. The modulation frequency should be be high enough to ensure that 1 // noise is 
not the limiting factor. However, we have observed NMR in H0F3 without modulation 
(except a small frequency modulation to lock the microwave transmitter frequency to 
the resonator): NMR was clearly seen both by dispersion (resonator frequency shift) 
and by microwave absorption (change in the resonator Q). 

Fundamental to the design of a CW NMR experiment is the decision whether to use a 
resonant or non-resonant cell. This question is discussed in section |3.3.2 . The use of 



a non-resonant cell is technically much easier but suffers a severe penalty, in terms of 
detection sensitivity, associated with an effective Q ~ 1. The most common detector is 
a microwave diode, although several other options are available and are briefly discussed 



in section 3.3.2 



42 



3.3.1 NMR line shape and derivative spectroscopy 

In this section we compare Lorentzian and Gaussian line shapes and their derivatives. 
Lorentzian lines are typical of NMR in liquids and are implicit in the Bloch equations; 
Gaussian lines are typical for solids with microscopic inhomogeneous broadening. In 
general, the NMR line shape will be neither Lorentzian nor Gaussian, but the simple 
analytic forms are convenient for calculations. The qualitative results obtained in this 
section do not depend critically on the the line shape, and we assume the Lorentzian 
line shape for convenience. 

We express line shape functions as absorption curves, normalised such that g(0) = 1, 
in terms of x which has units of the fractional half line width as defined earlier (so 
g(±l) = 1/2). We may write the Lorentzian line shape as 

*(*) = TT^' (3 - 17) 

and the Gaussian as 



g g (x)=e- l ^ x . (3.18) 



Figure 3.2 shows these functions and the first, second and third derivatives with respect 
to x. Higher derivatives of the Lorentzian line tend to get sharper whilst those of 
the Gaussian remain of similar width to the original function. However the general 
character of the derivatives is the same for both line shapes. 

Modulation of either the frequency or field may, with suitable scaling, be represented as 
a modulation of the parameter x. In the limiting case of small modulation (Ax m <C 1), 
the n th harmonic response will give the n th derivative. In this work we refer to the 
fundamental as the 1 st harmonic, twice the fundamental frequency as the 2 nd harmonic 
and so on. Note however, that the amplitude of the n th harmonic signal is proportional 
to (Ax m ) n , where Ax m is the modulation width. 
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Figure 3.2: NMR line shape and derivatives for Lorentzian and Gaussian lines. The 
solid lines are from the Gaussian g g (x) and the broken lines are from the Lorentzian 
gi(x). The upper graph shows the functions g g (x) and gi(x); below are the first, second 
and third derivatives with respect to x, denoted g', g" and g'" respectively. 
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3.3.2 Sample cells 



The most common CW configuration for radio-frequency operation is the 'marginal 
oscillator' which uses a resonant circuit containing the sample to form part of the 
oscillator. This is typically achieved by putting the sample inside the inductor of an 
LC tuned circuit. Resonant absorption causes a small reduction in the quality factor 
of the circuit, resulting in a change in the output amplitude. Modulation of the RF or 
of the applied field is used to overcome noise and drift. 

At microwave frequencies it is not possible to use lumped circuits; instead distributed 



elements are used (waveguides, coaxial resonators etc.). Figure 3.2 shows two possible 



experimental configurations. Coaxial resonators operating as reflection cells were used 



to obtain all the data reported in this work. As shown in figure 3.3, a circulator is 
used to direct the incident microwave power to the reflection cell and the reflected 
power to the detector. Cavity resonators operating as transmission cells with separate 



drive and detection ports have been used for CW NMR by Bleaney et al [11]. Cavity 
resonators and a circulator are also used in the Manchester pulsed spectrometer ]2l| ], 
albeit with resonator Q limited by the maximum acceptable ring-time. A significant 
advantage of reflection cells is they have only one coupling to a transmission line which 
makes adjustment of the coupling easier than for transmission cells where there are two 
couplings. 

A non-resonant transmission cell may be useful when large amounts of material are 
available. The best approach is to pack a section of transmission line or waveguide 
with sample, partially replacing the dielectric. NMR is detected by the decrease in 



transmitted power at resonance. Cowan and Cha 26] have studied terbium metal 
using a transmission cell with a temperature sweep at fixed frequency. Ferromagnetic 
HoFe2 has been studied by Ross |37j] as a powder in a transmission cell using frequency 
sweep. 



Cowan and Cha [26| used finely ground Tb powder to maximise the effective sample 
volume which is otherwise limited by the skin depth. They managed to measure the hy- 
perfine splittings at temperatures up to 160 K, exceptionally high for rare-earth NMR. 
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Figure 3.3: Two possible CW NMR configurations: a) reflection cell and b) transmis- 
sion cell. In both configurations the cell may be either resonant or non-resonant. 

The line widths observed ranged from 10-50 MHz which gave well resolved quadrupole 
structure; the quadrupole splitting of ferromagnetic Tb at 4.2 K is > 600 MHz with 
the central line at « 3.1 GHz. Initially, they used frequency sweeps but found prob- 
lems with power and backgound fluctuations. They found that the background was 
reduced by over an order of magnitude by sweeping the temperature whilst keeping the 
frequency constant. Sweeping the temperature effectively sweeps the field at the nuclei 
because the thermally averaged electronic moment changes. In Tb the lifetime of the 
electronic states at temperatures above 80 K is short compared to the period of the 
nuclear precession. This means that the nucleus experiences an effective hyperfine field 
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proportional to the thermal average of the electronic moment, which decreases with 
increasing temperature. In the present work, which is concerned with insulating rare- 
earth compounds at liquid helium temperatures, the lifetime of the electronic states is 
long compared to the period of nuclear precession and the effective hyperfine field is 
determined by the electronic ground state (see Bunbury et al 

A major disadvantage of non-resonant over resonant cells is that of reduced sensitivity. 
To achieve equal sensitivity, a completely filled non-resonant cell with length of order 
qQX is required to compete with a resonant cell of quality factor Q and filling factor q 
(A is the wavelength of the microwaves in the cell). The required sample volume is often 
unacceptably large. Powdered samples often used and generally the particles are ran- 
domly oriented. Thus, any externally applied field is randomly oriented with respect to 
the crystal axes of the individual particles. In general, this will result in different NMR 
spectra from the individual particles and hence in gross macroscopic inhomogeneous 
broadening. Thus, powder samples are usually used only for spontaneously magnetized 
materials. Non-resonant cells are clearly unsuitable for small single crystal samples. In 
that case a resonator with small effective cell volume and high Q provides the greatest 
sensitivity. 

Whilst offering high sensitivity, resonant cells have associated technical difficulties, 
especially when the Q is high. First, the microwave frequency must be kept close to 
the resonator frequency. This may require some form of servo system either for the 
resonator tuning or the microwave frequency. Also, if frequency modulation is used the 
response of the resonator may dominate the received signal. However, these problems 



are surmountable: see section 3.3.8 



3.3.3 Detection systems 

Microwave diode detectors produce an output voltage approximately proportional to 
the input power for low input power. As the input power increases they gradually 
change from 'square law' to 'linear law' with the output power proportional to the 
input power, ie. with an output voltage proportional to the square root of the input 
power (assuming a constant load impedance). Diode detectors are reliable and easy 
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to use. A more sensitive method of detection is to use a low-temperature bolometer. 
Si and InSb bolometers are commercially available with high sensitivity at helium 
temperatures. However, their response is slow and would limit the modulation rate to 
about 500 Hz. 

It is also possible to detect NMR by the temperature rise of the sample caused by 
resonant absorption. The temperature can be measured by mounting the sample on 
a thermistor. At helium temperatures the temperature coefficient can be very high so 
it is possible to detect very small temperature changes. Sensitivity to NMR obviously 
depends on the strength of resonance and on the spin-lattice relaxation time T±. 

Whatever the detector, it is desirable to use a lock-in amplifier referenced to the mod- 
ulation frequency or some harmonic. Clearly, if the response time of the detection 
system is long then the modulation frequency must be limited. In this work we have 
used a diode detector, and detector noise was not the sensitivity limiting factor. 



3.3.4 Quarter-wave coaxial resonator 

All the NMR data in this thesis were taken using tunable truncated quarter-wave 
coaxial resonators. First, consider the simplified section through a (non-truncated) 



quarter-wave coaxial resonator shown in figure 3.4 




Central conductor 
Optimal sample position 



z=0 



Figure 3.4: Simplified section through a quarter-wave coaxial resonator. 



48 



Ignoring end effects we may write down the current in the central conductor as a 
function of position z for the lowest frequency mode (I = A/4): 

J(s)=I cos(^) (3.19) 

and hence the magnetic field distribution: 

5(r ' z) = ^ /oCOS (~r) (3 - 20) 

We immediately see that the best position for a small sample is on the surface of 
the central conductor at the closed end of the resonator. In the previous work of 



McMorrow [51] and of Carboni et al [21], metallic single crystals have been used as 
the central conductor. With a metallic sample the active sample volume is severely 
limited by the skin effect. At GHz frequencies the skin depth is typically 1 fim at low 
temperatures. Using the sample as the central conductor places the active part at the 
ideal position in the cavity. 

To make a resonator for frequency uq the central conductor length is given by 



iK g (3.21) 



where c is the microwave propagation speed in the resonator. To a good approximation, 
if the resonator is full of liquid helium then \i rj //q and e ~ eo so c as cq = 1/ JJIq^q. 



For resonance at 5 GHz, equation 3.21 gives a central conductor length of ~ 15 mm. 
Note, however, that sample resonance can cause a significant change in the effective \x 
averaged over the resonator volume which changes the resonant frequency. 

The resonator can be made tunable by adding a piston at the 'open' end which provides 
a capacitance at the end of the central conductor. We now call it a truncated quarter- 
wave coaxial resonator. If the separation d between the end of the central conductor 
and the piston is small compared to the radius of the end of the central conductor then, 
ignoring end effects, 

d = eAZ Lj tan , (3-22) 
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where A is the end area of the central conductor; and 



Z =l[Eln(-\. (3.2:!) 



2ir V e \° 

The expression for d is due to Carboni who plots theoretical curves against mea- 
sured values showing reasonable agreement. 

It is convenient to use the concept of the 'line width' of the resonator which we simply 
define as 

Au r = (3.24) 

where ujq is the resonant frequency and Q is the quality factor. The Q and hence the 
line width Aoj r depend on the coupling to the resonator and on any resonant absorption 
inside the resonator. We use the term intrinsic Q to refer to the Q of the resonator 
in isolation; and loaded Q to refer to the Q of the resonator coupled to a transmission 
line. Both of these quantities are defined in the absence of NMR and we subsequently 
consider any resonant absorption as a perturbation. 

3.3.5 Parallel-LCR resonator model 

To model NMR in a resonator, it is convenient to consider an electrical analogue of the 



microwave resonator. Figure 3.5 shows a parallel-LCR circuit where the sample NMR 



can be modelled by a change in the value of L: see section 3.3.7. 



The resonator must be coupled to the microwave source and to the detector. This can 
be achieved with one or two transmission lines terminated with 'loops' protruding into 
the resonator so that the current is coupled to the RF magnetic field in the resonator. 
In this work we have used a reflection cell so we consider the case of a single coupling 
loop. The coupling of the transmission line to the resonator provides an impedance 
conversion. This is accounted for by transforming the impedance of the resonator 
analogue such that we can consider it as if it were directly connected to the transmission 
line. 



50 









Ffcoj 


R 


C 




■4 [ 

















L=L Q (1+ %) 



Figure 3.5: Parellel-LCR analogue of a resonator. 

To define the resonator parameters, consider the system in the absence of NMR. Any 
non-resonant losses in the sample are included in the quality factor Q and are effectively 
lumped into R. The inductance, L = Lq 'off-NMR' and the resonant frequency is simply 

■* = vhu- (3 - 25> 

3.3.6 Coupling to a resonator 

We first consider just the resonator when the sample is well off NMR resonance. If the 
resonance frequency of the resonator is ujq; the intrinsic Q is Qo] and the inductance 
in the equivalent parallel-LCK circuit is Lq then, on resonance, the admittance is 



Y = YM = -j-r- . (3.26) 

If the characterisctic admittance of the transmission line is Y t , we define the coupling 
coefficient a as, 

- - | (3.27) 
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The regime a < 1 is referred to as under-coupling; a = 1 as critical coupling; and a > 1 
as overcoupling. If two transmission lines are used the transfer characteristic will be 
more complex than the case of a single transmission line, involving both couplings. Here 



we restrict discussion to single port resonators as used for this work. From equation |3.27 
the voltage reflection coefficient p(uj) follows as 



p(u) 



aY + Y{co)' 



(3.28) 




■4-3-2-10 1 2 3 4 
Detuning parameter, Q (w - w)/'u? B 



Figure 3.6: Power reflection coefficients | p 2 
function of the detuning. 



for different degrees of coupling 
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Figure 3.6 shows the power reflection coefficient | p 2 (co) | as a function of the detuning. 
The loaded resonator line width is increased by over-coupling but is decreased by under- 
coupling, relative to critical coupling. The Q of a critically coupled resonator is half 
that of the same resonator unloaded (ie. in the limit infinite under-coupling a — > 0, 
when Q — > Qq). 

In practice the coupling depends critically on both the position and orientation of the 
coupling loop. Adjustment is achieved by sliding and twisting the transmission line 
with the coupling loop at the end. However, it is hard to know what the coupling 
is except in the region of critical coupling where the reflected power falls to zero on 
resonance. It is often not even easy to see whether the system is slightly under- or 
over-coupled. 



3.3.7 NMR in a resonator 



This section discusses the effect of NMR on the response of a resonant reflection cell. 
The transverse RF nuclear susceptibility of the sample Xs, can be represented by a 
change in L: 



L = L (l + q Xs ), (3.29) 

where q is the filling factor and Xs = x's ~ iXs- Expanding the inductance to an 
explicitly complex form gives 

L = L (l + qx' s -iqXs)- (3-30) 



The imaginary part of the susceptibility x's 1S chosen with the negative sign so that 
positive x's represents absorption. This can be seen by considering the impedance of 
L: 



Z L = uL oqX 's + iwLo(l + qx's), (3-31) 

where positive x's gi yes a positive resistance. It is convenient to transform this series 
impedance into a parallel admittance. Assuming | qx s I "C 1 and working to first order 
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in x' s an d Xs gives 

Y L « ^ - fc^l . (3.32) 

This is equivalent to a resistance (ojLq) / (qx") in parallel with an inductance Lq{1 + qx' s )- 
Having performed this transformation the admittance of the resonator follows simply, 
and dividing by Yq gives 

rM„f 1 + q***a\ +ieh (<J'-<4<i-rf)). (3.33) 



To find the transfer charactersic of the system we must consider this admittance con- 
nected to the transmission line. The voltage reflection coefficient is given by equa- 
tion 3.28 and | p 2 (u>) | is the power reflection coefficient. In an experiment using a 
resonator and square law detector, the output voltage is proportion to | p 2 {u>) | pro- 
vided the input power remains constant. 



Figures [3.7| and 3.£ show the effect of a NMR resonance that is narrow relative to the 
resonator line width. The intrinsic Q of the resonator Q r is taken as 100 which gives 
a loaded Q ~ 50 around critical coupling. The sample line width Auj s = lo/Q s where 
it is convenient to consider the sample 'quality factor' Q s for comparison with that of 
the resonator. We have taken Q s = 1000 = 10Q r have chosen the resonance frequency 
of the resonator as uq/2ti = 5000 MHz so that the frequencies are of the correct order 
for rare-earth NMR. The strength of the NMR is deliberately exaggerated. 

The optimal value of the coupling coefficient a depends on the type of detector. If the 
detector is 'linear' then the NMR signal is maximised for critical coupling (a = 1); 
if it is 'square-law' then a = 2 ± \/3 is optimal. In practice we have used coupling 
fairly close to critical, approximately 0.5 < a < 2. Here we take a = 0.5 and a = 2 



as examples. Figure 3.7 is for an under-coupled resonator, and figure 3.8 is for an 
over-coupled resonator. In both cases the solid line is the resonator resonance curve 
without NMR (ujq/2tt = 5000); three other curves are shown for NMR resonances 
at uj s /2it = 4980, 5000 and 5010 MHz. Although the curves for over- and under- 
coupling are significantly different, they both show that the NMR resonance creates a 
sharp and, unless uj s = ujq, asymmetric change in the reflection coefficient. Consider a 
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Figure 3.7: Power reflection coefficients as a function of frequency for a cavity with 
intrinsic Q r = 100 and a sample Q s = 1000 with the cavity under-coupled (a = 0.5). 



frequency modulation experiment where the carrier frequency is somehow kept tuned 
to the resonator resonance. It is clear that varying amounts of harmonics will be 
generated from a sinusoidal frequency modulation as the sample resonance frequency 



moves relative to the resonator frequency: see section 3.3.8. Although this case usefully 



illustrates the harmonic generation, in practice we have used resonators with Q r higher 
than Q s for the systems studied. 



Figures 3.£ and 3.1C| show the more realistic situation of Q r = 1000 and Q s = 100. 
Except for the swapping the line widths of the resonator and NMR, the other param- 



eters are as for figures 3.7 and 3.5. With the NMR resonance much broader than the 
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Figure 3.8: Power reflection coefficients as a function of frequency for a resonator with 
intrinsic Q r = 100 and a sample Q s = 1000 with the cavity over-coupled (a = 2). 

resonator resonance we cannot expect sharp variations in the reflection coefficient to 
be caused by NMR. However, there are definite changes in the magnitude of the reflec- 
tion coefficient and in the resonance frequency (discussed below). As before, the power 
reflection coefficient is increased for under-coupling but decreased for over-coupling. 
There is also a small asymmetry that will generate odd harmonics from a sinusoidal 



frequency modulation about the centre of the combined resonance: see section 3.3.8 



From figures 3.9 and 3.10 there is clearly a shift in the cavity resonance frequency when 
the sample resonance is close to it. This results primarily from the dispersive part of 
the sample susceptibility. The the cavity resonance occurs when the admittances of 
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Figure 3.9: Power reflection coefficients as a function of frequency for a resonator with 
intrinsic Q r = 1000 and a sample Q s = 100 with the cavity under-coupled (a = 0.5). 

the inductance and the capacitance cancel. 

To first order in x' s f ne cavity resonance frequency is: 

This expression is useful only when the sample line width is greater than the resonator 
line width so that x's is approximately constant over the centre portion of the resonator 
resonance. 
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Figure 3.10: Power reflection coefficients as a function of frequency for a resonator 
with intrinsic Q r = 1000 and a sample Q s = 100 with the cavity over-coupled (a = 2). 



Recalling equation 3.10, around the sample resonance we expect the resonator frequency 
to be increased when the uj s > ujq; equal the natural resonant frequency when uj r = 
uj s = loq; and be reduced when u s < loq- These inequalities are not affected by the 



strength of the coupling. We have used this method to detect NMR, see section 4.7.3. 



3.3.8 Frequency modulation with a resonator 

The use of frequency modulation with a resonator is unusual but has given good results 
(see chapter 5). If the NMR line width is small compared to the resonator line width 
it is easy to see how harmonics of the sinusoidal frequency modulation are generated 
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by the NMR resonance. If the modulation amplitude is comparable to the NMR line 
width it will differentiate the 'sharp' NMR line from the broader resonator resonance. 
Although this situation could be engineered by designing a low Q resonator, the low Q 
will reduce the sensitivity of the spectrometer. In our experiments we have opted for 
a high Q resonator so that the resonator line width is much less than the NMR line 
width. In that case it is not practical to use a modulation amplitude comparable to 
the NMR line width. 

If we keep the carrier frequency at the 'centre' of the combined resonator and sample 
resonance, in the absence of NMR the symmetric resonator line will generate only even 
harmonics of the modulation and the fundemental will be supressed. The harmonic 
response of the resonator is essentially like the derivatives of a Lorentzian NMR absorp- 
tion (see figure j3.2|) . As the NMR resonance is swept though the resonator resonance, 
or conversely, the overall line shape will be asymmetric and even and odd harmonics 
will be generated. However, because the response of the system is dominated by the 
roughly symmetric resonator resonance which generates mainly even harmonics, we 
choose to detect NMR using an odd harmonic of the modulation frequency. 

In practice, it is necessary to use feedback to keep the microwave frequency accurately 
to the combined resonator and sample resonance. In the absence of NMR, all odd 
harmonics of the resonator response are zero at the resonant frequency if we assume a 
symmetric resonator resonance. Thus, using the first harmonic as an error signal, the 
carrier frequency of the microwave oscillator can be locked to the combined resonance. 
NMR will cause a shift in the frequency at which the first harmonic signal is zero 



(see figures 3.9 and 3.10) and thus in the locked carrier frequency. Obviously, this 
means that the first harmonic cannot be used to detect NMR. However, the NMR also 
creates an asymmetry in the combined resonator and sample response which generates 
higher harmonics that can be used to detect NMR. In the regime Q r ^> Q s it is not 
practicable to use a modulation deviation that is comparable to the NMR line width. 
The resonator line width puts a limit on the maximum practical modulation deviation. 
Thus Aw m <C Aiv s , which imposes a significant sensitivity penalty on the use of higher 
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harmonics: approximately, 



sensitivity oc 



Au m 
Alu s 



) 



n 



(3.35) 



where Aw m is the frequency modulation amplitude, Au s is the NMR line width, and 
n is the detection harmonic. Thus third harmonic detection is the best option when 

Qr Qs' 

Ideally the third harmonic signal will be zero when off the NMR resonance. However, 
this is not the case in practice because of imperfect locking to the resonator and asym- 
metries in the overall transmission function of the microwave system. Fluctuations in 
the 'background' third harmonic signal are a significant contribution to the total noise. 
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4. The CW microwave NMR 
spectrometer 



Some of the requirements of a CW spectrometer for NMR spectroscopy of rare-earth 
ions in solids have been discussed in chapter 3. In particular, the relative merits of 
'sweep' and modulation in the field and frequency domains have been considered. The 
spectrometer has been built to allow any combination of sweep and modulation. Also, 
there are facilities to lock the carrier frequency to the combined resonator and sample 
resonance, and to allow harmonic detection. 

This chapter describes the spectrometer as built. Although we concentrate on field 
sweep experiments, the spectometer has been designed with flexibility in mind. The 
spectrometer architecture and arrangements for field and frequency modulation exper- 



iments is described in section 4.1. The microwave system and circuits used to set up 
the spectrometer are described in section [L^. Computer control has been implemented 
from the start and has influenced the organisation of the system. In particular, parts of 



the system controlled by the computer are grouped in the 'interface unit' (section |4.3| ). 

The current operating frequency range of the spectrometer is 4-8 GHz. However, it 
would be straightforward to extend the frequency range to 2-8 GHz: simple modifi- 
cations to the microwave system and to the software would be required. These are 
described in the appropriate sections. 
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The spectrometer has been used successfully only with field sweeps (section [4.7.1|) . 
Frequency sweeps have been attempted on several occasions but the background has 



proved very much larger than any expected signal. Section 4.7.2 discusses the frequency 
sweep method a little further but otherwise we consider only the field sweep method. 
Also, and as a counter-example to the accepted wisdom that modulation is required, 



section 4.7.3 reports detection of NMR by the shift in the resonator frequency and by 
the change in reflected power, both without modulation. 



4.1 Spectrometer architecture 



First, consider the experimental apparatus. Figure 4.1 shows the main parts of the 
apparatus and the principal interconnections. The computer is an IBM PC compatible 
with two interface cards: an IEEE488 bus interface (Brain Boxes Ltd., 'Elite' card) 
and a parallel input/output (PIO) card. The PIO card is used to drive the instrument 
interface via a buffering and signal conditioning card. The interface, the instrument 



unit and the cards it contains are described in section 4.3. The microwave system, 



'tuning circuit' and 'frequency locking circuit' are described in section 4.2 



The tuning circuit provides the facility to sweep the transmitter frequency about a 
nominal centre frequency set by the computer. The sweep waveform and reflected power 
signal (detector DC output) are used to provide an oscilloscope display of reflected 
power against frequency. Typically, the display shows a dip in the reflected power 
corresponding to the resonator resonance. The resonator can be critically coupled by 
adjusting the coupling until the reflected power goes to zero on resonance. In the actual 
NMR experiment the tuning circuit has no function except to connect the reflected 
power signal to the lock-in amplifier. 

In all experimental configurations (frequency modulation, field modulation and 'no' 
modulation) we have used some frequency modulation. Where this is not used to 
detect NMR, only a small deviation (typically 0.5 MHz) is used. 
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Figure 4.1: Block diagram of the CW NMR spectrometer hardware (shown arranged 
for a frequency modulation experiment). 
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The frequency locking circuit is used to select 'tuning' mode or 'normal' mode. In 
normal mode the transmitter frequency is locked to the resonator frequency by nulling 
the first harmonic signal generated by the frequency modulation. This is done by using 
the output of the first harmonic lock-in (fsig) as an error signal which controls the 
microwave frequency via the oscillator control card (YIG-tuned oscillator control card, 



section 4.3.1). In tuning mode the frequency locking circuit simply connects the sweep 



signal from the tuning circuit to the sweep input of the oscillator control card. 

In the following two sections we outline the frequency and field modulation techniques 
used for this work. The arrangement 'without' modulation is essentially a simplification 
of the arrangement for either of these techniques: just the frequency locking circuit and 
DC detection. 



4.1.1 Frequency modulation technique 

Frequency modulation is easy to implement and has previously been used for CW NMR 
at Manchester by Ross |p7| . However, the earlier experiments used transmission cells 
with reasonably flat frequency response. In this work we adopted the unusual combi- 
nation of frequency modulation and a resonant cell. This gives the benefit of Q but 
problems with background from the resonator response (section 3.3. 2| ) . In particular, 



the experiment must be performed close to the resonant frequency of the resonator. 
This requirement becomes more stringent with increasing resonator Q. 

For maximum response the frequency deviation of the modulation should be comparable 
to line width of the NMR. However, there is no point in using a frequency deviation 
larger than the resonator line width because there will be little excitation when off 
resonance. This gives the rough condition that the modulation deviation should not 
be larger than either the resonator or NMR line widths. 



Figure |4.2| shows the two 'circuits' that are fundamental to this configuration. Firstly, 
the frequency locking loop uses a lock-in amplifier at the modulation frequency to tune 
the microwave transmitter so that the first harmonic signal is zero. Assuming that 
the resonator line shape is symmetric, at least over the width of the modulation, this 
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will lock the source to the centre of the resonance. This is very important with high Q 
resonators, particularly if the field sweep causes a shift in the resonator frequency. With 
a paramagnetic sample the dominant shift is due to the fall in differential susceptibility 
as the magnetization is saturated. Shifts of up to 20 MHz have been seen in H0F3 (0— 
8 T) . We have also found that the feedback loop can improve the frequency stability of 



the transmitter; see section 4.5.2. For the frequency locking loop to work properly, the 
frequency excursion should be smaller than the resonator line width; however we have 
already argued that this is compatible with efficient NMR excitation so no additional 
constraints are imposed. 

Secondly, there is the measurement circuit. The reference at n times the modulation 
frequency is used to detect the n th harmonic with a 2-phase lock-in amplifier. It is not 
strictly necessary to have a 2-phase lock-in amplifier, but the phase of the harmonic 
signal was found to vary significantly with time and the microwave frequency, and 
taking just the magnitude (rms of the 2-phases) gave the best results. This lock-in 
amplifier is equipped with an IEEE488 interface and readings are recorded by the 
computer. 

In operation, the spectrometer is set up with the transmitter frequency locked to the 
resonator resonance by the frequency locking circuit, and the applied field is swept 
through some range. As the field changes, the computer records the field, frequency, 
harmonic output and any other parameters requested (that can be monitored either di- 
rectly; with an additional IEEE488 equipped DMM; or using one of the spare analogue- 



to-digital converter inputs). The software is described in section 4.4 . Not shown in the 
schematic is the levelling loop with which the microwave power output of the source is 
kept constant. 
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Figure 4.2: Schematic of the experimental configuration for a frequency modulation 
CW NMR experiment. 
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4.1.2 Field modulation technique 



Field modulation may be produced using either the main solenoid or an additional field 
coil. In practice, it is not possible to modulate the main solenoid current at more than 
a few hertz because of its high inductance, ~10 H. Modulation of 1 mT at 10 Hz would 
produce ~6.6 V back-emf and the magnet power supply cuts out if the voltage across 
the magnet terminals exceeds 5 V. Using the main solenoid would also mean that all 
the metalwork of the resonator, support system and solenoid former would be within 
the modulation field and subject to eddy-current heating. An additional field coil has 
therefore been wound directly onto the resonator, with a much smaller area enclosed 



and hence a much smaller inductance. The field coil is described in section 4.5.1. It has 
an inductance of ~0.12 mH and the back-emf is not a significant problem. The factor 
limiting the modulation field and frequency is eddy-current heating of the resonator. 
Typically, modulations of 1-10 mT at ~100 Hz were attainable without unacceptable 
heating. 



Figure 4.5 shows the experimental configuration schematically. The microwave trans- 
mitter and resonator frequencies are not stable enough over time or with field for the 
transmitter to stay tuned to the resonance without feedback. The frequency locking 
loop described in the preceeding section was therefore used, but with a frequency ex- 
cursion small (typically ±0.5 MHz) compared to the resonator and NMR line widths. 
The frequency of this modulation was typically 100 kHz and no problems of interference 
with the low frequency field modulation were encountered. 

The field modulation coil is driven by a power amplifier designed to provide up to ±3 A 
into a low-impedance inductive load. The modulation signal and square-wave reference 
are generated by an audio frequency oscillator (Thandar type TG102) and the response 
measured using the 2-phase lock-in amplifier. The in-phase signal gave good derivative 
spectra indicating little phase shift in the system at the operating frequency. As in the 
frequency modulation experiments, the applied field was swept through some range 
with the computer storing the readings. 



67 



□ 
PC 

IEEE488 bus and Instrument interface 



I 



Reference 
oscillator 



Sine/square 
oscillator 



Phase 




Loop 
filter 


locked osc. 







Frequency locking loop 



FM 



Set frequency 



Microwave 
transmitter 



Frequency 
meter 



Power 
amplifier 



Field 
monitor 



to field modulation 
coils around resonator 



from magnet 
power supply 





2-phase 
lock-in 



signal 



ref 





Lock-in 




n 



ref 



signal 



Circulator Detector 



Cryostat containing 
resonator, sample and 
magnet 



Figure 4.3: Schematic of the experimental configuration for a field modulation CW 
NMR experiment. 
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4.2 Microwave system 



The microwave system is very simple: it consists of a transmitter, circulator, resonator 



containing the sample and a detector (figure 4.4). The quarter-wave coaxial resonators 
are described, in section |4.5| . Semi-rigid coaxial cable (0.141" diameter with PTFE 
dielectric and SMA connectors) is used for microwave connections in the source and 
from the cryostat top-plate down to the resonator. The other microwave connections 
use flexible coaxial cables (Suhner Ltd., 'Sucoflex'). 
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Figure 4.4: CW NMR configuration using a circulator and resonator. 

This section also includes descriptions of the 'tuning circuit' and the 'frequency locking 
circuit' which are not part of the microwave system per se. They are, however, very 
important in tuning the system. 

4.2.1 Microwave transmitter 



Figure 4.5 shows schematically the arrangement of microwave components comprising 
the transmitter. The oscillator output is divided between two arms by a 3 dB split- 
ter. Half of the power is used to provide an output for frequency measurement. The 
additional -10 dB directional coupler in this arm was used to provide an output for 
a prescaler that has not been implemented. The other arm uses an electronic atten- 
uator (Hewlet Packard type 33142A) and an amplifier (described below) to provide 
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variable output power from <0.1 mW to ~130 mW. The directional detector (Krytar 
type 1211S, 1-12.4 GHz) allows the output power to be monitored. Usually the output 
power is levelled, keeping the directional detector output (level out) constant by feed- 
back to the level control input. However, the levelling loop control card (section |4.3.2 ) 
provides the facility of levelling the power at another detector. For example, it might be 
best to level the transmitted power in an experiment using a non-resonant transmission 
cell. 
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Figure 4.5: Schematic of the microwave transmitter. 

The oscillator (Philips/Sivers Lab. type PM 7020C) is based on a transistor coupled 
to an yttrium-iron-garnet (YIG) resonator. The YIG resonator has high Q and is 
tuned by changing the applied field. The oscillator is tunable over the range 4-8 GHz 
with power output varying between 22 and 39 mW over that range (measured). Other 
devices are available to cover the range 2-8 GHz. In addition to the tuning coil, the 
oscillator has another small coil to allow frequency modulation (FM). This provides up 
to ±50 MHz deviation with a bandwidth from DC to 100 kHz (-3 dB). 

The noise output of YIG-tuned oscillators is low, typically 50 dB down on the fun- 
demental; however there is a much larger harmonic content. The oscillator used is 
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specified to have total harmonic content more than 12 dB down on the fundemental. 
The current system operates over just one octave (4-8 GHz) and harmonics will be 
outside the bandwidth of the amplifier and considerably attenuated. Microwave noise 
has not been a significant problem. However, harmonics could be more of a problem 
in a 2-8 GHz system where, between 2 and 2.7 GHz, both second and third harmonics 
would be within the system bandwidth. 

The amplifier (Avantek Inc. type AWT-8052) is a 2-8 GHz device providing «17 dB 
gain with a maximum output power of ~130 mW (21 dBm). Modifying the source to 
cover 2-8 GHz would simply require replacement of the oscillator. 



4.2.2 Circulator 



Circulators are available covering octave bands and the spectrometer uses a 4-8 GHz 
device. Extension of the spectrometer operating range to cover 2-8 GHz would require 
two circulators which could either be manually exchanged or switched in and out using 
microwave relays. Ideally, circulators should not be placed in magnetic fields but no 
problems have been experienced with the circulator only 1.5 m from the 8 T magnet 
(roughly on axis). 



4.2.3 Detector 



The microwave detector uses a GaAs 'planar-doped barrier' diode matched to the 50 
transmission line by a miniature thin film circuit. Details are given below: 

Hewlet Packard, type HP33334C. 

Frequency range 0.01-33 GHz, response flat to ±0.45 dB over 0.01-26.5 GHz. 
Maximum continuous input power 200 mW (23 dBm). 



The detector and a pre-amplifier are built into a shielded box (figure 4.6). Separately 



buffered AC and DC coupled outputs are provided. The DC coupled output is used 
for tuning and to monitor operation of the frequency locking circuit. The AC coupled 
output is connected to the 2-phase lock-in amplifier. 
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Figure 4.6: Schematic of the microwave detector and pre-amplifier. 
4.2.4 Tuning circuit 

This circuit is required to set up the spectrometer prior to NMR measurement, and is 
shown schematically in figure \i.7[ The computer sets the nominal microwave frequency 
and in 'tune' mode the tuning circuit provides a symmetric frequency sweep about this 
frequency. A calibrated dial allows the sweep range to be adjusted between and 
1 GHz (±500 MHz). The YIG sweep output controls the microwave frequency via 
external sweep input of the oscillator control card. A second output provides a ±1 V 
triangular voltage waveform to control the oscilloscope X-deflection. The detector DC 
output (reflected microwave power) is buffered to provide the Y-deflection, and there 
is a straight-through connection to the first harmonic lock-in. 

The usual method of operation is initially to set the sweep to a wide frequency range 
to find the resonator resonance. When the resonance is located, the coupling can be 
adjusted roughly for low reflected power on resonance. Then the resonator can be 
tuned to the desired frequency and the sweep range gradually turned down for greater 
resolution. This facility is useful for checking the resonator Q and for fine adjustment of 
the coupling. Once the source is centred on the resonance the frequency locking circuit 
can be switched in. Without adjusting the oscilloscope, the state of the frequency 
locking loop is shown by the reflected power trace; see the following section. The 
X-deflection provides a convenient timebase without having to adjust the oscilloscope. 
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Figure 4.7: Schematic of the tuning circuit. 
4.2.5 Frequency locking circuit 

This circuit uses the first harmonic signal from the lock-in amplifier as an error signal 



which is amplified, filtered and fed back to control the microwave frequency. Figure 4.5 
shows the circuit schematically. The principle of operation is that the first harmonic 
in-phase signal will approximate to the derivative of the resonator line shape. This 
will be zero at the centre of the resonance, so using feedback to adjust the microwave 
frequency until the first harmonic signal is zero will lock the system to the resonator 
frequency. 

When the microwave frequency is within the resonator line width the frequency locking 
circuit can be switched in (normal). We have found it best to start with the loop gain 
turned right down (~0) and then slowly increase the gain (slowly with respect to the 
time constant of the lock-in). Usually the loop will lock at some point, indicated by 
the lock indicator and by a steady trace on the oscilloscope. Further increase of the 
gain usually improves the stability up to a point where the microwave frequency starts 
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to oscillate. The gain can then be turned down a little to stabilise the loop. If the 
frequency 'jumps' out of the resonance the procedure can be repeated by turning the 
loop gain down again. 

It is not practicable to predict the gain required for stability of the loop. The loop 
gain is affected by the resonator line width, the microwave power, the coupling to the 
resonator and the gain and time constant of the lock-in amplifier. The following 'typical' 
values were taken for the gains of each part of the loop: YIG-tuned oscillator tuning 
sensitivity « 4x10 s Hz/V, resonator /detector /lock- in frequency to voltage conversion 
factor « 4x 10~ 7 V/Hz (resonator Q « 800, lock-in gain adjusted to give ±5 V peak first 
harmonic output) and frequency locking circuit amplifier gain ~ 200. The component 
values were chosen: Rl = 10 kf2, R2 = 470 Q and C = 1 fiF, giving the loop a unity- 
gain frequency of f«6 kHz with the filter 'zero' at «300 Hz for a safe phase-margin 
(with 30 ms lock-in time constant). 
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Figure 4.8: Schematic of the frequency locking circuit. 
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4.3 Interface unit 



We use the term 'interface unit' to describe the unit that contains all the computer 
controlled interfacing and control electronics built for the spectrometer. The individual 
parts of the interface unit are build as 'cards'; grouped in one unit because this allows 
the computer interface and power supplies to be shared. First, we describe the shared 
parts of interface unit, followed by descriptions of each card that it contains. 

The interface unit is connected to the control computer via the home-built instrument 
interface. It would be preferable to connect the interface unit to the IEEE488 bus, along 
with the other instruments. However, the IEEE488 bus is sufficiently complicated that 
an interface to it would have to be controlled by a microprocessor and hence take a 
considerable amount of time to design and build. Such an interface could be retro-fitted 
in place of the current interface. 

The instrument interface is controlled from the computer via a PIO card (based on the 
8255 integrated circuit). Another circuit provides signal conditioning and buffering 
for the inputs and outputs of the PIO card. Buffering was implemented to allow the 
interface to operate over a long cable (10 m at present; could be longer if required) 
so that computer could be physically separated from the experiment if interference 
was a problem. The signal conditioning uses simple RC filters to slow the transitions 
to « 1 fis, to reduce crosstalk and radiated interference. In practice, no interference 
problems were experienced with the PC about 3 m from the cryostat and microwave 
components. 

The instrument interface is a simple parallel interface. There are 8 input lines, 8 output 
lines and 8 control signals (3 outputs, 5 inputs; directions given with respect to the 
computer). The 8 output lines use are used to specify either a 6 bit address and the 
direction of transfer, or 7 bits of output data. Two other lines are connected to allow 
the software to tell if the instrument is connected and if so, if the power is switched on. 
The software can thus respond with sensible error messages in these circumstances. 
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Inside the interface unit, the instrument interface is connected to the instrument bus. 
Physically, this is implemented as a backplane with wire-wrapped DIN41612 indirect 
edge-connectors. These provide 60 pins and 4 coaxial connections for each card. Mul- 
tiple pins are used for separate digital and anolgue power supply connections, 24 pins 
are used for the bus and 12 are unused. The addresses are decoded by individual cards; 
the bus connections are the same to all cards. However, the positions of the cards 
are fixed by point-to-point coaxial connections made on the backplane. There are also 
special power supply connections for the YIG-oscillator control card. 

Cards are built on either prototyping boards or specially designed printed circuit 



boards. Figure 4.9 shows a view of the instrument box with the lid off. The in- 
strument interface is in the back-left corner and power supplies are to the back-right. 
The instrument bus runs across the centre (blue wires) with the individual cards to the 
front. 




Figure 4.9: Photograph of the interface unit with lid removed showing internal organ- 
isation and contruction (the unit is ~480 mm wide). 
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4.3.1 YIG-tuned oscillator control card 



The YIG-tuned oscillator in the microwave transmitter has stringent power supply 
requirements which have been carefully adhered to. The power supplies have been kept 
separate from the supplies to the other cards. Operational amplifier based regulators 
are used because monolithic regulator integrated circuits are not stable enough. The 
regulators are located on heatsinks at the back of the intrument box with the control 
card providing connection to the microwave transmitter via the same cable as the 
tuning and modulation signals. 



DAC 



Buffer 



INT SWEEP 



Buffer 

EXT SWEEP ^ |X 



Buffer 



INT FM 



EXT FM 



Power supply 




Bus 

interface 



Buffer 



Instrument 
bus | 
I 

1 

I 



to YIG-oscillator 
^. TUNE 



High-power 
Gain 10 Buffer 



100 £2 




FM 



Power 



Figure 4.10: Schematic of the YIG-tuned oscillator control card 

A 14 bit digital-to-analogue converter (DAC) provides the tuning signal for the oscil- 
lator, figure 4.10| . This gives a frequency resolution of better than 300 kHz over the 



4-8 GHz tuning range. The tuning linearity is limited by the oscillator and not the 
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DAC. Linearity is not very important because the oscillator tuning can be adjusted 
until the measured frequency is as desired. The DAC output voltage and two sweep 
input voltages are summed and buffered to provide the tuning voltage. The two sweep 
inputs can be enabled or disabled from software using the analogue switches. 

The frequency modulation coil in the oscillator gives a frequency deviation of ±50 MHz 
for ±150 mA coil drive (the maximum allowed). The coil has an inductance of <2 /iH 
and a resistance of <1 £1. By putting 100 f2 in series with the coil the combined 
impedance is almost resistive over the frequency range 0-100 kHz. At 100 kHz the 
phase shift is ~0.7°. This simplifies the drive amplifier and an EL2001 buffer IC 
was used (Elantec Ltd. available from RS Ltd.). An output swing of ±10 V provides 
±33 MHz deviation and, with the buffer amplifier operating from ±15 V power supplies, 
it is impossible to over-drive the modulation coil. The buffer amplifier is preceeded by 



a gain 10 amplifier and FM input selection switches as shown in figure 4.10. The 
measured performance is: ±33 MHz deviation for ±1 V input (saturation at ±1.2 V) 
and frequency response from DC to >200 kHz (5% drop in deviation relative to DC). 

The sweep inputs are used for the frequency locking loop and for tuning the resonator. 
The tuning circuit provides a frequency sweep adjustable up to ±500 MHz. The FM 
inputs are used for frequency modulation, driven from either the modulation card or 
an external oscillator. There are internal (on the backplane) and external (on the front 
panel) connections for both the FM and sweep inputs. 



4.3.2 Level control card 



The microwave power was levelled at the output of the microwave source for the work 
reported in this thesis. To do this, output of the directional detector is compared with 
a reference voltage and the error signal is used to drive the electronic attenuator in the 
transmitter, completing the feedback loop. For this particular arrangement, it would 
be best to have the levelling loop built into the transmitter, with only the reference set 
by the computer. However, in the interests of flexibility, the error amplifier reference 
and switching circuits have been built as a card in the interface unit. 
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The level control card is shown schematically in figure 4.11 . The input is selected 
from one of three connections on the front panel, or one connection on the backplane. 
The high-gain input (I NO) is for the directional detector in the transmitter (output 
0-250 mV) and the medium-gain input (INI) is for a detector receiving most of the 
oscillator output, for example the detector of a transmission cell experiment (typically 
0-1 V). The other two inputs, IN2 and IN3 provide auxiliary inputs on the front panel 
and backplane respectively. The software selected input is compared with a reference 
voltage set by DACO. The amplified error voltage is filtered by a simple low-pass filter 
(r ~ 0.4 s) which averages variations in amplitude caused by frequency modulation (in 
the range 1-100 kHz). To level the power at the modulation frequency would require 
a fast levelling loop with appropriate phase compensation. The arrangement described 
here is not suitable for that purpose. A detection circuit senses whether the error 
amplifier is saturated, indicating unlevelled output. The detection circuit has an LED 
indicator and the output can also be read by software. 

The output to control the electronic attenuator in the microwave transmitter (via the 
driver circuit) is selected from software. Either the filtered error-amplifier output, or a 
DC control voltage from a second DAC (DAC1) can be selected. Direct control using 
DAC1 allows unlevelled operation. 



4.3.3 Modulation oscillator and detection reference cards 

For harmonic detection, a sinusoidal modulation signal and a reference square-wave 
at the harmonic frequency are required. The reference at frequency nf must have a 
stable phase with respect to the modulation sine-wave at frequency /. To do this a 
crystal oscillator provides a 10 MHz digital reference which is divided by 2m to give 
frequency nf (where m = ^j)- The nf signal is then divided by n to give a reference 
at frequency /. 

The modulation sinusoid is produced by phase-locking an oscillator to the digital refer- 
ence at frequency /. The phase-locked loop (PLL) has a fairly long time constant which 
gives it good noise immunity; the phase jitter relative to the reference is < 1°. However, 
when the reference frequency is changed it may take several seconds for the PLL to 
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Figure 4.11: Schematic of the level control card. 

capture and lock onto the new frequency. There are two frequency ranges for the PLL, 
selected by an electronically switched capacitance in the oscillator. Once locked, the 
PLL is stable over the ranges 800 Hz to 19 kHz on the low range, and 6 kHz to 140 kHz 
on the high range. The output amplitude (0-1 V) is controlled by a calibrated dial. 
When used to provide frequency modulation via the YIG-tuned oscillator control card 
this corresponds to 0-33.3 MHz peak deviation (0-66.7 MHz modulation width). 

The digital reference generation and the analogue/digital PLL are implemented on 
separate cards. The reference generator is controlled by software to set m in the range 
1-2048 giving a frequency range of 2.44 kHz to 5 MHz. The resolution goes as 1/m so 
the frequencies available are in < 2Hz steps at «2.5 kHz but in 2 kHz steps at 100 kHz. 
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Conversion of frequency to division ratio is done in software, selecting the available 
frequency closest to that requested. The harmonic n may be selected from the range 
1-32. However, the lowest PLL frequesncy is 800 Hz which allows only third harmonic 
detection with the lowest detection frequency of 2.44 kHz. 



4.3.4 Magnet control card 

The magnet control card not only controls the current through the magnet solenoid but 
also allows it to be monitored. This is achieved by amplification and 14 bit analogue- 
to-digital conversion of the voltage developed across a series resistor in the magnet 
power supply (see section fO|). The field is determined using the calibration factor 
supplied with the magnet (0.0957 T/A). The amplifier and ADC were calibrated by 
simultaneously recording the ADC reading and the voltage across the series resistor for 
a range of currents (using a Keithley type 197 DMM ±0.02 % accuracy). 

The magnet current may be controlled either manually, using the controls on the mag- 
net power supply, or from software via the magnet control card. A 12 bit digital- 
to-analogue converter sets the reference for the magnet power supply. The rate of 
stepping can be kept within safe limits by monitoring the back-emf developed across 
the solenoid from software. In practice, it was found more convenient to control the 
magnet manually, just using the computer interface to monitor and record the field. 



4.3.5 Analogue-to-digital conversion card 

A schematic diagram of the analogue-to-digital conversion (ADC) card is shown in 



figure 4.12. With ±10 V input range the 14 bit ADC provides a resolution of ~1.2 mV. 
Testing has shown that the accuracy approaches the resolution (a few mV), but noise 
has proved more of a problem. Noise on the readings is about ±4 counts giving a 
single-reading accuracy equivalent to 11-12 bits. This has been overcome by taking 
multiple readings in software whenever greater accuracy is required. For example, when 
measuring the solenoid current to find the applied magnetic field, the software averages 
10 samples reducing the noise to 1 or 2 counts. An amplifier with electronically switched 
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gain setting resistors allows the input range to be set to either ±1 V or ±10 V from 
software. Although there is only one ADC, one of eight input signals can be selected 
by software control of the 8-to-l input multiplexer. Four of the inputs are connected 
to front-panel connectors, the other four are used by other cards via point-to-point 
coaxial connections on the backplane. Differential inputs avoid the need to connect the 
shield to the instrument earth, so preventing the formation of earth loops. The inputs 
are protected from out-of-range signals of up to ±25 V. 
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Figure 4.12: Schematic of the ADC card. 
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4.4 Software 



The current version of the control software (NMR1-3) involves over 250 kbytes of source 
code and is not listed here. However, the structure and operation of the software 
is outlined below. It is written in C using Borland Turbo C to run on an IBM PC 
compatible under MSDOS (v5.0). Interface routines for the IEEE488 bus (supplied by 
Brain Boxes Ltd.) allow file-mode control from C. Interface routines to communicate 
with the instrument over the home-built instrument interface were also written in 
C, without recourse to assembly language. The structure of the software is shown 



schematically in figure 4.13 



When deciding how the control software should work, several principles were used for 
guidance: 



• Conversational (keyboard based) style of operation as opposed to menus. Where 
possible commands should be close to English, and unambiguous abbreviations 
should be accepted. For example, the command to start a sweep is SWEEP, followed 
by the parameter to sweep, the starting value, the ending value and the step size. 
SW, SWE and SWEE are allowed. 

• Data should be saved automatically. The file names are generated automatically 
using a prefix and a sequential index. If the data are not required then the file / s 
may be deleted later by intentional user action. 

• Data should be displayed graphically as they are collected. The software allows 
axes to be scaled, data to be binned in groups of n samples, and data from 
different parameters to displayed without interrupting operation. 
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Figure 4.13: Schematic showing the structure of the spectrometer control software. 
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4.5 Resonators 



A variety of quarter- wave coaxial resonators have been used with the Manchester pulsed 
spectrometer for many years. The designs used for this work have been developed from 



those of Carboni [JT^, 20], McMorrow [51] and Beck H. Previous designs for pulsed 
NMR have been designed with Q's of <500 to limit the ring time. Although fundemen- 
tally of similar design, the resonators for CW NMR were more carefully constructed in 
order to raise the Q. The resulting resonator Q's ranged from ~700 to 2500 depending 
on precise configuration and care of assembly. Brass was used for all the resonators. 
It would be possible to reduce losses and hence increase the Q by using copper, but 
copper is much more difficult to machine than brass. Probably a better approach would 
be to gold plate a brass resonator. 



Figure |4.14| shows the tunable resonator used for the most of this work. One design 
feature instrumental in raising the Q is the slight rake on the ends of the central barrel 
section (~ 1°, shown exaggerated). The top- loading version (shown) has significantly 
lower Q than one where there is only the coupling loop entering the resonator and the 
sample is glued directly to the central conductor (with Stycast 1266), typically 700- 



1000 instead of up to 2500. Carboni [19] combined the ideals of small effective resonator 



volume and wide tuning range by using a tapered central conductor. The sample can be 
placed close to a central conductor of small radius giving the small effective resonator 
volume, whilst maintaining a large end area of the central conductor for wide tuning 
range. 

A similar resonator was used for the field modulation experiments and is shown in 



figure 4.15. The internal dimensions are the same but the wall of the barrel section 



was made much thinner and a superconducting modulation coil wound around it. Fig- 



ure 4.16 is a photograph of this resonator complete with brass base and support rods. 
The rod to the left of the photograph, with a brass coupling section, moves the tuning 
piston via a screw drive and a lever in the base. The stainless steel pipes brazed into 
the top of the resonator are guide pipes for the microwave coax and for the sample rod 
(one can be seen in the photograph, the other is obscured by the support rod in the 
foreground). All the stainless steel tubes are non-magnetic. 
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Figure 4.14: Schematic section through the tunable quarter- wave coaxial resonator 
used for frequency modulated CW NMR. The design shown allows top-loading of the 
sample, alternatively the sample may be glued directly to the central conductor. 

4.5.1 Field modulation coil 

In designing the modulation coil, eddy currents induced in the brass former /resonator 
were the main consideration. Eddy currents not only result in heating, they also reduce 
the field inside the resonator. For effective field modulation it was estimated that 
~1-10 mT amplitude would be required at the sample. 
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Figure 4.15: Schematic section through the tunable quarter- wave coaxial resonator 
used for field modulated CW NMR. The design shown allows top-loading of the sample, 
alternatively the sample may be glued directly to the central conductor. 
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Figure 4.16: Photograph of the resonator with field modulation coil. The brass base 
section is «70 mm in diameter. 

The strategy adopted was to minimise the physical size of the modulation coil by 
winding it directly onto a thin- walled brass resonator. The dimensions of the resonator 
are determined by the resonance frequency required and by the need to allow access 
for the coupling loop and sample rod. The internal dimensions were chosen as 12 mm 
diameter by 10 mm long, the same as for the resonator without a modulation coil, 
but with 0.5 mm wall thickness. Taking the resistivity of brass as ~ 4 x 10 -8 Sim at 
4.2 K, and assuming a ~ 10 mm length of the resonator wall to be linked by all the 
flux through the the coil, gives an expression for the average power dissipation due to 
eddy currents: 



where P e ddy is the average power dissipation (W); v mo d is the modulation frequency 




(4.1) 
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(Hz); and (B moc [) is the RMS of the modulation field (T). This is probably an underes- 
timate because it does not allow for other metalwork in the vacinity, in particular the 
top and base-plate of the resonator. However, using equation ^1] as an estimate, 10 mT 
RMS modulation at 10 Hz (or 1 mT at 100 Hz) gives ~ 20 fiW power dissipation which 
is not significant. For 10 mT at 100 Hz the dissipation is ~ 2 mW which is significant, 
but not prohibitive. We calculated that a 5 mm long coil of 100 turns wound directly 
onto the former would produce ~ 6.5 mT/A at the sample. 

The actual coil consists of 111 turns of 0.38 mm diameter 'Niomax' (TiNb in copper) 
wire. The wire has a critical field of > 8 T, an important consideration as it is inside 
the bore of the main magnet. The estimated inductance of the coil was 230 /^H without 
the former. Using this inductance, a modulation frequency of 100 Hz and a modulation 
current of 1 A rms the back emf is ~0.15 V. A simple amplifier, using an operation 
amplifier with a complimentary darlington output stage, was built to provide up to 
±3 A drive for the modulation coil. 
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Figure 4.17: Equivalent circuit used to model the field modulation coil and induced 
currents in the former. 



From the equivalent circuit shown in figure |4.17| , we may write down an expression for 
the effective impedance of the modulation coil, 
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„ „ uj 2 M 2 R e ( T co 2 M 2 L e \ 

Zc = Rc + m + um +W [ L *- m + uj 2 l 2 ' ' (42) 



where i? c , L c , R e and L e are the resistances and inductances of the coil (drive) and 
eddy circuits respectively. The coupling between the circuits is modelled by the mutual 
inductance M. The values of R e , L e and M are not independent in terms of their effect 



on the impedance of the drive circuit. In order to parametrize equation 4.2 we make 
the substitutions 

m = —— and I = — -. (4.3) 

Re Re 

The measured value of R c , including the copper leads, was 66 mf2 with the resonator 
at 4.2 K. The coil itself was superconducting so this resistance is essentially that of the 
leads from the cryostat top-plate. Using the measured value of R c , the parameters L c , 
m and I were determined by fitting the magnitude and phase of Z c to experimental 
data from 20 Hz to 2 kHz. The fits are shown in figure fOJ: L c = 1.24xl0 -4 H, 



m = 4.63 xl0~ 8 B 2 /n and I = 3.98 xl0~ 4 H/fi. The coupling k M between the drive 
and eddy cicuits is given by 

, L c L e L c l 

= ~~M 2 ~ = ~m 2 ' (44) 

The parameters L c , m and I are not very accurately determined from the fits; small 
changes in the fitting method or in the data points produce significantly different re- 
sults. However, ku does appear to be reliably determined from the fits. From all fits, 
where the parameters varied by as much as 30 %, the result Um ~ 0.94 was obtained 
with <1% variation. This confirms that the frequency modulation coil is strongly 
coupled to the eddy circuit, as expected for a coil wound onto a brass resonator. 

The resonator for field modulation is a prototype. It should be possible to reduce the 
power dissipation significantly by thinning the barrel section further and reducing the 
amount of other metal in the vicinity. A better approach would be to use an insulating 
former, plated with gold on the inside. Also, the leads from the top of the cryostat to 
the coil 'rattled' at high drive currents. This could be simply fixed by taping the wires 
to one of the support rods. 
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Figure 4.18: Phase and magnitude components of the modulation coil impedance. 
The points represent measured values and the line is a fit to the data. 
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4.5.2 Resonator frequency stability 



At temperatures where the helium bath is superfluid (T\ ~ 2.19 K) the frequency 
stability of the resonator is very good. Using the frequency locking loop to lock the 
oscillator frequency onto the resonance significantly improves on the stability of the 
oscillator alone. Typically, the frequency jitter is <20 kHz (~0.001 %) over a period 
of seconds. This is increased with large FM deviation (>5 MHz) but still remains 
<300 kHz. There is, however, a much larger long term drift, up to to ~1 MHz per 
hour, and sometimes the frequency 'jumps', typically by 100-200 kHz. The drift is 
attributed to slow movement of the tuning components and the jumps are probably 
caused by knocking the apparatus. Jumps are more frequent when moving around close 
to the spectrometer than when it is left alone. 

The formation of helium bubbles in the resonator was the most significant problem 
when working at 4.2 K. Various attempts were made to alleviate it: drilling holes in 
the top of the resonator to allow bubbles to escape, drilling holes in both the top and 
the bottom of the resonator, and working at low power. Drilling holes had very little 
effect, presumably because the resonator already has one or two fairly large holes in the 
top where the coaxial cable and sample rod fit loosly. Reducing the microwave power 
slows down bubble formation and reduces the frequency excursion (smaller bubbles) 



as shown in figure 4.1S , The field modulation coil also heats the helium inside the 
resonator and causes similar behaviour. Modulation coil drive of 2 A rms and 0.1 mW 
microwave input produces bubbling similar to 0.5 mW microwave input alone (no field 
modulation) . 

The result of this bubbling was that it proved impossible to work with high power at 
4.2 K. Typically, microwave powers of 1 mW or less were used. The bubbling disappears 
as the helium bath is cooled through the superfluid transition. Presumably there would 
be no bubbling problem in a continuous flow system where there is helium gas in the 
resonator. A similar design of resonator has worked with high stability when evacuated 
at nitrogen temperatures (on a test experiment using EPR of DPPH). 
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Figure 4.19: Graphs of resonator resonance frequency variation with time for various 
microwave input powers. The lower power traces clearly show the resonance frequency 
rise as a bubble forms, and then fall suddenly as it escapes. 
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4.6 Cryogenics and magnet 



Samples in resonant cells are suspended in the bore of a superconducting magnet at 
liquid helium temperatures. The magnet, magnet power supply (PS120A) and cryostat 
were supplied by Cryogenic Consultants Ltd. At 4.2 K the maximum field available is 
8 T. By pumping on the helium bath containing both the magnet and the experimental 
resonator the temperature could be reduced to «1. 3 K with available pumps. Reducing 
the temperature to 1.5 K or less allows fields of up to 10 T to attained. Cooling the 
resonator and sample uses a significant proportion of the helium in the cryostat because 
the magnet is also cooled. The cryostat holds «21 1 of helium, of which «2.5 1 are 
required to cover the magnet (essential to avoid quenching). Cooling from 4.2 K down 
to 1.5 K uses «10 1 and takes 2-3 hours. Fortunately, when cooled, the helium boil-off 
rate is not significantly greater than when at 4.2 K (0.3 1 per hour with no current in 
the magnet leads and no microwave input). The magnet is equipped with a persistent- 
mode switch that allows a constant field to be maintained without any current in the 
magnet leads. At high fields this gives a significant reduction in the boil-off rate. 

To find the sample temperature, the pressure above the helium bath is measured using 
a McLeod gauge. The temperature error is ±0.05 K at 1.5 K, provided that localised 
sample heating is not significant. The field is determined by measuring the current 
through the magnet soleniod and using a field-current calibration supplied by the man- 
ufacturer. The field is homogeneous to within 0.18% in a cylinder of diameter 10 mm 
and height 10 mm at the centre of the solenoid. 
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4.7 Techniques 



In this section we summarize experience of the successful field sweep technique, the 
unsuccessful frequency sweep technique, and results obtained 'without' modulation. 
Representative spectra from the field sweep technique can be found in chapter 5; none 
are given here. 



4.7.1 Field sweeps 

The field sweep technique works very well, principally because the frequency response 
of the microwave system, excluding the sample, is unaffected by the applied field. The 
sample may cause a shift in the resonant frequency of the resonator (section 4.7.3| ) but 



this is a relatively small and 'smooth' effect, and has not proved to be a problem. 

We have obtained good NMR spectra from both field modulation and frequency mod- 
ulation experiments. The principal advantage of field modulation is that the spectra 
have straightforward derivative character and it is thus easy to identify the centres 
of the NMR lines. Frequency modulation with third (or higher) harmonic detection 
gives more confusing line shapes. In some cases, it is the uncertainty in identifying the 
central feature of the line shape that dominates the measurement uncertainty. As ex- 
pected, we found that no advantage was gained from using higher than third harmonic 
detection. 

In order to sweep the field, the current in the superconducting magnet must be ramped 
up and down. This means that it is not possible to use the magnet in persistent mode, 
which imposes a penalty in terms of helium boil-off. 



4.7.2 Frequency sweeps 

Frequency swept systems require either a broadband transmission cell with the most 
level frequency response attainable; or a resonant cell whose resonant frequency is made 
to track that of the source, or conversely. The former method has been implemented by 
Cha [24| and, in this department, by Ross 67]. The latter strategy has been investigated 
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as part of this work. 

We have arranged the spectrometer so that the frequency of the microwave transmitter 
tracks the resonant frequency of the cell. Thus frequency sweeps can be performed by 
manual or motorized tuning the cell. Manual tuning causes far too much disturbance to 
the system, giving a background signal that is large, discontinuous and un-repeatable. 
Motorized tuning gives a much more repeatable background. In all attempts the back- 
ground was much larger than any expected NMR signal. However, we note that the 
background was significantly reduced with field as opposed to frequency modulation. 
We conclude that this technique is inappropiate for systems where the transfer charac- 
teristics depend strongly on frequency. 



4.7.3 Shift in the resonator frequency due to NMR 

As part of this work we have detected NMR by the shift in the frequency of the 
combined resonator and sample resonance. This unusual detection method was feasible 
only because of the exceptionally large enhanced nuclear susceptibility of H0F3. 



In section |3 . 3 . 7| , we have shown that in the regime where the resonator line width 
is much smaller than the NMR line width, the combined resonance frequency will 
respond to changes in the sample susceptibility. To a very good approximation, we 
expect a shift in frequency is proportional to — x' s , the dispersive part of the transverse 
RF susceptibility. The shift does not depend on the resonator Q but the frequency 
stability of the system is improved by high Q so, overall, the technique benfits from 
high Q. 

When performing initial field sweep experiments on H0F3 the shift in resonator res- 
onance frequency resulting from the field dependece of the differential electronic sus- 
ceptibility was very large, ~ 20 MHz. This suggested that it might be possible to see 



NMR by monitoring the frequency of the 'locked' system. Figure 4.20 shows a partially 
resolved NMR spectrum obtained in this way. The data were taken at 4.2 K using a 
single crystal sample with the a axis parallel to the field (4.2 x 1.9 x 1.6 mm). The 
resonator Q was ~ 2500 and the frequency jitter was less than 100 kHz. 
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We have fitted a 7-line holmium spectrum to the frequency shift data: the result is 



shown as the solid line in figure 4.20. The model has 7 parameters: the position of the 
central line, the quadrupole splitting (in field), the resonance amplitude (the same for 
all 7 lines), the line width (the same for all 7 lines), and 3 parameters to allow for a 
quadratic background. The fit does not reproduce the structure that can be seen in 
the experimental data. By constraining the line width in the model, it is possible to 
reproduce a similar structure to that of the experimental data. However, we have been 
unable to get accurate agreement. This is not unreasonable because we have calculated 
the transverse RF susceptibilty from the Bloch equations which assume a Lorentzian 
line shape. In H0F3 the holmium resonances are inhomogeneously broadened and the 
resulting line shape is not Lorentzian. 



Figure 4.21 illustrates how the simple model of 7 resonances can produce a structure 
similar to that of the experimental spectrum. The figure shows 7 individual dispersion 
curves and the sum which represents the expected frequency shift. The sum clearly 
shows structure similar to that the experimental frequency shift data. 



Figure 4.20 also showns the DC reflected power as a function of the applied field. There 
is a significant change in the reflected power resulting from NMR absorption and some 
structure is evident. However, the overall shape is rather asymmetric and no attempt 
has been made to model it. 
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Figure 4.20: NMR of holmium in H0F3 seen by frequency shift of the resonator 
resonance (upper) and by DC power absorption (lower). 
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Figure 4.21: Form of resonator frequency shift predicted from the susceptibility change 
caused by holmium NMR. The 7 components resulting from the individual NMR lines 
are shown in conjuction with the resultant which is simply the sum of the component 
shifts. 
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5. Holmium Trifluoride 



This chapter is concerned with the low temperature magnetic and spectroscopic prop- 
erties of H0F3. The magnetic properties are dominated by the holmium ions and we 
use the Hamiltonian described in chapter 2 to model the ground J-manifold. H0F3 



orders antiferromagnetically at Tjy = 0.53 K [13] but here we consider only the param- 
agnetic phase. Although the term Van Vleck paramagnet is often applied to H0F3, the 
susceptibility is not temperature independent [11]. This is because of the large nuclear 



contribution and the need to reduce the temperature below ~2 K before thermal pop- 
ulation effects become insignificant. Therefore it is probably more helpful to consider 
H0F3 simply as a singlet electronic ground state system. In addition to magnetometry 
and NMR measurements on H0F3 , we have also made NMR measurements on holmium 
as a dilute substituent in YF3: l%Ho:YF3. 

The first part of this chapter surveys the crystal structure determinations for H0F3 
and YF3 in the literature. From the structure some important properties that affect 
magnetic measurements are deduced. When discussing the structure we refer to the 
'holmium sites' throughout. In the dilute material these will be mostly yttrium sites 
with occasional holmium occupancy. 

Measurements of the electronic energy levels of the Ho 3+ ions from several sources are 
compared with one another. Sharma et al [[H]] have determined crystal-field parameters 
from their measurements of the electronic energy levels, and we use these as the basis 
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of our calculations. 



We have studied the magnetization of H0F3 for applied fields in the range 0-10 T along 
the principal crystallographic directions. The results show highly anisotropic magneti- 



zation, which accords with the earlier measurements of Bleaney et al [11]. The magne- 
tization data are compared with calculations from the crystal-field parameters. Both 
pulsed and CW NMR have also been used to study H0F3 and l%Ho:YF3. Whereas the 
magnetization measurements give thermal and spatial averages, the NMR frequencies 
reflect the properties of the ground state of the holmium ions. Thus, NMR provides a 
complementary test of the crystal-field parameters. We have also used the crystal-field 
parameters to calculate some properties that have not been measured, and suggest 
some potentially interesting experiments. 

The samples studied in this work were prepared at the Clarendon Laboratory, Oxford. 



Methods of preparing rare-earth trifluorides are reviewed by Carlson and Schmidt 22]. 
The H0F3 crystals are pale pink while the l%Ho:YF3 crystals are clear. Although 
we refer to the Ho in YF3 crystals as l%Ho, the concentration is nominal and not 
important to this work. We only require the holmium ions to be 'dilute' so that the 
dipolar and exchange fields are insignificant. 



5.1 Crystal structure 



The structure of rare-earth trifluorides was first studied by Zalkin and Templeton 78] 
using x-ray powder diffraction. They found TbF3, DyF3, H0F3, ErF3, TmF3 and YbF3 
to be isostructural with YF3. The unit cell is orthorhombic, space group 
and contains four formula units. Following convention, the unit cell axes are labelled 



a, b and c where the lengths are related by b > a > c. Table 5.1 gives the unit cell 



dimensions of YF3, H0F3 and TbF3 JT8|. The dimensions for H0F3 and TbF3 were 



verified by Piotrowski et al [57] 



Zalkin and Templeton [78] also found some of the trifluorides with the hexagonal LaF3 



structure but this is not confirmed in later work on H0F3 and TbF3 by Piotrowski 



et al [p7f[. Gries and Haschke [38] compare data from trifluorides with and without 
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Compound Unit cell dimensions Unit cell volume 

a (nm) b (nm) c (nm) ((nm) 3 ) 

YF 3 0.6353 0.6850 0.4393 0.1912 

HoF 3 0.6404 0.6875 0.4379 0.1928 

TbF 3 0.6513 0.6949 0.4384 0.1984 

Table 5.1: Unit cell dimensions of YF3, H0F3 and TbF3 from powder x-ray diffraction. 

significant amounts of oxide fluorides and conclude that pure H0F3, DVF3 and TbF3 
exhibit only the Pnma phase. The other rare-earth trifluorides and YF3 exhibit differ- 
ent phases at high temperatures. 

The volume per rare-earth ion in the trifluorides is ~1.5 that in the rare-earth metals. 
In H0F3 for example, the volume per holmium ion is 0.0482 (nm) 3 (table |5.1| ) whereas in 
metallic holmium the volume per holmium ion is 0.0312 (nm) 3 p[. Thus the saturation 
magnetization of H0F3 is ~65% that of metallic holmium. 



±(r a ,r b ,r c ) ±(i + r a ,i-r 6 , \-r c ) 

±(-r , 3 + n, -r c ) ±(| - r a ,-r b , \ + r c ) 

The notation (r Q , r^, r c ) is used in preference to the more usual (x,y,z) to avoid confusion with 
the coordinates used later. 



Table 5.2: Eightfold general ion positions in YF3 and H0F3 (space group Pnma), 
expressed as fractions of the unit cell dimensions. 

The positions of the ions have been determined using x-ray diffraction by several work- 



ers. The eightfold general positions are given in table 5.2. The parameters r a , ri, and 



r c required to find the positions of the 4 rare-earth ions and 12 fluorine ions in the unit 



cell are given in table 5.3. Zalkin and Templeton |7q] determined positions only for 



YF3, noting that H0F3 is isostructural with YF3. The determinations of Bukvetskii 



and Garashina [16, 13] and of Piotrowski et al 57] were both for H0F3. The two deter- 
minations for H0F3 agree fairly well although no uncertainties are given by Piotrowski 
et al. There are some much larger differences between the H0F3 and YF3 determina- 
tions. In this work the positions determined by Bukvetskii and Garashina are used. 



The positions of the four holmium ions in the unit cell are given in table 5.4, together 
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with the numerical labels which we use throughout this work. A unit cell is shown in 
figure 5.1 . 





Site & 




Zalkin 1 


Piotrowski 


Bukvetskii 3 


Ion 


symmetry 


Component 


YF 3 


HoF 3 


H0F3 


Y/Ho 


4c (m) 


r a 


0.367 


0.367 


0.3679(1) 






n 


0.25 


0.25 


0.25 






r c 


0.058 


0.061 


0.0614(1) 


F' 


4c (m) 


r a 


0.528 


0.522 


0.5224(9) 4 






n 


0.25 


0.25 


0.25 






r c 


0.601 


0.584 


0.5860(10) 4 


F" 


8d (1) 


r a 


0.165 


0.165 


0.1662(6) 






n 


0.060 


0.066 


0.0645(6) 








0.363 


0.384 


0.3797(9) 



1 Zalkin and Templeton [78 

2 Piotrowski et al 57 1. 



3 Bukvetskii and Garashina fll6f , uncertainties in parenthesis. 

4 Components transformed to conform with usual convention, (r a 



and r r 



Table 5.3: Ion positions in YF3 and H0F3 from powder x-ray data, expressed as 
fractions of the unit cell dimensions. 



Ion Position in Position in Position in 

site a directionf b directionf c directionf 

~0 0.3679 025 0.0614 

1 0.6321 0.75 0.9386 

2 0.8679 0.25 0.4386 

3 0.1321 0.75 0.5614 

f Bukvetskii and Garashina [EgI. 



Table 5.4: Positions of the four holmium ions in the H0F3 unit cell, expressed as 
fractions of the unit cell dimensions. 
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Figure 5.1: Ion positions in the H0F3 unit cell with the axes marked in fractions of 
the lattice parameters: a = 0.6404 nm, b = 0.6875 nm and c = 0.4379 nm. The filled 
(smaller) circles represent Ho 3+ ions and the shaded circles represent F~ ions. 
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5.1.1 Holmium site symmetry 



All four holmium sites in the H0F3 unit cell have C\h symmetry. That is, the only 
symmetry element is a reflection plane which, in this case, is the a-c plane. This 
makes the b direction unique and it is convenient to identify it with the z axis of 
the coordinate system (x, y, z) conventionally used to define crystal-field operators. 
Because the crystal axes are orthorhombic, we can conveniently adopt the correspon- 
dence (x,y,z) <-> (c,a,b); from now on we shall use the (x,y,z) system exclusively 
when discussing the holmium sites. 

Reflection symmetry in the x-y plane means that the lattice is invariant under the 
transformation z —> — z relative to a holmium site. Thus, in any sum over the lattice, 
terms in xz and yz must vanish. In the case of a dipole sum this gives: 

B xz = B zx = and B yz = B zy = 0, (5.1) 

where B xz means the component of the field in the x direction due to moments in the 
z direction, etc.. Similar considerations apply to the electric field gradient tensor: 



V xz = V zx = and V yz = V zy = 0, (5.2) 

where V xz means q z q x , etc.. The crystal-field parameters B\ are related to the compo- 
nents of the electric field gradient tensor by 

B° 2 oc V zz 



Bf 1 oc (V xz ± iV yz ) 



Bt 2 



oc (V 3 



ax i ^yy) 



(5.3) 



Thus the parameters B^ 1 and, in fact, all parameters B\ (k = 2, 4 and 6) with odd q 
vanish. This means that no more than 15 real parameters are required to specify 
the crystal-field interaction, whereas 27 would be required in the absence of reflection 
symmetry. 
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5.1.2 Holmium site inequivalence 



Consideration of the holmium and fluorine ion positions in H0F3 shows that each of the 
four holmium ions in the unit cell has different surroundings. However, the surroundings 
of the different sites are related by simple transformations. If we take holmium site as 
a reference, the surroundings of site 1 differ by reflection in both the x-z and y-z planes. 
Similarly, the surroundings of site 2 differ by reflection in y-z and the surroundings of 
site 3 differ by reflection in x-z. 

Consider the interactions of the holmium ions with their surroundings, crystal-field and 
Zeeman (including dipolar and exchange). These interactions are unaffected by inver- 
sion of the surroundings with respect to the ionic coordinate system; or, equivalently, 
inversion of the ionic coordinate system with respect to the surroundings. The inver- 
sion symmetry of the electronic crystal-field interaction can be seen by noting that all 
terms are of even order in the components of J. We have already established that the 
surroundings of all four holmium sites are invariant under reflection in the x-y plane. 
Combining this invariance with the coordinate inversion we see that the interactions 
will be unaffected by combined reflection in both the x-z and y-z planes (since an addi- 
tional reflection in x-y would complete an inversion of the surroundings). Hence, sites 
and 1 are magnetically equivalent. Similarly, sites 2 and 3 are magnetically equivalent. 

Now consider a magnetic field applied within the x-z plane. The inequivalent site 
surroundings are related by reflection in the x-z plane. The strengths of interactions 
on the ion are unaffected by the chirality of the surroundings, so although the response 
(the expectation value of J y for example) may differ in direction, the energy levels of 
the ion will be the same for both sites. Thus, for a field in the x-z plane, the sites will 
be indistinguishable by any technique depends only on the ionic energy eigenvalues, as 
opposed to the eigenstates. 

Similarly, we may consider a magnetic field applied within the y-z plane and note that 
the inequivalent ion sites are alternatively related by reflection in the y-z plane. Hence 
the sites will appear to be equivalent, in terms of the energy levels, if a field is applied 
in the y-z plane. There is however, no symmetry argument to support site equivalence 
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for a field applied in the x-y plane or, indeed, for any field with non-vanishing x and y 
components. Furthermore, consideration of the dipolar or crystal fields leads us to 
conclude that the sites will be inequivalent for such a field. 



5.1.3 Dipole sums 

Bleaney et al [11] give the results of a dipole sum performed over a sphere of 3.0 nm 



radius about holmium site 0. Their results are obviously subject to some numerical 
errors as the sum of the diagonal terms B xx + B yy + B zz is not quite zero as required 
by Laplace's equation (the discrepancy is ~0.6% of the largest term). This cannot be 
explained by the imperfect convergence of the sum at 3.0 nm radius, because Laplace's 
equation must hold for each term in the sum (see section 2.5.1| ). 



Table 5J3 shows the dipole sum results of Bleaney et al [|ll]] and of this work, both 



calculated using the ion positions given by Zalkin and Templeton [78]. The results are 
given in terms of the field at holmium site due to a moment of 1 [Ib on the holmium 
ions at each of the other ion sites. If the sum is performed about the other holmium 
sites, then identical results are obtained for the diagonal terms (B xx etc.). For the cross 
term B xy , the results for different sites are: -0.0022 T//xb at site due to moments at 



site 1 (as in table |5.5|) , -0.0022 T/fiB at site 1 due to moments at site 0, 0.0022 T/hb 



at site 2 due to moments at site 3, 0.0022 T//xg at site 3 due to moments at site 2 and 
zero for all other combinations. The terms B xz and B yz must vanish because of the 
site symmetry and the calculation does in fact give zero to the numerical accuracy of 
the computer program (rounding errors < 10~ 15 T//xb). 

We have also calculated the dipole sum using the ion positions of Bukvetskii and 
Garashina (T^[. These results are compared with those obtained using the parameters 
of Zalkin and Templeton |78| in table |5.6|. Although the different ion positions give 



slightly different results the two dipole sums are in good agreement. Figures 5.2 and 5.2 
show how the dipole sum converges with increasing radius. The figures show the totals 
of the components at site due to moments at all sites. Totals of the components at 



site due to equal moments at all sites are given in tables |5.5| and 5.£. These must 



be treated with extreme caution because they are physically meaningless unless the 
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moments on the holmium ions at all sites are the same. 



Ion B xx (T/(i B ) Ryy ( t /Mb) b ^z (T/a*b) ^>xy (T//xb) 





t 


i 


t 


t 


t 


t 


t 


t 





0.0327 


0.0330 


-0.0135 


-0.0133 


-0.0192 


-0.0195 


0.0000 


< 0.002 


1 


-0.0260 


-0.0260 


-0.0149 


-0.0148 


0.0408 


0.0410 


-0.0022 


< 0.002 


2 


-0.0052 


-0.0050 


0.0642 


0.0642 


-0.0590 


-0.0590 


0.0000 


< 0.002 


3 


-0.0115 


-0.0118 


-0.0214 


-0.0212 


0.0329 


0.0330 


0.0000 


< 0.002 


total 


-0.0099 


-0.0098 


0.0144 


0.0149 


-0.0045 


-0.0045 


-0.0022 





f this work, summation over a sphere of 100 nm radius. The terms fluctuate by < 0.00001 T/ [Ib 
for radii between 50 nm and 100 nm. 

% Bleaney et al pJjj, summation over a sphere of 3.0 nm radius. 



Table 5.5: Comparison of dipole sum results for holmium site in H0F3, calculated 
using the ion positions of Zalkin and Templeton. 



Ion B xx (T/fi B ) B yy (T/fi B ) B zz (T/ M b) (T/h b ) 
t t f t t t t t 

0.0327 0.0327 -0.0135 -0.0135 -0.0192 -0.0192 0.0000 0.0000 

1 -0.0260 -0.0258 -0.0149 -0.0152 0.0410 0.0410 -0.0022 -0.0023 

2 -0.0052 -0.0057 0.0642 0.0649 -0.0590 -0.0592 0.0000 0.0000 

3 -0.0115 -0.0115 -0.0214 -0.0211 0.0329 0.0327 0.0000 0.0000 
total -0.0099 -0.0103 0.0144 0.0150 -0.0045 -0.0047 -0.0022 -0.0023 

f Ion positions of Zalkin and Templeton ]7q |, summation over a sphere of radius 100 nm. The 
terms fluctuate by < 0.00001 T/^b for radii between 50 nm and 100 nm. 
I Ion positions of Bukvetskii and Garashina fig] , summation over a sphere of radius 100 nm. 
The terms fluctuate by < 0.00001 T//xg for radii between 50 nm and 100 nm. 



Table 5.6: Comparison of dipole sum results for holmium site in H0F3, calculated 
using the ion positions of Zalkin and Templeton and of Bukvetskii and Garashina. 
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Figure 5.2: Graph showing diagonal components of the dipole sum as a function of 
the radius of the sphere over which the summation was performed. 
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Figure 5.3: Enlarged graph to show convergence of the B aa component of the dipole 
sum with increasing radius of the sphere over which the summation was performed. 
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5.2 Determinations of the electronic energy levels 



There are no complete determinations of the electronic levels of the ground 5 /g manifold 
of Ho 3+ in H0F3. Sharma et al [71| studied H0F3 by optical spectroscopy and measured 
transitions from the ground multiplet 5 /§ to the excited multiplets 5 Cr4 and 3 Kg. From 
the transitions Sharma et al determined the separation of the positions of the first 
two excited levels relative to the ground level. Using these levels and the 26 levels in 
the 5 C?4 and 3 Kq excited multiplets they obtained a set of crystal field parameters by 
least-squares fitting (see section |5.3[ ). The measured levels and those calculated from 
the crystal field parameters are given in table |5.7| . The calculated values are given in 
both cm -1 and K units for comparison with optical and heat capacity results. 



Sharma et al 1 Kraus et al 2 Bleaney et a/ 3 Crystal field calculation 



(cm 1 ) 


(cm 1 ) 


(cm 1 ) 


(cm 1 ) 


(K) 


0.00 


0.00 


0.00 


0.00 


0.00 


5.90 




6.59 


6.86 


9.87 


37.99 


38 




42.96 


61.81 




81.5 




93.45 


134.46 




105 




96.65 


139.06 




165 




174.72 


251.39 








180.53 


259.75 




170 




193.43 


278.32 








217.95 


313.59 








242.21 


348.49 




253 




242.76 


349.28 








264.23 


380.94 








275.67 


396.63 








281.44 


404.94 




358 




350.89 


504.86 








370.58 


533.19 




401 




377.27 


542.82 



1 Sharma et al [[7_jJ from optical spectroscopy at temperatures from 1.6 K to 40 K. 

2 Kraus et al |47j] from Raman and infrared transmission spectroscopy, T«2 K. 

3 Bleaney et al [[ll[ from optical spectroscopy at temperatures from 1.5 K to 4.2 K. 



Table 5.7: Electronic energy levels of the 5 /g ground manifold of Ho 3+ in H0F3 

Sharma et al j7l]] quote an uncertainty of < 0.1 cm -1 on their level measurements. 
The measurements were made over the range 1.6 K to 40 K; those at temperatures 
above 1.6 K were used to identify the first two excited levels of the 5 Is manifold by 
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considering population effects. 

Kraus et al report several levels in the ground manifold, determined by infrared 
and Raman spectroscopy at 2 K. They suggest that agreement with the levels calcu- 
lated using the crystal field parameters of Sharma et al [71] is acceptable for all but 



the level at 105 cm 1 . Both sets are given in table 5/7. Bleaney et al [|ll|] report a 



splitting of 6.59 cm -1 between the ground and first excited states, measured by optical 
spectroscopy between 1.4 K and 4.2 K. The same work reports heat capacity measure- 
ments that show a Schottky-type anomaly with a maximum at ~3.5 K. This fits fairly 
well with a separation of 5.90 cm -1 which should give a peak at 3.4 K. A separation of 
6.59 cm -1 should give a peak at 3.8 K, not in good agreement with the observed value. 
Bleaney et al note that the analysis is complicated by a steeply rising lattice contribu- 
tion. However, that would only shift the observed peak to a higher temperature and so 
cannot explain the lower temperature of the observed peak. At the temperature of the 
observed peak the contribution to the specific heat resulting from the second excited 
state at 38 cm -1 (~55 K) is less than 2 x 10~ 5 of the contribution from the first excited 
state, so that cannot explain the discrepancy either. 



5.3 Crystal field 

Sharma et al obtained crystal-field parameters by fitting to the electronic energy 
levels measured in zero field (section |5.2|). Corrections to the parameters are given 



by Ram and Sharma [65] and the corrected parameters are reproduced in table 5.8. 
Sharma et al and Ram and Sharma use the WYBOURNE convention. The parameters 
have been converted to MORRISON form which is more convenient when considering 
coordinate rotations. Details of these conversions are given in Appendix A. 



Sharma et al |71[ give a second set of parameters related by 30° rotation about the 
z axis. Subject to the corrections of Ram and Sharma |65| , the two sets have been 
verified to transform into each other by ±30° rotations, to the precision given. Both 
sets have an rms deviation of 2.5 cm" 1 from the experimental data. 

At a site of C\h symmetry, reflection in the y-z plane is equivalent to a rotation of 180° 
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Sharma et al\ Sharma et al\ 



(WYBDURNE form) (MORRISON form) 





(cm 1 ) 


(K) 




(K) 





78.5 


112.95 




112.95 


B\ 


-199.6 


-287.18 


Re(-Bf) 


-287.18 


—2 


272.7 


392.36 


Im(B%) 


-392.36 


u 


17.9 


25.75 


Bo 


25.75 


B$ 

A 


-101.7 


-146.33 


Re(B%) 


-146.33 


B\ 

— 2 


-33.4 


-48.06 


Im(Bf) 


48.06 


Bi 

4 


265.8 


382.43 


Re(Bj) 


382.43 


—4 


-95.9 


-137.98 


Im(5f) 


137.98 


U 


224.6 


323.15 


S 6 


323.15 


S| 


156.6 


225.32 


Re(5f) 


225.32 




-552.3 


-794.65 


Imjsf) 


794.65 


R6 
"4 


58.1 


83.59 


Re(-Bf) 


83.59 


R6 

-D_4 


-48.12 


-69.24 


Im(Bf) 


69.24 


R 6 


-405.9 


-584.01 


Mb!) 


-584.01 


R 6 

•°-6 


31.3 


45.03 


Im(Bl) 


-45.03 



t Sharma et al [ fry from fitting to electronic energy levels determined by optical spectroscopy. 
X As f, but convention changed to MORRISON so that coordinate rotations can be calculated 
easily. Crystal-field parameter conventions are discussed in appendix A, coordinate rotations 
are discussed in appendix B. 



Table 5.8: Crystal field parameters for H0F3. 

about the y axis. The rotation can be represented by reflections in the y-z and x-y 
planes, but the lattice is invariant under reflection in x-y plane so the rotation is simply 
equivalent to reflection in the y-z plane. In the notation of appendix B, a rotation of 



180° about y is represented by D(0, 7r,0) = d(7r). Equation B.13 gives 



= (-l) W 4tJ(0), (5-4) 



and noting that d^(0) = 5 q > q , this simplifies to 

dg*,(ir) = {-l) k - q '5 q , q . (5.5) 

With crystal field parameters in MORRISON form (see appendix A), the transformed 
parameters M' k are given by equation B.5: 



k 

q'=-k 
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This simplifies to 



M ,k = (_!)fe-9 M fe (5.7) 



and we note that M k q = (-l) q (M k )* (see section A. 4). Also, (-l) fc = 1 as k is even 
and the final transformation is 



Mf = (Af*)*. (5.8) 
Rotation of 180° about the x axis is represented by D(— §,7r, |). Combining equa- 



tions |5.8| and [B.9| gives an expression for the transformed parameters M" k : 



Mf = e i «i(M a fc )*e- i "i (5.9) 



= «)*. (5.10) 

Rotation of 180° about both the x and y axes is equivalent to the reflections in the y-z 
and x-z planes. In C\h symmetry this transformation leaves the parameters unchanged: 

M' q " k = {{M k )*y = M k q . (5.11) 

We have now established the transformations required to find the crystal field param- 
eters for all four holmium sites, given a set a parameters for one. Most importantly we 
see that although there are four sites, there are only two sets of crystal field param- 
eters. Specifically, sites and 1 have the same crystal field parameters (say M k for 
site and thus M'' lk = M k for site 1); sites 2 and 3 share a related set of parameters 
{M' k = M'^ k = (M k )*). This accords with the argument given in section 15.1.21. 
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5.4 Magnetometry 



We have made magnetometry measurements on H0F3 at 5 K with magnetic fields 
applied parallel to the principal crystallographic directions. These measurements were 
made using the Manchester vibrating sample magnetometer (VSM) with applied fields 
in the range to 10 T. The temperature of 5 K was chosen as the lowest easily attainable 
with the continuous-flow cooling system. We estimate the sample temperature to have 
been within ±1 K of 5 K, and the crystal alignment to have been within 5° of the 
nominal directions. 

The dimensions of the crystal were 3.2x1.98x1.66 mm parallel to the a, b and c direc- 
tions respectively, and the mass was 74.8±0.1 mg. The magnetization data have been 



scaled to give the moment per holmium ion and are shown in figure 5.4 



From the data in figure 5.4, the susceptibilities in vanishing applied field are 6.0, 



0.48 and 0.88 \xb atom 1 T 1 for fields along the a, b and c axes respectively. By 



extrapolating the data of Bleaney et al [11] to 5 K, the a axis susceptibility is ~ 
6.7 atom -1 T _1 . Agreement is satisfactory given that the uncertainty in the tem- 
perature of the sample is ±1 K. 



Figure 5.5 shows expectation value of the ionic moment calculated from the magnetiza- 



tion data and the results of crystal-field calculations for each of the principal crystallo- 



graphic directions. The calculations and parameters used are discussed in section 5.4.1. 
Here we use the notation (J r ) s for the expectation value of J r when the applied field is 
along the s direction, and («7} s for the expectation value of J when the field is applied 
along the s direction. Thus, | («7) s | is the magnitude of the ionic moment for a field 
applied along the s direction and {J s ) s is the colinear component of the moment for a 
field applied along the s direction. 



The particularly striking feature of figure 5J3 is the difference between ( J c ) c and | {J 



For 8 T applied along the c-axis, the calculation gives moments canted at 64° to the 
field. The moments at sites and 1 are canted in the opposite direction those at sites 
2 and 3. The non-colinear part of the moment cancels when averaged over all four 
holmium sites so the bulk magnetization is determined by ( J c ) c . Although not in good 
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Applied field, B (T) 



9 10 



Figure 5.4: Magnetization curves for a H0F3 crystal with the field applied parallel to 
the a, b and c axes. Magnetization is given in terms of the moment per holmium ion 
parallel to the applied field. The symbols represent magnetometry measurements and 
the lines are calculated from the crystal-field parameters. The calculated expectation 
values have been scaled by gj = 1.2417. The experimental data were taken with the 
sample at 5 K. 
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Figure 5.5: Magnetic moment per holmium ion calculated for fields applied along the 
principal crystallographic directions. For each direction the magnitude of the moment 
and the magnitude of the moment parallel to the field are plotted. In the case of the 
field parallel to the 6-axis, these are identical as is expected from the site symmetry. 
The data shown are thermal averages calculated for T = 5 K. 
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agreement with the magnetometry data, ( J c ) c is much closer to the measured value than 
| (J) c |- The NMR frequencies, on the other hand, are determined by | {J) c I- So, an 
NMR experiment with the field applied along the c axis would provide an independent 
test of the crystal-field parameters. 



Figure |5.5| shows the thermally averaged moments calculated for 5 K as would be 
measured by magnetometry. However, in rare-earth insulators at liquid helium tem- 
peratures we expect the NMR frequencies to be determined by the electronic ground 



state moment. For example, Bunbury et al IS] have shown this to be the case in 
Ho(OH)3. In practice the values are very similar at moderate to high fields (< 0.4 % 
difference between 5 K thermal average and ground state expectation values at 3 T) 
because only the ground state is appreciably populated. 

5.4.1 Crystal-field calculations 

Crystal-field calculations were performed using the parameter conversion and calcu- 
lation programs described in Appendix C. We have used the crystal-field parameters 



of Sharma et al JFl]] (table |5^), and the ionic parameters given in table 2.1. The 
data shown were calculated using a purely electronic Hamiltonian which makes the 
calculations very much faster because the Hamiltonian is only 17 by 17 rather than 
136 by 136 if the nucleus is included. To check that the effect of the hyperfine inter- 
action was not significant a few additional calculations using the combined electronic 
and nuclear Hamiltonian were also performed. They gave very similar results so we 
have confidence in the purely electronic calculations. We have neglected the dipolar 
and exchange interactions which give corrections of < 0.2 T at the ion with 8 T ap- 
plied along the a axis. Such corrections do not significantly affect the agreement of the 
experimental and calculated data. 
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5.5 NMR 



Bleaney et al report field-swept NMR measurements on a 'self-aligned' sphere of 
H0F3. By applying a 1.5 T field prior to the experiment the sphere was aligned with 
the crystallographic a-axis parallel to the field. At 3.52 GHz they found a spectrum 
with marginally resolved quadrupole structure (at T = 1.61 K). The centre of their 
spectrum corresponds to an applied field of 0.559 T. To obtain the field experienced 
by the Ho ions, dipolar field and exchange field corrections were included. Prom their 
susceptibility measurements the moment per Ho ion is 5.01//^ so, from the dipole sum, 
Bdip = 0.072 T. The exchange field was estimated from susceptibility measurements, 
giving B ex = 0.094 T. This gives a total field at the ion of B = 0.725 T. 



From our crystal- field calculations (see section 5. 5.3] ), we obtain an NMR spectrum 



centred at 3.145 GHz in a field of 0.725 T. This is about 12% lower than the frequency 



given by Bleaney et al [11]. To give the central line at 3.52 GHz from our crystal-field 
calculation would require an ionic field of 0.875 T, which suggests the need for better 
crystal-field parameters. 



Bleaney et al [11] report a quadrupole splitting of 146(4) MHz which is in fair agree- 
ment with our calculated value of 132 MHz at 0.725 T. The quadrupole structure of 
their spectrum is only just resolved, with a line width of ~ 0.04 T. This line width in 
field corresponds to a line width in frequency of ~ 130 MHz (from the calculated ionic 
moment as a function of field). If the line width is assumed to be independent of fre- 
quency, then their observed line width disagrees with that of our NMR data, which was 
taken at higher fields and gives a significantly smaller line width in frequency. However, 
our data show very similar line widths for concentrated and dilute samples, suggesting 
that the dominant inhomogeneous broadening mechanism is crystal imperfection. If 
that is the case then we would expect the line width to be increased by the greater 
variation in quenching at low fields. 

In this work we have used CW and pulsed NMR to study H0F3 and l%Ho:YF3. The 
crystals have been both 'self-aligned' (H0F3 only) and oriented with the applied field 
nominally along the a direction. Most of the data have been taken by CW NMR with 
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the sample and helium bath at 1.5 K to avoid frequency instability caused by helium 
bubbling. The pulsed NMR data have been used to confirm that certain features of 
the CW spectra are not instrumental artifacts. 



5.5.1 Self-aligned sample 

We have taken NMR spectra from a small, nearly spherical, H0F3 sample that was 
allowed to move freely in the sample holder. By applying a field of 8 T the sample was 
aligned along the crystallographic a axis (the easy direction). 

Although the sample was only about 1 mm in diameter, reasonable spectra were ob- 
tained; a representative spectrum is shown in figure |BT6| (upper graph). This spectrum 
illustrates some of the problems of frequency modulation and third harmonic detec- 
tion: there is a large background and the resonances are hard to interpret. The largest 
uncertainty in the resonance positions is due to resonance shape, and is estimated to 



be ±0.05 T. The lower graph of figure 5.6 shows that there is also a significant shift in 
the cavity resonance during the sweep, and some frequency instability toward the high 
field end of the sweep. 



Figure 5.7 summarises the spectroscopic data. Between 4.95 GHz and 5.15 GHz it was 
not possible to resolve individual resonance lines so only the estimated centres of the 
spectra are plotted. The 'collapse' of the quadrupole splitting is attributed to the fact 
that the intra- and extra-ionic contributions to the quadrupole splitting, which have 
opposite signs, cancel each other at a field of the order of 2 T. The quadrupole splitting 
is predominantly extra-ionic at low fields and predominantly intro-ionic at high- fields. 
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Figure 5.6: A 7 line holmium NMR spectrum from a small self-aligned H0F3 sam- 
ple. The vertical lines on the upper graph show the inferred line centres plotted to 
the nearest 0.05 T. The data were taken using frequency modulated CW NMR with 
third harmonic detection at 50 kHz, 1 mW incident microwave power, and ±1.7 MHz 
frequency modulation (considerably smaller than the NMR line width). The lower 
graph shows variation of the resonator frequency during the field sweep. The sample 
temperature was 1.5 ± 0.1 K. 
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Figure 5.7: Graph showing NMR resonance frequencies against applied field for a 
self-aligned H0F3 sample. Several sweeps were taken at some frequencies and in these 
cases the average resonance positions have been plotted. The data were taken using 
frequency modulated CW NMR and third harmonic detection. The sample temperature 
was 1.5 ±0.1 K. 
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5.5.2 Two spectra 



All the NMR data except for those from the self-aligned sample show two distinct but 



overlapping spectra. Figure 5.5 shows two particularly well resolved spectra taken by 
frequency modulated CW NMR at high field. Figure |5.9| shows two marginally resolved 
spectra taken from the same sample using spin-echo NMR. That two spectra are seen by 
both the CW and pulsed NMR provides strong evidence that they are not an artefact 
of either technique. Although we have obtained resolved spectra from both CW and 
pulsed NMR, the resolution of our CW NMR data is superior. Since these spectra 
are obtained from the magnetically dilute sample (l%Ho:YF3) we cannot attribute the 
presence of two spectra to different dipolar fields at different sites. 



Figure 5,10| shows three sets of resonance positions for l%Ho:YF3. Sets a) and b) were 



taken using frequency modulated CW NMR and set c) was taken using spin-echo NMR. 
The difference between sets a) and b) corresponds to a slight change in the orientation 
of the sample such that in set b) the two spectra are less well separated. The agreement 
of sets b) and c) (between which there was no mechanical movement) indicates that 
the effect is not an artefact of the CW technique. The uncertainties shown arise from 
the difficulty in locating the centres of the resonance lines. 

The situation is confusing in some cases where there spectra appear to show 8 lines. 
This is attributed to there being two 7 line spectra shifted by an amount equal to the 



quadrupole splitting. Figure 5.11 shows one such spectrum, in this case taken by field 



modulated CW NMR for which the line shape is a simple derivative. 



Figure 5.12 shows the calculated values of | (J) | for a misaligned crystal. We consider 
the ground state only. The coordinate system is chosen such that 9 = 0° corresponds to 
the a axis and 9 = 90°, cp = 0° corresponds to the b axis. For a field applied along the 
a axis (9 = 0°), | {J) |= 6.6338 and the sites are equivalent. This confirms invariance 
under he C\h symmetry operation, ie. reflection in the a-c plane, which corresponds 
to <p — > 180 — <p. Also, reflection in the a-b plane exchanges sites and 1 with sites 
2 and 3, and corresponds to (p — > —(p. Taking 9 = 5°, the difference between the two 
sites is maximum for (p = 90° and gives /S. s i tes | (J) | = 0.0085 which is equivalent to 
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Figure 5.8: Two clearly resolved 7 line holmium spectra from l%Ho:YF3 (upper 
graph). The sample was aligned with the crystallographic a axis nominally parallel to 
the applied field. The data were taken using frequency modulated CW NMR with third 
harmonic detection at 50 kHz, 1 mW incident microwave power, ±1.7 MHz frequency 
modulation, and loaded cavity Q ps 500. The lower graph shows variation of the 
resonator frequency during the field sweep. The sample temperature was 1.5 ± 0.1 K. 
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Figure 5.9: Two resolved 7 line holmium spectra from l%Ho:YF3 taken using spin- 
echo NMR at 5.300 GHz. The sample was aligned with the crystallographic a axis 
nominally parallel to the applied field. The sample temperature was 1.5 K. 
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Figure 5.10: Graph showing the positions of the resonances in field for three successive 
NMR experiments performed at the same frequency and on the same sample. Between 
taking set a) and set b) the sample was removed and replaced. The frequency shifts 
are attributed to a small change in the orientation of the sample. The two experiments 
were otherwise performed under identical conditions. Set c) was taken by spin-echo 
NMR without moving the sample from set b). In all cases the frequency was 5.3000 db 
0.0005 GHz and the sample temperature was 1.5 ± 0.1 K. 

A s it es B ap p ~ 0.26 T. This is larger than the observed field shift between any of the 
spectra. We may conclude that the observation of two distinct spectra is adequately 
explained by a misalignment of < 5°. 



5.5.3 Crystal-field calculations 

Crystal-field calculations were performed use the parameter conversion and calculation 
programs described in Appendix C. As before, we have used the crystal-field parameters 



of Sharma et al [71| (table |5.8|) , and the electronic and nuclear parameters given in 



table 2.1. In particular, the extra-ionic quadrupole parameter Q ex t = 9.66 x 10 4 was 



derived from the ratio of the nuclear to electronic antishielding factors of Carboni |2C] 
(table [5.8D for holmium ethylsulphate. (The nuclear and electronic antishielding factors 
vary considerably between different compounds and are not know for H0F3.) We have 
used the ratio for holmium ethylsulphate because it gives reasonable agreement with the 
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Figure 5.11: Graph showing an '8 line' holmium spectrum. The two 7 line holmium 
spectra are shifted by an amount equal to the quadrupole splitting so there appears to 
be just a single 8 line spectrum. These data were taken by field modulated CW NMR. 
The modulation frequency was 107 Hz with ±6 mT deviation and 1 mW incident 
microwave power. A small frequency modulation at 100 kHz was used to track the 
resonator resonance; the sample temperature was 1.5 K. 
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Figure 5.12: Graph showing the calculated variation of | (J) | as a function of mis- 
alignment relative to the crystallographic a direction. Curves are given for 6 = 1° and 
6 = 5°. In each case there are two curves: one for holmium sites and 1, the other for 
sites 2 and 3. In the case 9 = 0°, tp undefined and | (J) | is slightly higher than the 
mean value for 6 = 1°. 
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observed quadrupole splittings, much better than that for holmium hydroxide (also in 
table [O]) for example. The collapse of the quadrupole splitting is particularly sensitive 
to the antishielding ratio. If there were better agreement between the experimental 
and calculated data, it would be possible to determine the ratio of antishielding factors 
for H0F3 by fitting. However, given the discrepancy between the data and and the 
crystal-field calculation, we take the holmium ethylsulphate value as a reasonable first 
approximation to that for H0F3. 



Figure summarizes our CW NMR data for l%Ho:YF3 and shows the NMR fre- 



quencies predicted from our crystal-field calculations. The dimensions of the crystal 
were 8.0 mm parallel to the a direction and 2.0 x 1.8 mm parallel to the b and c directions 
(order unknown). The sample was mounted with the a direction nominally parallel to 
the applied field. Because the sample is magnetically dilute (1% holmium) we do not 
expect any significant corrections for dipolar field, exchange field or sample shape. It 
is thus reasonable to make direct comparison with calculated NMR frequencies shown 
in figure [5.13 , Although the qualitative agreement is good, there are significant dis- 



crepancies. The experimental and calculated dipolar splittings differ by 2-3%, and the 
calculated quadrupole splittings are ~ 20% greater than the experimental values. How- 
ever, we have already noted that the calculated quadrupole splittings depend strongly 
on the antishielding ratio, which is not known with any accuracy. 

An interesting aspect of the collapse in the quadrupole splitting is that we observed 
spectra consisting of just two, partially resolved, lines between 5.0 and 5.1 GHz. We 
attribute the two lines to the two inequivalent sites, where the quadrupole structure 
within each is not resolved. 
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Figure 5.13: Graph showing NMR resonance frequencies against applied field for 
l%Ho:YF3. The points (joined with dotted lines as a guide to the eye) are experimental 
points, and the solid lines are from crystal-field calculation. The sample was oriented 
with the crystallographic a axis nominally parallel to the applied field. The data were 
taken using frequency modulated CW NMR and third harmonic detection, The sample 
temperature was 1.5 ± 0.1 K. 
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5.5.4 Dipolar and exchange fields 



There is a significant difference between the spectra from the concentrated (H0F3) and 
dilute (l%Ho:YF3) samples. We attribute this to the dipolar and exchange fields in 
the concentrated sample which will be negligible in the dilute sample. The difference 



is illustrated by figure |5.14| where the spectroscopic data from both samples is plotted 
together. 



In addition to the CW NMR data plotted in figure 5.14 , C. Carboni of this depart- 
ment has made spin-echo measurements on the same H0F3 sample. These are in good 
agreement with the results of this work. However, his data show only a single holmium 
spectrum which we attribute to better alignment of the a axis with the applied field. 
The quadrupole splitting vanishes at 1.65 T and 5.05 GHz in the spin-echo data. 

At 5.2GHz, the centre of the spectrum from the concentrated sample is about 0.45 T 
lower in field than the spectrum from the dilute sample. Assuming that both the 
exchange and dipolar fields are negligible in the dilute sample, we conclude that B^ip + 
Bex = 0.45 T in H0F3 (under the stated conditions). In order to calculate B^p, we 
take the ionic moment as 8.0 fiB, corresponding to the NMR frequency of 5.2 GHz. 
Ignoring a small cosine error (< 5%) associated with the canted spin structure gives 
a bulk magnetization M ~ 1.54 x 10 6 Am-1. The dimensions of the crystal were 
8.0 x 1.98 x 1.66 mm parallel to the a, b and c directions respectively. Treating the 
crystal and an ellipsoid, with principal axes proportional to the dimensions, we obtain 
an estimated demagnetizing factor, N ~ 0.065, by interpolation from the tables of 
Cronemeyer p7||. We conclude that the Lorentz and demagnetization contributions to 
the dipolar field are Bl + B^ m = (~ — 0.065)/ioM 0.52 T. Using the dipole sums 
given in section 2J5, with B yy = B aa = 0.015 T/x^ 1 (see table 5J3), and ignoring effects 



due to canting, we estimate Bi nt ~ 0.12 T (the internal field or contribution to the 
dipolar field due to ions within the Lorentz sphere). Thus B^ip ~ 0.64 T and hence 
B ex = -0.19 T. The fact that B ex is opposite and approximately equal to B^ n % accords 
well with our observation that the NMR frequencies from the almost spherical (self- 
aligned) H0F3 sample closely match those of the dilute material. However, our value 



for the exchange field does not agree with the results of Bleaney et al [11|. Scaling 
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Figure 5.14: Graph showing NMR resonance spectra for both H0F3 ('concentrated', 
solid lines) and l%Ho:YF3 ('dilute', dotted lines). All the data were take using fre- 
quency modulated CW NMR at 1.5 ± 0.1 K. 
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their value of B ex = 0.094 T for an ionic moment of 5.01 \xb gives B ex ~ 0.15 T for an 
ionic moment of 8.0 [is- The origin of this discrepancy is unclear. It is unlikely that it 
can be explained by Bleaney's neglect (and ours) of the effect of canting on the dipolar 
field. 



5.5.5 Other predictions from the crystal-field parameters 

We have made a set of crystal-field calculations for fields applied in the a-c plane. The 
crystal-field axes were chosen such that the 6-axis corresponds to 6 = 0° and hence 
for 9 = 90°, ip covers the a-c plane. The ground-state values of | (J) | are shown in 



figure |5.15| as a function of the field direction ip, where ip = 0° corresponds to the field 
along the c axis. The two curves represent the two different pairs of sites (the solid 
line for sites and 1). This confirms that the calculation does not distinguish between 
the sites for fields along the c and a axes ((p = 0° and 90°). We also see that the 
responses from the different sites are related by reflection in the a-b plane {<p — ► —(p) 



or the c-b plane (ip — > 180 — (p). However, the most striking features of figure 5.15| 



are the sharp dips in | (J) | at <p ~ 20° and 160°. To find out whether this result is 
sensitive to the crystal-field parameters, several more calculations were performed with 
modified parameters. Qualitatively the results appears to be fairly insensitive to the 
parameters. For example, a similar dip is produced even if M| is set to zero. We have 



no experimental evidence for or against the behaviour predicted in figure 5.15 



It is interesting to note that our crystal-field calculations do not give the crystallo- 
graphic a axis as the easy axis for either site alone. This point is illustrated buy 
figure 5.15| . However, the calculations do suggest that the a-axis should be the easy 



axis for a crystal composed of holmium with equal numbers of holmium ions on the two 
inequivalent sites. The 'averaged' easy axis is then in agreement with our experimental 



data and that of Bleaney et al [11]. 
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Figure 5.15: Magnitude of the ground state expectation value of J, | (J) |, as a 
function of the direction of the field of 8 T applied in the a-c plane. The angle ip covers 
the a-c plane, and tp = 0° corresponds to the field applied along the c direction. The 
two lines are for the inequivalent holmium sites: one for sites and 1, the other for 
sites 2 and 3. 

5.5.6 Nuclear relaxation 

From spin-echo measurements on l%Ho:YF3 taken at 5.300 GHz and 4.482 T we have 
estimated the spin-spin relaxation time to be T2 = 7 ± 2 fis. For this we have assumed 
exponential relaxation which gives a fairly good fit to the experimental data. We have 
attempted to measure 1\, which turns out to be very long compared to that in other 
rare-earth insulators, such as the hydroxides. This caused problems because the pulse 
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spectrometer has a maximum pulse separation if 1 s. However, by observing multiple 
stimulated echoes following a single spin-echo experiment, we estimate T\ to be of the 
order of 2 s. No relaxation measurements were made on H0F3. 



5.6 Summary 

From analysis of the crystal structure of H0F3 we have deduced that the holmium 
sites are not magnetically equivalent, unless the applied field lies in either the a-b 
or b-c plane. The a direction is the easy direction, and NMR results from a self- 
aligned H0F3 show a single, well resolved, 7-line holmium spectrum. In support of our 
theoretical prediction that there are two magnetically inequivalent sites, we have NMR 
results which show two clearly resolved holmium spectra. However, these results are 
from samples aligned with the a direction nominally parallel to the applied field. To 
reconcile the experimental results with our theoretical prediction, we have considered 



wc 



a small misalignment of the sample from the nominal a direction. In section p.5.2 
showed that a misalignment of < 5° can give the observed separation between the two 
holmium spectra. Further supporting the idea of random misalignment we have shown 
different spectra from a single sample and taken in the same conditions except that the 
sample has been removed and replaced between experiments. 



The crystal-field parameters of Sharma et al [71] provide good qualitative agreement 
with our NMR and magnetometry results. The quantitative agreement is fair, with 
a discrepancy of 2-3% between the calculated and measured dipolar splittings in 



yo 



Ho:YF3. This suggests the need for more accurate crystal-field parameters. 



5.7 Further work 



The most useful additional information would be a more accurate set of crystal-field 
parameters. These would be best derived from a set of electronic levels that included 
the complete ground manifold. The parameters of Sharma et al (71]] are derived from 
the lowest 3 levels of the 17- level ground manifold and levels from excited multiplets. 
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It is reasonable to expect that parameters derived from levels including more from the 
ground manifold would give a better descriptions of NMR and magnetometry measure- 
ments. 

More spectroscopic information for the concentrated and dilute samples would allow a 
better analysis of the dipolar and exchange fields. The use of spherical samples would 
remove the need for demagnetizing field corrections which we have shown to be sub- 
stantial. To test the predictions of out crystal-field calculations, an NMR experiment 
with the applied field along the c axis would be interesting. The NMR frequencies 
would immediately support of disprove the prediction that the moments are strongly 
canted away from the c axis. However, H0F3 is clearly an extremely anisotropic system, 
so it is reasonable to expect that the transverse enhancement factor will be anisotropic. 
It may be that the enhancement factor is high for fields along the a axis, which would 
help explain the strong signals seen in this work. But, it may be low in other directions 
making NMR more technically demanding. It would be desirable to use the crystal-field 
model to estimate enhancement factors for other directions before attempting NMR. 

As part of this work, we have made several attempts to perform an NMR experiment 
with a H0F3 sample mounted such that it could be rotated, putting the applied field 
anywhere in the a-c plane. All attempts have failed either because of mechanical fouling 
of the rotation device, or because the crystal broke or shattered in an applied field. In 
spite of the technical difficulties, an investigation of either magnetization or NMR as a 
function of orientation of the field in the a-c plane would be very useful. Our crystal 
field calculations predict 'hard' directions which could be checked by either technique. 
We also note that NMR with an applied field in the a-c plane will, in general, give 
distinct spectra for the inequivalent holmium sites. 
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A. Parameter conventions for 
crystal fields 



There is a great deal of notational confusion in the literature on crystal fields. Much of 
the problem stems from the magnetism community, which has been concerned almost 
exclusively with the lowest J-manifold of rare-earth ions in situations of high sym- 
metry. The simplicity of such systems has encouraged the use of ad hoc crystal-field 
parametrization schemes. These, though adequate for their immediate purpose, cannot 
be generalised in a consistent manner to situations of greater complexity. Conventions 
which are adequate to describe the lowest J-manifold of a system described by 2, 4 or 
even 6 crystal-field parameters tend to break down for low-symmetry systems, which 
may require up to 27 crystal-field parameters, especially if excited manifolds are to be 
included in the calculation. For these, a more coherent approach is required. 

In this appendix we review the development of crystal-field notation, starting with the 



still widely used operator equivalents of Stevens 74]. Other commonly used notations 
are then described, followed by a brief note on computation using the 3-j symbols. 
Finally, we present a detailed explanation of the inter-relationships of parameters in 
the notations described. 

All notations discussed here use operators with symmetries related to those of the 
spherical harmonics. Parametrizing the interactions of 4/ electrons with the crystal 
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field using such operators is convenient because it is relatively easy to write down the 
matrix elements between the free-ion electron states used as a basis for calculations. If 
B\ represents the coefficient of an operator related to the spherical harmonic Yk q (as 
defined in equation A.22| ) then the ranges of k and q are limited as follows (see, for 



example, [f>I|| ): 

• For a k th order operator, —k < q < k. 

• As / = 3 for 4/ electrons and only terms for which k < 21 are required, k < 6. 

• Time-reversal invariance requires k to be even. 

• To represent a real potential, B^ q = (—l) q (B^,)* which implies that the must 
be real and that we only need (2k + 1) real parameters for a given k. 

• The Bq term need not be included as it only produces a shift in all the energies, 
not a change in the splittings. 

• Finally, all coefficients B^ with q ^ are subject to a phase factor which depends 
on the choice of x and y axes. If there are no crystallographic grounds to pick 
particular x and y axes, then the phase factor is usually chosen to set Im(l?2) or 
Im(S|) to zero. 

These conditions result in a maximum of 27 parameters (only 26 of which are indepen- 
dent: one may be set to zero by choice of axes) which may be expressed in terms of 3 
real and 12 complex parameters: 

rO r 1 r2 
JD 2 i JD 2 i JD 2 ) 

RO d1 r2 R 3 r4 

Bl BlBl Bl,Bt Bl and S|. 
The B® are real and the B^ (q ^ 0) are complex. Some notations use wholly real 
parameters by including Bf, for negative q (—k < q < k) with operators related to the 
tesseral harmonics. 

It is not usually the case that all of these parameters are required. The number of 
parameters is limited by the symmetry of the site of interest since the crystal-field 
Hamiltonian must be invariant under the operations of the symmetry group. However, 
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the symmetry operations are tied to a set of axes so the number of parameters will de- 
pend on the choice of axes. For example, the surroundings of a site with C\h symmetry 
are invariant under reflection in the x-y plane. If the x and y axes are not chosen to 
lie in this plane of reflection symmetry then, in general, 26 rather than 14 parameters 



will be required. The effect of the site symmetry is illustrated in table A.l, where the 
parameters required for all 32 point groups are listed. 



Site symmetryf 



Class 



Crystal-field parameters required 



C3h, 



O h , T d , O, 
T h , T 

D&h, D 3 h, D, 

Dih, D2d, C^v, D4 
D3, D^d, C 3v 

C4, C^h, 5*4 

D2, D 2 h, C2v 



C2, C 2 h 
Cm {C s , Cu 



C\, Ci 



Cubic| 

Hexagonal 

Tetragonal 
Trigonal 

Tetragonal 

Trigonal 

Orthorhombic 

Mono clinic 

Triclinic 



B°2 
B°2 

B% 

R° 
°2 

B% 

B°2 
B°2 

Bl 
R° 

B°2 

Bl 
R° 

B°2 

Bl 
Bl 



Re(Bj), B, 



R° R° 



Re(£, 

Mb!). 



Bl, Re(Bj), B° 6 , Re(Bl 

Bl Re(Bf), 
Re(Bl), Re(Bl). 

Bl Re(Bj), 
Bl 

Bl Re(Bf), 
Bl Bl 



Re(B 2 



2h 



■1J- 



Re(Bj), Re(B, 
Re(Bi), Re (,R|), Re(5 6 6 ). 

Re(flf), 
Bl, B A , 
Bl Bl Bl 

Re( J B 2 1 ), Bl 

nl n2 d3 d4 

r> 4 , r> 4 , 234, JD4, 

B l Bl, Bl Bl Bl Bl 



f Point groups are used to describe the site symmetry, for an explanation of the labels see 
Sachs J6^] for example. 

X The four parameters for cubic symmetry are not all independent (see section |A.3| ) . 

Parameters B^, are complex for q > 0, where only the real part is required it is written explicitly 
as Ke(B'ji). In this tabulation the phase factor determined by the orientation of the a; and y 
axes is chosen to give the minimum number of parameters, for example in the triclinic class B\ 



is made real. Tabulation after Walter 76 



Table A.l: Parameters required to describe the crystal-field interaction in the 32 point 
group symmetries. 
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A.l Stevens operator-equivalent notation 



In Stevens's original paper ||74|, four potential functions were considered. These are 



reproduced in table A. 2; the summations remind us that crystal-field theory was based 



on the idea of independent electrons in a non-spherically symmetric potential. 

y 2 ° = £(3z 2 -r 2 ) 

y 4 ° = ]T(35z 4 -30rV + 3r 4 ) 

y 6 ° = J](231z 6 - 315rV + 105r 4 z 2 - 5r 6 ) 

V$ = J2(x 6 ~ 15xV + 15*V - y 6 ) 

Table A. 2: Potential functions considered by Stevens. Summations are over the 
coordinates of all electrons. 

Each V™ may be factorized and written in the form f n (r)Y™(9, (p) so as to transform 
according to an irreducible representation of the rotation group. Stevens notes that 
there are simple relations between the matrix elements of potential operators and ap- 
propriate angular momentum operators inside a manifold of constant J (a result that 
follows from the Wigner-Eckart theorem, see, for example, Brink and Satchler |l2|). By 
comparing the matrix elements of the potential operators and of angular momentum 
operators, Stevens determined equivalent angular momentum operators. In generating 
such operator equivalents allowance must be made for the non-commutation of J x , J y 
and J z . Operator equivalents that have been symmetrized to allow for non-commutation 



are given and are reproduced in table A3. The factors a, (3 and 7 are referred to as 



operator- equivalent coefficients and relate the angular momentum operators to the po- 
tential operators. Stevens 0] lists the values of the operator-equivalent coefficients (a, 
(3 and 7), together with the matrix elements of V®, V® and Vq for the ground manifolds 
of all the rare-earth ions (in the L-S coupling scheme). Values of a, (3 and 7 corrected 
for intermediate coupling (a significant factor in many rare-earth ions) are also avail- 



able: see, for example, Wybourne [77]. The factors (r n ) are expectation values of r r 



for the 4/ electrons and thus depend on both the ion and the configuration. Freeman 



and Desclaux [37] tabulate values of (r n ) from relativistic calculation. 
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V 2 °^a(r 2 )[3J 2 z - J(J + 1)] 

V~ 4 ° = /3(r 4 )[354 - 30J(J + 1) + 25 J 2 - 6J(J + 1) + 3 J 2 (J + l) 2 ] 
V 6 ° = 7 (r 6 )[231J 2 6 - 315J(J + 1)4 + 735J 4 + 105J 2 (J + l) 2 J 2 - 525J(J + l)jf 
+294 J 2 - 5 J 3 (J + l) 3 + 40 J 2 (J + l) 2 - 60 J (J + 1)] 



The expressions for V 2 °, V4 and Vg° are due to Stevens Q; that for Vg 6 is due to Elliott and 
Stevens il. 



Table A. 3: Operator equivalents inside a manifold of constant J. 

In applying Stevens' method to cerium ethyl sulphate, Elliott and Stevens |3^] give a 
more general expansion of the crystal-field potential in terms of spherical harmonics: 

V cf= E E A%(r n )Y™(e,<p). (A.l) 

n=2,4,6m=0...n 

where the A™ are numerical parameters and the Y™ (9, ip) are spherical harmonics. In 
the case of / electrons in C^h symmetry this is reduces the four terms discussed by 



Stevens [74] but Elliott and Stevens note that the Vg term is a combination of Yg and 

y-6 



Elliott and Stevens |35| later considered J-mixing, the admixture of other J-manifolds 
by the crystal field. They calculated operator-equivalent coefficients between states 
with different total angular momenta, still however, assuming L-S coupling. They 
also note similarities between the operator-equivalent coefficients and Racah's coeffi- 



cients [61] and give an expression for the operator equivalent of Vg (table |A.3j ). In a 



later paper [34] they tabulate the matrix elements of Vg 6 . In order to extend crystal- 
field calculations to sites with C^ v symmetry Judd [^5|] added operator equivalents for 
Vf and Vq together with tables of their matrix elements within a constant J manifold. 

A significant shift in notation was introduced by Baker, Bleaney and Hayes Q (here- 
inafter referred to as BBH), who recast the crystal-field interaction in the form of an 
explicitly quantum-mechanical Hamiltonian, 7i c f. 
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Kef = 12 22 KO™ (A.2) 

n=2,4,6m=0...n 

within a manifold of given J, where 

B n = A™{r n )(J || a n \\ J). (A.3) 

The operator-equivalent coefficients a, (3 and 7 are rewritten in the more general form 
(J || a n || J) where n = 2,4,6 respectively. The A™ are the same as those of equa- 
tion |A.1| . The operator equivalents do not include the operator-equivalent coeffi- 
cients or the radial factors (r n ). Both factors are included in the parameters B™, which 
restricts application of the -B™ a single J-manifold. The parameters A™ are (in princi- 
ple) common to an entire L5-term. BBH tabulate the matrix elements of 0\ and 0\ for 
application to tetragonal symmetry. They also obtain a complete 'spin Hamiltonian' 
for the manifold by adding a Zeeman term to Hqf- Using the same notation, Jones 
et al [^4| tabulate the operator equivalents 0\, 0\ and 0\ and their matrix elements, 
thereby extending the notation to orthorhombic symmetry. Corrections to some of the 



matrix elements 03, 44, WM are given by Hutchings H3]. A complete tabulation of 



for n = 2,4,6 and m = 0, 1 . . .n (equation |A.2| ) is given by Orbach |55] and is 



reproduced in table A. 4. 



A similar notation is used by Scott [fF(]] and by Abragam and Bleaney |T|] (but they 
are not entirely consistent, sometimes using the BBH form). The difference is that the 
operator-equivalent coefficients, (J || a n \\ J), are not included in the crystal-field 
parameters. To avoid confusion with BBH, we use S 1 ™ to denote the A™(r n ) (the 'B™ 



of Scott |70[). Then, by equation A.3 



5™ = (J || a n || J)S™. (A.4) 
The crystal-field Hamiltonian may now be rewritten as 

H*= E E S™(J\\a n \\J)0™. (A.5) 



n=2,4,6m=0...n 
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O 2 = 3 J 2 - J ( J + 1) 

0\ = \[Jz(J+ + J-) + (J+ + J-)Jz] 

Ol = \[Jl + Jl] 

Ol = 35 J 4 - [30 J (J + 1) - 25] J 2 - 6J(J + 1) + 3 J 2 (J + l) 2 

0\ = \[(7Jl - 3J(J + \)J Z - J Z )(J+ + J_) + (J+ + J_)(7J 3 - 3J(J + 1)J 2 - J z )\ 

0\ = \[(7J 2 - J(J + 1) - 5)(J 2 + J 2 ) + (J 2 + J 2 )(7J 2 2 - J(J + 1) - 5)] 

of = + Jl) + (4 + j 3 -)j z ] 
ot = \[Ji + Ji] 

Ol = 231J 2 6 - 105[3J(J + 1) - 7] J 4 + [105J 2 (J + l) 2 - 525J(J + 1) + 294] J 2 2 
-5J 3 (J + l) 3 + 40J 2 (J + l) 2 - 60J(J + 1) 

0\ = ^[(33J 2 5 - (30J(J + 1) - 15) J 2 3 + (5J 2 (J + l) 2 - 10J(J + 1) + 12))J 2 (J + + J_) 
+(J+ + J-)(33J 2 5 - (30J(J + 1) - 15) J 2 3 + (5J 2 (J + l) 2 - 10J(J + 1) + 12)) J z ] 

0\ = ^[(33J 4 - (18J(J + 1) + 123)J 2 + J 2 (J + l) 2 + 10J(J + 1) + 102)(J 2 + J 2 ) 
+(</+ + J^)(33J 2 4 - (18J(J + 1) + 123) J 2 + J 2 (J + l) 2 + WJ(J + 1) + 102)] 

Of = J[(11J 3 - 3J(J + 1) J 2 - 59J 2 )(J 3 + J 3 ) 
+(J 3 + J 3 _)(1U 2 3 - 3J(J + 1) J 2 - 59J 2 )] 

O 4 = i[(HJ 2 - J(J + 1) - 38)(J 4 + Jl) + (J 4 + Jl)(lU 2 - J(J + 1) - 38)] 

o 6 5 = J[j,(4 + j^) + (4 + ^)j,] 

Ol = \[J e + + J 6 ] 

Table A. 4: Stevens operator equivalents, O™ (J± = (J x ± iJ y )). 
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A. 2 Generalization of the Stevens operator equivalents 



All of the early notations used symmetrized operator equivalents with m > 0, to simplify 
the notation. This was permissible in the high-symmetry systems being studied but is 
not in more general situations. Unfortunately there was no standard generalisation of 
the Stevens method; instead several slightly different notations evolved. 

Bleaney |J defined operators 'O^ 1 ' and 'O^ 2 ' with an unusual normalization. Baker 
and Williams || introduced operators @i m and 0j m of a tesseral nature with notes 



on their relation to the O™ of BBH. Prather [58] also defined a similar tesseral set of 



operators and Hutchings M2] notes both the Stevens and tesseral types. Buckmaster 



[14] used two sets of operators of a similar form to those of Stevens but unsymmetrized 



and with different normalizations. Buckmaster's '0[™' are the of BBH where I 
corresponds to n; here we write Buckmaster's 'OJ 71 ' as O™ (for — n < m < n), defined 
as 

O™ = 0° form = 

= \ (6+ m + 6~ m ) for m = 1 . . . n. (A.6) 



The operators O™ are listed in table A. 5. Unfortunately, it is not easy to devise a 



convenient extension of BBH notation with m < using these operators. 

In an attempt to provide a sensible extension to negative m of the Stevens notation, 
Rudowicz [^] defines a set of oparators 'O^', which we shall denote where n *-* k 
and m *-* q. These are defined for —n < m < n such that for m > 0, O™ = O™. In 
terms of the 6% (the 'Of 1 ' of Buckmaster Q): 

O™ = O° n = 6 n form = 

= 0™ = i(0+ m + 0~ m ) for m > 

= ii(0+ m -Q- m ) form<0 (A.7) 

In terms of these the crystal-field Hamiltonian, W c /, becomes 



Hcf= E E (A. 

n=2,4,6 m=—n...n 
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where B™ = B™ for m > 0. Thus this new, tesseral, form can be directly applied to 
parameters using the BBH notation, simply by setting B™ = for m < 0. 

An alternative approach is to maintain the BBH notation but to allow the B™ (except 
for m = 0) to be complex, retaining the condition m > 0. In terms of Buckmaster's 
operators O™ we may write the extended set of 'Stevens' operators as 



0™ = 0° for m = 

= ±6™ for 777 > 

= i(-l) m O™ for 777 < (A.9) 



and Tiqf becomes 

n *f= E E B nO™ (A.10) 
n=2,4,6 m=—n...n 

where the B™ are complex for m > and we define B~ m = (— l) m i?™* to generate 
B™ for 777 < 0. This condition is also true for m = as the B® are real. 



Alternatively, we may rewrite 7i c f as 



Hrf= E 

n=2,4,6 L m=l...n 



B° n d° n + 2 {2Re(i?™)(0™ + 0- m ) + 2iIm( J B™)(0™-0- m )} 

(A.ll 



which, on comparison with equations A.. 6 and A.9, gives 



Hcf= E E Be(TC+ E E 27lm( J B-)(0™-0- m ). (A.12) 

n=2,4,6m=0...n n=2,4,6m=l...n 

This convention has the advantage of reducing to the original BBH form when the B™ 
are real, while permitting extension to more general situations. It is also useful as it 
allows easy conversion to tensor-based notations where the use of complex parameters 
and summation over m = —n ... 77 is natural. 
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6 2 ° = [3 4 - J( J + 1)] 
Of 1 = \[J Z J± + J±J Z ] 

of = 4 

61 = 354 - [30J(J + 1) - 25] 4 - 6J(J + 1) + 3J 2 (J + l) 2 

Of = ~ 3J(J + " + ^±(74 - 3J(J + 1)J 2 - J z )\ 

Of = \\iyJl ~ J(J + 1) " 5)4 + 4(?4 - J(J + 1) - 5)] 

6^ = ^4 + 4^] 
of 4 = 4 

6g = 231J 2 6 - 105[3J(J + 1) - 7]4 + [105J 2 (J + l) 2 - 525J(J + 1) + 294]4 
-5J 3 (J + l) 3 + 40J 2 (J + l) 2 - 60J(J + 1) 

Of = -[(334 - (30J(J + 1) - 15)J 3 + (5J 2 (J + l) 2 - 10J(J + 1) + 12)) J 2 J ± 
+J±(33J 2 5 - (30J(J + 1) - 15)4 + (5J 2 (J + l) 2 - 10J(J + 1) + 12)) J z ] 

0± 2 = ^[(334 - (18J(J + 1) + 123)4 + J 2 (J + !) 2 + 10J(J + 1) + 102)4 
+4(334 - (18J(J + 1) + 123)4 + «/ 2 («/ + I) 2 + 10J(J + 1) + 102)] 

Of = >14 - 3J( J + 1) J, - 59J 2 )4 + 4(H4 - 3J( J + 1) J. - 59J 2) ] 

6 ±4 = ^[(H4 - 4«/ + 1) - 38)4 + 4(H4 - J(J + 1) - 38)] 

6 6 ±5 = 1[j 2 4 + 4j 2 ] 
of = 4 

Table A. 5: The operators of Buckmaster, O™ for n = 2, 4, 6 and m = . . . ± n. 
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A. 3 The notation of Lea, Leask and Wolf 



In cubic symmetry the crystal-field interaction is specified by only two independent pa- 
rameters. It can be expressed most simply if the crystal-field axes (x, y, z) are chosen to 
be along the local fourfold symmetry axes. In which case the crystal-field Hamiltonian 
is 

H cf = B A {Ol + 50f) + Be(O 6 - 210|) (A.13) 
where the O™ are Stevens operator equivalents in the notation of BBH. From equation 



A.13 it is clear that the B^ and Bq are identical with the B® and B® of BBH. Lea, Leask 
and Wolf [48] rewrite TL c f as 



H cf = 5 4 F(4)0 4 + B 6 F(6)0 6 (A.14) 
where 

_ 02 + 5Q| Q 6 ° + 21Q| 

4 = ir( 4 ) and ° 6 = ir (6 ) • ^ A - 15 ) 

The scaling factors, -F(4) and F(6), were chosen to keep the eigenvalues of the fourth 
and sixth-order operators in the same numerical range. This significantly assisted 
computation, given the limited computing resources available at the time. Values of 
F(4) and F(6) for different J-manifolds are tabulated by Lea, Leask and Wolf |48T| . To 
separate the strength of the crystal field from the relative magnitudes of the fourth and 
sixth-order terms, they write 

B A F(4) = Wx and B 6 F(6) = W(l - \ x |) (A.16) 

where — 1 < x < +1. In terms of the new parameters, W and x, the Hamiltonian 
becomes 

H cf = W[xO i + (l-\x\)0 6 }. (A.17) 

This expression is valid only when the crystal-field axes coincide with the fourfold 
symmetry axes. However, it can readily be transformed to a system in which the z axis 
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is along a twofold or threefold symmtery axis. This notation has been extensively used 
for crystal-field calculations in cubic symmetry. 



Walter 76] extends the notation of Lea, Leask and Wolf to systems with lower sym- 
metry. Like the cubic parameters W and x, a magnitude parameter, W', is used in 
conjunction with a set of bounded weighting factors x' nm where | x' nm | < 1. It should 
be noted that the normalisation factors used by Walter in the cubic case are slightly 
different from those of Lea, Leask and Wolf. This notation may be of use for fairly 
high symmetry, such as C$h hexagonal symmetry where four parameters are required. 
However, for low symmetry, it becomes cumbersome. 



A. 4 Notations used by optical spectroscopists 



Morrison and Leavitt [52] summarise what is probably the most common notation 
among optical spectroscopists. They use a set of tensor operators C q (sometimes writ- 
ten l Cq ) with coefficients l B q \ referred to here as M k to avoid confusion with BBH. 
The crystal-field Hamiltonian takes the form 

n cf = £ £ M k q C k q (A.18) 

k=2,ifiq=-k...k 

where C q = (-l) q (C k )*. For H c f to be real, this implies that 

M k q = (-l) q (M k )*. (A.19) 

Thus it is necessary to specify the values of M k only for q = ... k. This condition 
also implies that the Mg must be real. The crystal-field parameters M k are common 
to manifolds of different J which is essential when interpreting optical data. It is the 
inclusion of operator-equivalent coefficients or reduced matrix elements within param- 
eters such as those of BBH that limits their application to a single J-manifold, see the 
conversions in section |A.6| . 

Wybourne [j?^] uses a scheme based on the same tensor operators, C k , but re- written 
to give a set of real parameters. The parameters are denoted l B g k ^ by Wybourne 77] 
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but here we will use Wj? to distinguish them from other notations. The crystal- field 
Hamiltonian becomes 



fe=2,4,6 q=-k...k 



(A.20) 



where the 2^ are defined in a similar way to the tesseral harmonics (see, for example, 
Hutchings p3]]), as shown below. 



(C* +(-l)«C; 



-i(C* - (-im fe ) 



for g = 
for g > 
for g < 



(A.21) 



Wybourne's notation is widely used, see for example |34|, 23]. Conversion between the 
parameters Mg of equation A.1S and the Wg of equation A.2C is simply a matter of 



multiplying some by —1, see section A.6 



A. 5 Computation using the 3-j symbols 



This section presents a notation based on that of Edmonds [ 30 1 , which provides a 
convenient intermediate for computation. However, we do not suggest that this form 
should be adopted in addition to the many others already used; that would only add 
to the current confusion. The notation of Morrison and Leavitt [^] is probably the 
most convenient for general use, including low symmetry systems. 



Let us define the spherical harmonic, Yj iq (9,ip) (Edmonds |30| equation 2.5.29, or Arf- 
ken as 



Y kq {e^) = {-iy 



(2k + l) (k-q)l 
47r (k + q)\ 



P fc 9 (cos 0)e 



(A.22) 



Expressions for associated Legendre functions Pf} (cos 9) are given by Arfken ||. The 
spherical harmonics Yk„{d,<p) are orthonormal over both 9 and <p. From these we may 
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define the solid harmonics, ([30] equation 5.1.2), as 



y k q {r) = r k Y kq (6^). 



(A.23) 



In the notation of spherical tensors, the k th rank tensorial set T k may be built from 
the (2q + 1) functions y k for q = —k ... A;. In particular, the tensorial sets produced 
are irreducible. This means that the irreducible tensor operators T k transform under 
rotation (r-transformation) in subsets T k (see appendix B). By application of the 
Wigner-Eckart theorem we may factorize the matrix elements of T k : 



( 7 JM | T„ fc | yj'M') = (-1 



\J-M 



J k J' 
-M q M' 



( 7 J || T k || 7 'J A 



(A.24) 



where index 7 denotes all other quantum numbers required to specify the J-manifold 

/ . . . 

31 32 33 



The 3-j symbols 



are defined by Racah |3S|, 30]. In this work we 

y mi ni2 ' 

consider the crystal- field interaction only within a J-manifold, so 7 = 7' and J = J' . 



Thus, equation A.24 simplifies to 



(JM I T fc I JM') 



•1) 



J-M 



J k J 



y -M q 



M' 



(J || T k || J), 



(A.25) 



where the redundant index 7 has been omitted. We may now express the crystal field 
interaction in terms of a set of parameters E k and the operators T k : 



n 



fc=2,4,6 q=-k...k 



(A.26) 



The matrix elements of 7i c f are 



(JM I 7i cf I JM') ^ E k {JM I T g fc I JM') 

fc=2,4,6g=-fe...fc 



(A.27) 



Substituting the expression for the matrix elements of T q (equation A.25| ) gives 



(jm 1 n C f I jm') = 53 53 £j(-i) J - M 

fc=2,4,6g=-fc...fc 



J k J 
-M a M' 



(J || T k || J). 
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(A.28) 



This form is useful for computation, because an explicit formula for the 3-j symbol 
is given by Edmonds |30| (equations 3.7.3 and 3.6.11). The reduced matrix elements 
(J || T k || J) are given by 



(J II T fc || J) = 1 



(2J + £; + !)! 



(2J-k)\ 



(A.29) 

A similar form for the crystal-field parameters has been used by Bishton and New- 
man H in the study of electron correlation effects. The tensor operators, Tjf, are listed 
by Smith and Tornley [73] for k = ... 6 and q = 0. . . ± k (labelled 'Og ) and by 



Buckmaster JTJ] for k = 2,4, 6 and q = . . . ± k (labelled 'Ofc± 9 ')- The operators for 
k = 2, 4, 6 and q = . . . ± k are reproduced in table [A.6|. These are not the same as 



the tensors given later by Buckmaster [15] 



The tensors operators, , used here are of the same form as the used before, but 
without the host-dependent normalization. The relationship is simply 

C\ = {J || a k || J)T*. (A.30) 

As noted earlier, the operator-equivalent coefficients (J || \\ J) depend on both the 
ion and the coupling scheme (L-S or intermediate). Thus these must be known to 
convert between notations. 
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To =\{ZJl-, J {J + l)\ 
Tl l =T^-[JzJ± + J±Jz\ 



J ±2 - —r J ± 



T 4 = -[354 - [30 J (J + 1) - 25] 4 - 6J(J + 1) + 3 J 2 (J + l) 2 ] 

T±i = =F^[(7J| - 3J(J + l)J g - J Z )J± + J±(7J 3 Z - 3J(J + l)J z - J z )\ 



ri, - - J{J + 1) - 5)4 + 4(7 4 - J(J + 1) - 5)] 



7±3 — T— 3— [^z4 + 4^2 



T 4 _ V70 4 



T 6 = ^[ 231 4 - 105[3J(J + 1) - 7] J 4 + [105 J 2 (J + l) 2 - 525 J( J + 1) + 294] J 2 
-5J 3 (J + l) 3 + 40J 2 (J + l) 2 - 60J(J + 1)] 



^1 = =F^[(33J| - (30J(J + 1) - 15) 4 + (5J 2 (J + l) 2 - WJ(J + 1) + 12)) J Z J± 

+J±(33J| - (30J(J + 1) - 15) J 3 + (5J 2 (J + l) 2 - WJ(J + 1) + 12)) J z ] 
\/105 

T ±2 = ^j-[( 3 34 " ( 18J (^ + 1) + 123) Jf + J 2 (J + l) 2 + WJ{J + 1) + 102)4 
+4(33J Z 4 - (18J(J + 1) + 123)4 + J 2 {J + !) 2 + 10J(J + 1) + 102)] 

7±3 = ^^-[( n 4 - 3J ( J + ^ - 59J*)4 + 4(H4 - 3J(J + - 59J,)] 
3±4 = ^ [(114 - ■/(■/ + 1) - 38)4 + 4(H4 - 4^ + 1) - 38)] 

T i5 = ^[J,4 + 44 
6 _ V231 6 

J±6 " "32"^ 



Table A. 6: Irreducible tensor operators for g = 2, 4, 6 and q = 0...±k. 
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A. 6 Relationships between crystal-field conventions 



This section describes the relationships between the notations described in this ap- 
pendix. The task of finding parameters for one convention from those of another is 
made unnecessarily complicated by the fact that most authors use the same set of sym- 
bols (B?) to denote quite different parameters. We use the indices k and q throughout, 
where k «-> n and q «-> m as appropriate, and the following names denote the different 
conventions: 



BBH for the operator equivalent convention of Baker, Bleaney and Hayes ||]. In the 
present work Bf, denotes the 'B™ of the Baker, Bleaney and Hayes. The B\ are 



defined in equation |A.2| and extended to complex parameters in equation A. 11 
This form is also used by Abragam and Bleaney All in some places. 



SCOTT for the convention of Scott jT0| and Abragam and Bleaney |lj (not used con- 
sistently by Abragam and Bleaney). In the present work S\ denotes the l B™' of 
Scott and the i A^(r n y of Abragam and Bleaney. The are defined by equation 



A.5 



MORRISON for the convention of Morrison and Leavitt f53[|. In the present work 



denotes the l B^ of Morrison and Leavitt. The are defined by equation A. 18 



WYBOURNE for the convention of Wybourne [0] . In the present work Wt denotes the 
'B*' of the original work. The Wt are defined by equation [A.20| . 



EDMONDS for the notation based on the 3-j symbols which provides a useful intermediate 
for calculation. The Eg are defined by equation A. 28. An explicit formula for the 



3-j symbols is given by Edmonds [ 30 1 and avoids the need for tables of operators 
or matrix elements. 



The MORRISON and WYBOURNE notations are very similar as they are based on the same 
tensor operators, C*. However, the MORRISON parameters, are real for q = and 
complex for q = 1 . . . k whereas the WYBOURNE parameters, Wjf are real for q = —k...k. 



152 



MORRISON parameters are such that 
M fc real, k = 2, 4, 6 
Mq complex for q = 1 . . . k, k = 2, 4, 6 
M* = (_l)*(M* )*. 



WYBOURNE parameters are such that 

W~ real for k = 2, 4, 6 and q = —k . . . k. 



Conversion between MORRISON and WYBOURNE parameters simplifies to: 



Re(M*) 
Im(M g fc ) 



(-l)(9 + 1 )W^g 



for q = 1 . . . k 
for g = 1 ... fc. 



(A.31) 



Conversions between BBH, EDMONDS, SCOTT and MORRISON notations are a simple matter 



of multiplying factors as shown below, and diagramatically in figure A.l 



B\ = (J || a k || J)S* 



(A.32) 



fk\q\^k 



(A.33) 



^ = II «fc II 



(A.34) 



The /j.ui of equation |A.33| are listed in table [A.7| . They are also tabulated by Kass- 
man |^6| and by Buckmaster [14|; Wybourne's [77] tabulation contains many errors. 
The (J || «fc || J) depend on the ion and the coupling scheme assumed (eg. L-S or 
intermediate). An extensive tabulation of (J || a k \\ J) for the rare earths is given by 
Han fll 
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WYBOURNE 



(equations A. 31) 



SCOTT 



MORRISON 



fk\q\ 



{J || a k || J) 



(J II «fc II J) 



BBH 



EDMONDS 



lk\q\ 



Arrow direction corresponds to multiplication by the factor given. Conversions involving 
(J || a k || J) depend on the ion, the coupling scheme and J. 



Figure A.l: Relationship of crystal-field conventions and conversions between them. 



f: 



20 



J21 — 7T 

D 



122 



/40 
hi 
/42 
/43 



2VE 



5 

2yTo 
5 

2^/35 
35~~ 

4^70 
35 



/60 
/6i 

/62 
/63 
/64 
/65 
/66 



16 

4^/42 
21 - 

16\/i05 
105 

8\/l05 
105 

8-/I4 



21 

8jy/77 
23T 

16\/231 
231 



Table A. 7: Factors fu q i to convert crystal-field parameters from SCOTT to MORRISON 
notation. 
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B. Coordinate rotations for 
crystal fields 



This appendix describes the transformation of crystal field parameters under a general 
rotation of the coordinate system. For this we need to understand a little about the ten- 
sorial sets used to parametrize the crystal field. In what follows we use the language of 



Fano and Racah [ 36 1 . Tensorial sets are defined in such a way that a coordinate rotation 
corresponds to a linear transformation of the components of each tensor. Transforma- 
tions corresponding to coordinate rotations are called r-transformations. We may write 
a tensor as 



T = Y^EiTi (B.l) 

i 

where the Tj are unit tensors and the Ei are the coefficients or tensor components. 
In particular, irreducible tensorial sets are minimal subsets of tensors that transform 
without admixture from other subsets. If the index k identifies individual irreducible 
sets and q denotes the elements (tensors) of each set, we may rewrite T as 

T = J2T k = J2 E " T "- ( B - 2 ) 

k k,q 
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In three-dimensional cartesian space the spherical harmonics form irreducible sets. In 
quantum mechanics, irreducible tensor operators may be defined by the following com- 
mutation rules: 



J z ,T k = qT k and (B.3) 



J ± ,T k = J(k±q + l)(k^q)T k (B.4) 



Each fc-th degree^ tensor T fc has order (2k + 1). That is, T k has (2k + 1) components, 
the operators T k for q = —k ... ... A;. 

T is invariant under coordinate transformation; however both the components and the 
unit tensors depend on orientation. Let D k be the unitary matrix corresponding to a 
coordinate rotation of T k . The components in the rotated coordinate system, E' k are 
given by 

q= — k 

The summation is over q = —k . . . k; however to represent a real crystal field poten- 
tial we have the condition E k = (—l) q (E^_ )*. This condition is invariant under the 
transformation so there is no need to calculate E'J? for q' < 0. If, as is usual, the 
new coordinate system is re-labelled x, y and z then there is no need to transform the 
operators T k 



Equation B.2 is the same as equation A. 26 if T is identified with Tl c f- Equation B.2 



is also similar to equation A. 18 but the normalisation within each irreducible set is 
different. Thus, to perform a coordinate rotation of a set of crystal field parameters (in 
either MORRISON or EDMONDS form) we simply require the matrices D k of equation |B.5| . 
For application to crystal-field models with 4/ electrons we need consider only the cases 
k = 2, 4 and 6. 

A general coordinate rotation may be described by an angle £ and an axis u. However, 
a more common notation specifies a general rotation by three Euler angles ip, 9 and (p. 
The rotations are performed in sequence: ip about z, then 9 about y' and then if) about 



1 Many authors use the term rank instead of degree. 
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Figure B.l: Schematic of the general rotation from (x, y, z) to (x'" , y" 1 , z'") coordinates, 
specified by the Euler angles </?, 9 and ip. 



z" (see figure B.l). The corresponding standard r-transformation of J is then 



D^, 6, if) = e^'e^e^'. (B.6) 

If the axis z is chosen such that {JM' \ J z \ JM) = MSm'Mi then the transformation 
reduces to 

Di, q (^e, ( p) = e i ^'(e idJ y) q , q e i ^. (B.7) 

The i/j and (p rotations only involve multiplication of the parameters by phase factors 
and may calculated simply in most of the common crystal-field notations. However, the 
6 rotation will, in general, mix all possible terms within a given irreducible tensorial 
set. In the language of crystal fields, all parameters Eg for a given k may be admixed. 
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Even in high symmetry systems where the initial choice of axes results in few non-zero 
parameters, parameters for all < q < k may be generated if the crystal field does not 
have high symmetry about the new axes. Let us define the notation 



(Bi 



and hence rewrite equation B.7 as 



% q ^,e^)=e^'d k qlq {9)e^. 



(B.9) 



Fano and Racah (3^, appendix D] derive an expression for the matrix elements of d k : 

^{k + q')\{k-q')\{k + q)\{k-q)\ 



d 9 \(0) 



k+q 

£(-!) 

r=0 



(k — q' — r)\{k + q — r)!r!(r + q' — q)\ 



2 2 



(B.10) 



where terms involving factorials of negative numbers are omitted. This is equivalent 
to extending the defintion of n\ such that 1/n! = for all negative integral n 
Most texts are unclear on this point. Fano and Racah note some useful symmety 
properties: 



4«w 
4«w 



(-i) q '- q d k qgl (9) 

(-iy'-i6 k L q ,_ q {6) 
(-l) k -*'d k ,_ q (n-9) 



(B.ll) 
(B.12) 
(B.13) 



Equations B.£ and B.10 may be combined to given a single expression for a general 
rotation, 



D^(V,M 



D (q' ip+qip) 



k+q 



E(-D 



r=0 



. ^{k + q')\{k-q')\{k + q)\{k^qy. 
(k — q' — r)\(k + q — r)!r!(r + q' — q)\ 



x( C0 ^)( 2fe -<?'+9-2r)( sin |)(2r +g '-,)_ 



(B.14) 



As noted earlier, rotation about the z axis is straightforward and equation [B.14| sim- 
plifies to 



D k q , q (0,0,<p)=e iqv 8 q , g , 



(B.15) 
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as expected. In this case the matrices D k are diagonal so there is no mixing of compo- 
nents with different q. Thus crystal-field parameters of the conventions denoted BBH, 
SCOTT, MORRISON and EDMONDS will all transform correctly. However, the BBH and SCOTT 
conventions use operators which do not form tensorial sets and hence parameters of 
these conventions do not undergo general rotation with D k (ip, 6, ip). Parameters using 
tesseral harmonics (e.g. WYBOURNE) do not transform properly under either rotation. 

Coordinate rotations are discussed in many books on angular momentum but there is 



wide variation in notation. A good discussion is given by Rose [66, appendix II], giving 



in the same form as equation B.14 . A particularly succinct introduction is given 



by Brink and Satchler [12|. 
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C. Crystal- field software 



Two programs have been written in collaboration with R G Graham of this department. 
I have written a program to convert crystal-field parameters between conventions, to 
change the units in which they are expressed, and to allow rotation of the coordi- 
nate system. This program is based on the work described in appendices A and B. 
R G Graham has written a program to build and diagonalise the combined electronic- 
nuclear Hamiltonian for a single J-manifold of any rare-earth ion in a site of arbitrary 
symmetry. This appendix briefly describes these programs. 

The software has been tested by comparing the output with hand-calculated special 
cases and with the output of earlier programs for higher symmetry. We have also used 
it to reproduce published results from the crystal-field calculations of other workers. 

C.l Parameter conversion program: CFPCONV 

The CFPCONV program allows the crystal-field and other parameters of the Hamiltonian 
to be specified in a 'human-readable' form. Many constants are stored in CFPCONV so 
they only need to be specified if non-standard values are required. For example, by 
specifying the ion, CFPCONV assumes the latest published values for A, C, I, J, g n 
g 3 and the (J || || J) for the gound manifold in the intermediate coupling approx- 
imation. Crystal-field parameters can be converted between the notations discussed 
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in appendix A. After manipulation, the parameters may be saved either in the same 
format as they were input, or fomatted for the REION program (using the EDMONDS no- 
tation). In this way CFPCONV provides an interface to REION, the program to diagonalise 
the Hamiltonian and calculate the expectation values of operators of interest. 

The SCOTT, BBH, WYBOURNE and MORRISON notations are as given in appendix A. The 
intermediate form denoted EDMONDS is used by CFPCONV to create input files for REION. 
However, the normalisation is slightly different to that given earlier: 

E k q =E k g (J\\T k \\J), (C.l) 

where the E q are as given in appendix A. In terms of the parameters E q , the matrix 
elements of the crystal-field Hamiltonian are, 



(jm | n cf | jm') = E E K(~ l ) J ~ M I . I • (C2) 

fe=2,4,6 q=-k...k 



The parameters E are used in order to reduce the calculation to its simplest form: 




the multiplication of 3-j symbols by the parameters. Equation (12 can be straightfor- 
wardly implemented because all constants depending on the ion and coupling scheme 
are included in the parameters. 

General rotation of the coordinate system is achieved by repeated rotation about the 
z and y axes (rotations D((p, 0,0) = D(0, 0, <p) and D(0,#,0) in the notation of ap- 
pendix B). Rotation about the y axis is not supported in the SCOTT, BBH and WYBOURNE 
notations. In these cases the parameters must first be converted into MORRISON or 
EDMONDS notation, and then back again after transformation if desired. 



161 



C.l.l Example input file for CFPCONV 



Figure GA shows an example input file for the CFPCONV program. The crystal-field 



parameters are those of Sharma et al [[nj for H0F3, corrected by Ram and Sharma [65]. 
The first line is the title; comments can be appended to statement lines or added after 
the END statement. Crystal-field parameters not specified (B21 for example) are taken 
as zero by default. The parameters are in the WYBOURNE notation with units of cm -1 . 



#CPFs for HoF3 from [SHARMA+81] , set la (B6-4 from [RAM+85] ) 
FORMAT WYBOURNE 
UNITS cm-1 
ION Ho 



QEXT 


9 


66E- 


-4 




TEMP 


4 


2 






B20 = 


78 


5 






B22 = 


-199 


6 & 


272 


7 


B40 = 


17 


9 






B42 = 


-101 


7 & 


-33 


4 


B44 = 


265 


8 & 


-95 


9 


B60 = 


224 


6 






B62 = 


156 


6 & 


-552 


3 


B64 = 


58 


1 & 


-48 


12 


B66 = 


-405 


9 & 


31 


3 



END 



Figure C.l: Example input file for CFPCONV. 
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C.2 Calculation program: REION 



This program is written in Fortran77 on an IBM PC compatible using the Salford 
FTN77/386 compiler, which allows the use of more than 640 kBytes of contiguous 
memory (the usual limit under MSDOS). Taking the most complex example, the com- 
bined electronic-nuclear Hamiltonian for the 5 Ig ground manifold of 165 Ho (J = 8, 
I = |) requires a 136 by 136 complex matrix. Using double precision (64 bit) numbers, 
the Hamiltonian matrix occupies over 300 kBytes of memory. In practice several matri- 
ces of this size are used to calculate the complete Hamiltonian, so more than 1 MByte 
is required. 

The terms included in the Hamiltonian are described in chapter 2. Terms may be 
excluded from the calculation by setting the appropriate parameter to zero. Similarly, 
a purely electronic calculation is selected by setting / to zero; the calculation is, of 
course, much faster with the smaller Hamiltonian. REION calculates the energy levels, 
eigenstates, transition moments and expectation values of the electronic angular mo- 
mentum operators. Thermal averages are calculated for a temperature specified by the 
user; ground state values can be selected by specifying zero temperature. 

On an IBM PC compatible computer (33 MHz '486DX) the building and diagonalisation 
of a 136 by 136 Hamiltonian matrix takes ~ 5 min. We have also modified a version to 
run on Unix workstations. When running on a Hewlet Packard 720 workstation, this 
version is several times faster than the PC version. 



C.2.1 Example input file for REION 

Figure |C.2j shows an annotated example input file for the REION program. This was 



obtained by manipulation of the CFPCONV input file for H0F3 shown in figure C.l 



By rotating the coordinate system, the crystal-field z axis has been chosen along the 
crystallographic a axis, rather than the C\h axis (the b axis, which is the crystal-field 
z axis for the original data). This choice results in non-zero values for all the E which 
are specified as real part then imaginary part. All energies are in kelvin; angles in 
degrees; and fields in tesla. 



163 



251 

temp . out 
temp . swp 

CPFs for HoF3 from [SHARMA+81] 

1000.00 

3.50000 

8.00000 

1.17900 

1.24150 

0.03897430 

2.50998E-05 

0.00000 0.00000 0.00000 
0.00000 

0.00000 0.00000 0.00000 

0.00000 0.00000 0.00000 

0.00000 

0.00000 

9.66000E-04 

4.20000 

-177.873 0.00000 

0.00000 236.374 
-128.173 0.00000 
-184.913 0.00000 

0.00000 91.8854 
-198.755 0.00000 

0.00000 2.41033 
-99.4329 0.00000 
-397.949 0.00000 

0.00000 -233.904 
-1020.99 0.00000 

0.00000 817.491 
-101.958 0.00000 

0.00000 -940.023 
-57.5958 0.00000 
END 



Output control byte 
Output file name 
Secondary output file name 

set la (B6-4 from [RAM+85] ) 
Maximum energy for output 
/, nuclear angular momentum 
J, electronic angular momentum 
g T , nuclear ^-factor 
g d , electronic ^-factor 
A, dipolar hypferfine constant 
C, quadrupolar hyperfine constant 
Magnitude of applied field: initial, final, step. 
Magnitude of dipolar field 
9 a for the applied field: initial, final, step. 
(p a for the applied field: initial, final, step. 
9(1 for the dipolar field 
(fd for the dipolar field 

Qext 

T, temperature 
E (real) 

E (real) 
El 

E (real) 

% 



Figure C.2: Annotated example input file for REIQN. 
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